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effects  be  Investigated  at  lower  pressures  and  temperatures.  Some  pressure 
and  temperature  conditions  exist,  probably  Important  to  low  temperature 
plumes,  below  which  correlation  effects  In  the  4.3  y m C02  band  could 
become  Important, 
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mental data.  We  realized  that  these  correlation  effects  were  important  to  the  prob- 
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NOMENCLATURE 
(units  are  given  in  parentheses) 
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a 

b 

B 
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eu 


E 

Eq 

f(So,E),f(So) 


g 


g5 


h(S0,E,Az7,y) 

I0,I1 


K 

k 

k(a) 

L~,L 
v ’ 

1. 

1 

m 


m. 

N 


pe’Pe. 


Q 

Q 


path  coordinate,  amount  of  absorber  per  unit  cross  section 
(atm-cm,  or  molecules/cm^) 


parameters  for  function  in  Table  1 


beta  function  (see  Reference  11) 
velocity  of  light  (3  * 10*®  cm/sec) 

extinction  coefficient  at  wavenumber  v in  inverse  path  coordinate  units 
(atm-cm)”1  or  (molecules/cm^)”1 

lower  state  energy  (cm  *) 

lower  state  energy  at  SQ  = Se  (cm’*) 

joint  distribution  function  for  line  strength  and  lower  state  energy,  and  dis- 
tribution function  for  line  strength  (in  inverse  strength  units) 


l/S  jj—  (dimensionless) 

degeneracy,  Equation  21  only 

-27 

Planck’s  constant  (6.62  * 10  erg  sec) 

product  of  line  shape  function  and  S/S,  Equations  17  and  18  only 

Bessol  function  of  the  first  kind  with  imaginary  argument,  zeroth  and  first 
order 

rotational  quantum  number 


source  function  in  equation  of  transfer  (in  radiance  units) 
constant  for  distribution  function  (dimensionless) 

Boltzmann’s  constant,  Equations  13,  21,  24,  44  only  (1.38  x 10  ergs/K) 

_ 2 -1 

average  line  strength-to-spacing  ratio,  S/d  (atm-cm)  or  (molecules/cm  ) 

9 -1 

spectral  radiance  (W/cm''-sr-cm  ) 

2 -1 

spectral  radiance  of  blackbody  at  temperature  of  substance  k (W/cra  -sr-cm  ) 
th 

physical  length  of  i 11  element  of  path  of  n elements  (cm) 

temperature  exponential  for  strength  variation  with  temperature  (dimensionless) 
rotational  m-factor,  -j  for  P-branch  transitions,  j+1  for  R-branch  transitions 
number  of  lines  in  resolution  interval 

iu 

effective  absorbing  pressure  (atm),  of  i species  when  subscripted 

til 

local  and  average  effective  broadening  pressure  (atm)  of  r species  when 
subscripted 

tTt  / 

parameter,  2uw  /(I  + u) 

partition  function,  Equation  21  only  (dimensionless) 


II 
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q,q  local  and  path-averaged  value  respectively  of  the  temperature  ratio  T /T 

(dimensionless)  J 

Rjj  matrix  element  of  a generalized  transition  (dimensionless) 

S absorption  line  strength  (cm*  Vatm-cm) 

S line  strength  averaged  over  all  lines  in  a spectral  resolution  interval 

(cm“  Vato-cm) 

S iine  strength  value  (at  T0)  at  which  assumed  dependence  of  Sn  upon  E goes  to 

E = 

° -1 
SQ  Hne  strength  at  temperature  Tq  (cm  /atm-cm) 

T temperature  (K) 

Tq  reference  temperature  (296K) 

(T0/F)  path-averaged  value  of  the  temperature  ratio,  T0/T  (dimensionless) 

u dimensionless  optical  depth,  defined  in  Equation  29 

v parameter,  Ap/y 

equivalent  width  of  a line  (cm*^) 

average  equivalent  width  averaged  over  all  lines  in  a resolution  interval 
(cm'l) 

w parameter  S/S 

o e 

x another  form  of  dimensionless  optical  depth  (see  first  footnote  in  Section  2.7) 

y,y  local  and  average  value  of  second  band  model  parameter  defined  in  Equation 

29  (dimensionless) 

yQ  value  of  second  bandmodel  parameter  at  STP  (dimensionless) 

«aj  broadening  coefficient  for  species  j with  active  species  a (dimensionless) 

<x»,a  cell  window  and  sample  gas  absorptances  respectively  (dimensionless) 

© 

P second  dimensionless  band  model  parameter  in  another  form  (see  first  footnote 

in  Section  2.7) 

y,y  local  and  average  values  of  the  pressure-broadened  halfwidths  (cm  *) 

y * normalized  incomplete  gamma  function  (see  p.  260  of  Reference  11) 

T gamma  function 

~ -1 
v wive  number  (cm  ) 

~ -1 

At;  resolution  interval  (cm  ) throughout  most  of  the  report.  Also  used  as  a 

dummy  variable  of  integration  describing  spectral  distance  from  line  cent;  - 
in  Equations  12  through  30 

p reflectance  of  a single  solid-air  interface  of  the  cell  window  material 

(dimensionless) 

r transmittance  (dimensionless) 

t / transmittance  of  the  sample  in  the  white  cell,  and  of  the  empty  white  ceil 

g /c  respectively  (dimensionless) 
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SUMMARY 

An  experimental  and  theoretical  study  has  been  made  of  the  atmospheric  attenuation  of 
plume  radiation.  A particular  facet  of  this  problem  was  examined,  namely  the  effect  of  the 
correlation  of  the  water  vapor  and  C02  emission  lines  of  the  plume  with  the  corresponding 
absorption  lines  in  the  atmosphere.  These  correlations  can  produce  lower  effective  transmit- 
tance of  the  plume  radiation  than  the  usual  treatments  of  atmospheric  attenuation  predict. 

To  measure  the  magnitude  of  correlation  effects,  experimental  measurements  were  made 
of  the  radiance  of  a hot  cell,  filled  with  mixtures  of  C02,  H20,  N2  and  other  plume  gases. 

These  radiances  were  measured  both  unattenuated  and  attenuated  by  a long  cool  atmospheric 
path.  A band  model  method  to  handle  the  line  correlation  effects  was  also  developed. 

We  were  able  to  conclude  that  these  correlation  effects  are  severe  enough  In  the  2.7  u m 

water  vapor-C02  band  that  their  neglect  can  lead  to  an  error  of  about  50?,  to  a factor  of  two 

in  computed  radiances  and  radiant  Intensities  for  horizontal  paths  at  15  km  altitude.  Correl- 
ation effects  were  minor  in  the  4.3  p m band  except  In  the  region  of  the  band  head  where  the 

neglect  of  correlation  leads  to  about  a 20X  error  in  a 10  km  path,  and  in  the  band  center 

where  there  was  too  much  absorption  to  obtain  a good  signal -to-noise  ratio.  iTie  measurements 
did  not  extend  to  several  higher  simulated  altitudes.  Therefore,  it  was  not  possible  to  make 
a quantitative  statement  about  the  error  expected.  If  correlation  effects  were  neglected,  for 
more  realistic  situations,  such  as  slant  paths  to  space.  However,  correlation  effects  are 
expected  to  become  more  Important  at  lower  pressures  and  temperatures. 

We  were  also  able  to  show  that  the  set  of  emission-absorption  parameters  compiled  by 
Ludwig  ^ and  the  most  commonly  used  to  compute  C02  contributions  to  plume  radiances, 
do  not  predict  the  measured  emission-absorption  spectra  very  well,  while  the  water  predic- 
tions showed  adequate  agreement.  In  view  of  the  heavy  dependence  of  the  C02  parameters 
of  Reference  1 on  theory,  this  result  is  not  surprising;  it  Is  recommended  that  a program  be 
conducted  to  Improve  the  knowledge  of  the  emission-absorption  parameters  of  C02  so  that 
accurate  plume  radiant  Intensities  can  be  obtained. 

The  measured  emission-absorption  spectra  were  also  compared  to  the  extensive  line  tab- 
ulation of  McClatchey.  The  comparison  was  unsatisfactory  In  most  cases  because  the  line 
tabulation  does  not  include  many  of  the  emission-abosrption  lines  which  are  important  at  high 
temperatures.  Because  of  these  missing  lines,  the  line  tabulation  in  its  present  configuration 
cannot  be  expected  to  accurately  predict  the  emission  of  hot  gases  at  temperatures  over  800- 
1000°K. 


Ti  C.B.  Ludwig,  et  al.,  Handbook  of  Infrared  Radiation  from  Combustion  Gases,  NASA-SP- 

3080,  Scientific  and  Technical  Information  Office,  NASA,  Washington,  District  of  Columbia, 
1973. 
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Simple  band  model  methods  were  shown  to  yield  incorrect  results  when  substantial  line 
correlation  effects  were  present.  However  a band  model  method  was  developed  which  treats 
these  correlation  effects  in  a realistic  manner.  When  compared  with  a line-by-line  calcula- 
tion for  a given  realistic  random  array  of  lines,  this  new  band  model  does  a much  better  job 
of  predicting  the  correlation  effect  than  any  of  the  previously  available  methods.  Its  agree- 
ment with  the  measured  data  is  obscured  to  some  extent  by  the  disagreement  between  the 
available  emission -absorption  parameters  and  the  measured  spectra.  But  where  good  agree- 
ment does  occur,  the  new  band  model  predicts  the  line  correlation  effects  and  hence  the  effect 
of  atmospheric  attenuation  better  than  other  available  methods.  The  use  of  this  band  model, 
in  the  form  in  which  it  is  presented,  is  not  substantially  more  complicated  than  the  use  of 
others;  accordingly  we  recommend  this  band  model  be  further  Investigated  before  the  model 
Is  applied  to  routine  plume  calculations. 

Many  important  aspects  of  this  problem  have  yet  to  be  resolved;  hence  we  recommend 
that  the  following  tasks  be  undertaken. 

(1)  The  numerical  evaluation  of  the  functions  in  the  model  was  not  completed  in  this  pro- 
gram. It  is  important  that  this  be  done  before  widespread  use  is  made  of  the  model. 

(2)  A better  set  of  emission -absorption  parameters  is  needed,  particularly  for  COg. 

The  need  is  urgent  because  there  are  no  existing  parameters  which  accurately  treat 
the  emission  processes  of  the  ’’red  spike"  region  of  the  4.3  /im  C02  band.  The  ex- 
perimental apparatus  used  for  the  present  program  is  ideal  for  the  needed 
measurements. 

(3)  Correlation  effects  in  the  4.3  jim  band  will  become  more  important  as  the  plume 
temperature  is  reduced  from  1200°K  and  the  altitude  increased  above  15  km.  The 
experimental  conditions  used  in  this  program  did  not  treat  these  cases.  Hence,  line 
correlation  effects  in  this  band  should  be  investigated  at  lower  temperatures  (and 
pressures  where  possible)  in  a manner  similar  to  that  done  for  the  program. 

(4)  High-altitude  broadening  processes  were  not  treated  in  the  development  of  the  new 
band  model.  The  band  model  development  should  be  extended  to  include  these  Doppler 
broadening  processes  so  the  techniques  can  be  useful  at  high  altitudes. 


INTRODUCTION  j 

Hot  gas  emission  is  the  primary  source  of  infrared  radiation  from  aircraft  and  rocket  j 

plumes.  The  vibration -rotation  bands  of  the  infrared -active  species  in  the  plume  (primarily  j 

H20  and  C02)  are  the  most  prominent  features  of  the  plume  spectrum.  The  development  of  1 

means  to  accurately  compute  plume  emission  and  irradlance  at  a sensor  are  important  both  5 

for  sensor  effects  and  radiation  suppression  effects.  In  both  these  areas,  the  theory  of  radia-  j 

tive  transfer  provides  the  link  between  the  chemical  and  aerodynamic  processes  governing  j 

the  plume  physical  properties  and  the  observed  emission.  j 

I 

i 

1.1  BACKGROUND  j 

The  emission  processes  from  the  gaseous  species  involved  have  been  studied  in  many  j 

contexts  in  the  past.  The  spectral  details  of  molecular  emission  (and  absorption)  at  very  J 

high  spectral  resolution  have  been  used  to  study  molecular  structure.  Meteorologists  have  | 

long  known  that  the  vibration -rotation  bands  of  HgO  and  C02  molecules  play  an  important 
role  in  the  environment  of  our  planet.  Radiative  emission  from  hot  water  vapor  and  CC>2  has 
also  been  important  in  combustion  heat  transfer  studies.  ^ 

The  study  of  plume  emission  looks  at  these  emission  processes  in  a way  which  is  differ-  ! 

ent  from  the  disciplines  described  above,  although  it  has  some  aspects  in  common  with  each. 

Plume  studies  require  the  accurate  computation  of  the  magnitude  of  the  radiation  from  a given 
"cloud"  of  water  vapor  and  C02,  at  some  position  remote  to  that  "cloud,"  as  a function  of 
wavelength  or  frequency.  The  temperatures  of  these  sources  of  water  vapor  and  C02  are 
usually  high  with  respect  to  ambient  temperature,  ranging  from  a few  hundred  degrees  K to 
more  than  2000°K.  The  spectral  resolution  required  of  these  studies  depends  upon  the  device 
to  which  the  plume  emission  studies  are  to  be  related.  A spectral  resolution  which  is  fine 
with  respect  to  the  total  spectral  extent  of  the  vibration -rotation  bands,  but  coarse  with  re- 
spect to  the  width  or  spacing  of  individual  lines  within  those  bands,  is  usually  adequate  for 
quantities  measurable  with  field  equipment.  Thus  the  presence  of  great  ranges  of  temperature 
distinguish  plume  emission  studies  from  meteorological  studies.  The  expectation  that  the  ob- 
servation site  will  be  well  separated  from  the  cloud  of  gas  means  that  there  will  be  an  inter- 
vening path  also  containing  C02  and  water  vapor.  This  adds  a complexity  which  is  not  usually 
present  in  hot  gas  radiative  transfer  studies.  The  study  of  plume  emission  processes  is  there- 
fore a mixture  of  the  considerations — from  meteorological  studies — of  absorption  and  emis- 
sion of  atmospheric  species  (which  arise  from  the  intervening  atmospheric  path)  and  of  the 
study  of  the  radiative  transfer  of  hot  gases. 
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The  required  spectral  resolution  dictates  the  type  of  treatment  which  is  appropriate  for 
the  spectral  character  of  the  radiation.  For  most  plume  studies,  extremely  high  spectral 
resolution  is  not  required.  The  smallest  interva’  oi  interest  is  small  with  respect  to  the 
whole  band  but  large  with  respect  to  individual  line  widths  and  spacing.  This  implies  a reso- 
lution of  a few  cm”1  (say,  5 to  50),  whereas  the  vibration-rotation  bands  are  a few  hundred 

cm”1  in  width.  The  two  bands  of  primary  interest  to  plume  studies  are  (a)  the  2.7  f im  H00 

-1  1 

and  C00  band,  which  extends  roughly  from  3200  to  4100  cm  (2.4  to  3.1  |im),  and  (b)  the  4.3 

^ -1 
pm  band  of  C02  which  extends  from  about  2150  to  2400  cm  (4.17  to  4.65  pm).  The  2.7  pm 

HgO  band  is  a combination  of  two  fundamental  bands  and  an  overtone  band  (t^,  fundamentals 
and  2i/2);  the  4.3  pm  C02  band  is  the  fundamental  band.  The  C02  band  at  2.7  pm,  in  the 
same  region  as  the  water  vapor  band,  is  the  result  of  a number  of  combination  bands.  The 
typical  width  of  the  lines  in  these  bands,  at  normal  plume  temperatures  and  atmospheric  pres- 
sures, is  of  the  order  of  several  hundredths  of  a wavenumber.  In  the  very  finest  frequency 
scale,  one  would  expect  substantial  absorption  (and  emission)  variations  over  a fraction  of 
the  line  width.  Thus,  it  would  require  spectral  resolution  of  the  order  of  0.01  cm”1  to  see 
these  variations.  (This  resolution  is  beyond  the  capability  of  all  but  the  largest  of  laboratory 
instruments.) 

The  average  spacing  of  the  lines  in  these  bands  if  of  the  order  of  a wavenumber  if  only 
the  strongest  lines  are  considered,  and  substantially  less  if  the  weak  lines  are  included.  For 
example,  the  very  complete  tabulation  of  atmospheric  absorption  lines  prepared  by  McClatchey 
et  al.,^  contains  about  250  water  vapor  lines  in  a 40  cm”1  wide  interval  centered  about 
3550  cm”1  (2.82  pm). 

Thus  a treatment  of  the  spectral  dependence  in  which  any  given  resolution  interval  is  as- 
sumed to  contain  the  integrated  effect  of  a large  number  of  lines  is  appropriate.  Any  single 
resolution  interval  will  contain  many  lines,  and  the  total  effect  of  any  given  line  is  largely 
concentrated  within  the  resolution  element  in  which  it  falls.  Treatment  of  the  frequency  de- 
pendence in  this  manner  is  called  a "band  model"  treatment.  Band  model  techniques  have 
been  developed  as  a means  of  averaging  over  the  very  finest  spectral  structure  of  absorption 
and  emission.  Band  model  treatment  of  the  plume  emission  is  employed  throughout  this  study. 

Several  good  surveys  of  band  model  techniques  for  use  with  gases  near  ambient  tempera- 
ture are  available.  Goody ^ provides  a fundamental  description  of  the  processes  of  gaseous 
emission  and  absorption  leading  to  the  development  of  band  models.  Anding^  presents  a 

2.  R.A.  McClatchey,  et  al.,  AFCRL  Atmospheric  Absorption  Line  Parameters  Compilation, 
AFCRL-TR-73-006,  Air  Force  Cambridge  Research  Laboratories,  Bedford,  Massachusetts, 
January  1973. 

3.  R.M.  Goody,  Atmospheric  Radiation  I:  Theoretical  Basis,  Clarendon  Press,  Oxford, 
England,  196C 

4.  D.  Anding,  Band  Model  Methods  for  Computing  Atmospheric  Slant -Path  Molecular  Absorp- 
tion,  Report  No.  7142-21 -T,  Willow  Run  Laboratories,  Ann  Arbor,  Michigan,  February  1967. 
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detailed  comparison  of  various  band  models  with  available  experimental  data.  These  two  ref-  i 
erences  show  that  there  are  a number  of  available  band  model  techniques  which  predict  the  | 

transmittance  of  water  vapor  and  CO2  paths  near  atmospheric  temperature  and  pressure  | 

quite  well.  f 

Band  models  have  also  been  developed  specifically  for  the  treatment  of  hot  gas  radiation.  | 
Penner^  provides  a well-known  introduction  to  high  temperature  gas  emission  in  which  | 

techniques  closely  related  to  band  models  are  used.  Ludwig^  provides  a more  recent  treat-  | 

ment  in  which  sufficient  numerical  values  are  included  so  that  reasonably  accurate  hot  gas  j 

emission  calculations  can  be  made  without  reference  to  further  numerical  data.  The  numer-  j 

ical  values  included  in  Reference  1 provide  the  data  base  for  much  of  the  numerical  calcula-  j 

tions  made  on  this  program.  j 

J 

Neither  the  atmospheric  band  models  nor  the  hot  gas  band  models,  when  taken  alone,  ; 

provide  an  adequate  treatment  of  plume  emission  problems  in  the  presence  of  an  intervening  ! 

atmosphere.  The  reason  for  this  is  that  transmittance  calculations,  for  which  the  atmospheric 
band  models  are  well  suited,  assume  that  the  source  of  radiation  is  spectrally  continuous. 

That  is,  the  spectral  distribution  of  the  source  radiation  is  assumed  to  be  uniform  over  a 
spectral  interval  at  least  as  wide  as  the  spectral  resolution  of  the  band  model.  This  approx-  j 

imation  is  well  satisfied  for  most  opaque  thermal  emitters,  such  as  the  terrain  features 
which  are  the  most  significant  source  of  thermal  radiation  in  the  spectral  region  beyond  3 /im. 

In  contrast,  the  emission  from  plume  gases,  being  largely  due  to  the  same  vibration - 
rotation  bands  which  give  rise  to  atmospheric  absorption,  are  (spectrally)  highly  discontin- 
uous and  consist  of  spectral  line  emission — in  many  cases,  the  same  lines  in  which  atmo- 
spheric absorption  occurs.  The  emission  lines  are  thus  correlated  with  the  absorption  lines 
in  many  instances.  A band  model  treatment,  by  definition,  averages  over  the  fine  line  struc- 
ture. Hence,  It  is  important  that  any  band  model  procedure  which  treats  a situation  in  which 
the  emission  and  absorption  lines  are  correlated  be  structured  to  carry  out  this  averaging 
after  the  whole  path  is  treated.  In  particular,  procedures  which  average  over  the  line  struc- 
ture of  the  source  alone  and  then  over  the  structure  of  the  atmospheric  absorption  alone, 
without  any  consideration  of  the  effects  of  these  correlations,  will  yield  questionable  results 
where  line  correlation  effects  are  known  to  be  present. 

These  line  correlation  effects  are  the  primary  phenomenon  investigated  in  the  present 
program.  Although  correlation  effects  are  well  known,  they  are  often  overlooked  in  many 
computations  of  plume  emission.  The  reason  is  that  there  has  been  no  definitive  treatment 


5.  S.S.  Penner,  Quantitative  Molecular  Spectroscopy  and  Gas  Emissivities,  Addison - Wesley, 
Reading,  Massachusetts,  1959. 


17 


of  their  magnitude.  Burch  et  al.,^  made  limited  measurements,  but  provided  insufficient 
information  by  which  the  importance  of  the  demonstrated  effects  could  be  Judged. 

Calculations  have  beemmade  which  attempt  to  estimate  correlation  effects  in  the  atmo- 
spheric attenuation  of  hot  gas  emission.  At  least  for  water  vapor  the  predicted  effects  were 

F71 

large.  In  197Q,  band-model  calculations  were  made1  1 with  a conscious  effort  to  include  the 
atmospheric  path  together  with  a path  through  a plume  as  one  long,  non-isothermal  path.  The 
objective  was  to  avoid  averaging  over  the  line  structure  until  after  the  attenuation  of  the  at- 
mospheric path  had  been  included.  This  model  contained  the  implicit  assumption  that  the 
lines  were  perfectly  correlated  and,  hence,  that  the  attenuation  produced  by 'the  atmospheric 
path  was  much  more  severe  than  the  attenuation  of  continuum  or  greybody  emission.  Subse- 
quently, a band  model  treatment  was  developed^  in  which  the  non -homogeneous  character 
of  the  non-isothermal  path  was  considered  but  in  which  the  lines  were  still  assumed  to  be 
completely  correlated.  Young  extended  these  ideas^  but  still  did  not  account  for  partial 
correlation.  The  results  from  these  efforts  indicate  severe  correlation  effects.  At  the  time 
that  the  work  reported  in  References  7 and  8 was  performed,  it  was  not  apparent  to  the 
authors  that  the  assumption  of  complete  correlation  had  been  made.  During  the  course  of  the 
present  analysis  it  became  clear,  however,  that  complete  correlation  was  indeed  inadvertently 
assumed  in  the  above-mentioned  work. 

There  is  no  reason  to  believe  that  the  source  emission  lines  and  atmospheric  absorption 
lines  are  completely  correlated.  At  high  temperatures  energy  states  above  the  ground  state 
become  populated,  leading  to  many  lines  in  emission  which  are  not  significantly  present  in 
the  absorbing  gas.  No  work  had  been  performed  to  assess  the  degree  to  which  the  emission 
and  absorption  lines  correlate.  No  significant  measurements  of  gas  emission  as  viewed 
through  a cool  path  have  been  previously  made. 

1.2  OBJECTIVES  AND  APPROACH 

The  objectives  were  to  determine  the  magnitude  of  spectral  correlation  between  a target 
plume  and  the  atmosphere,  as  well  as  to  develop  methods  whereby  atmospheric  absorption, 
including  correlation  effects,  could  be  adequately  treated  in  techniques  to  compute  plume 
radiant  emission.  Such  methods  should  be  useful  to  calculate  the  irradiance  at  a sensor, 


6.  D.E.  Burch  and  D.A.  Gryvnak,  Infrared  Radiation  Emitted  by  Hot  Gases  and  Its  Trans- 
mission through  Synthetic  Atmospheres,  Report  No.  U -1929,  Research  Laboratory,  Aero- 
nutronic  Division,  Ford  Motor  Company,  Newport  Beach,  California,  October  1962. 

7.  F.S.  Simmons,  et  al.,  "Calculation  of  Radiation  from  Hot  H2O  and  CO2  Viewed  Through  a 
Cool  Intervening  Atmosphere,"  Applied  Optics,  Vol.  9,  p.  2792f  (1970). 

8.  G.H.  Lindquist  and  F.S.  Simmons,  "A  Band  Model  for  Very  Non-Uniform  Paths,"  JQSRT, 
Vol.  12,  p.  807 f (1972). 

9.  S.J.  Young,  Band  Model  Formulation  for  Inhomogeneous  Paths,  SAMSO-TR-74-255,  The 
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after  the  target  radiation  has  been  modified  by  the  atmosphere.  The  following  three  tasks 
were  undertaken  to  satisfy  this  objective. 

Task  I.  Band  models  were  developed  to  tfeat  the  nonhomogeneities  present  in  the  plume* 

atmospheric  path  combination  and  to  represent  line  correlation  effects  In  the  hot  and 
cold  gases. 

Task  II.  Measurements  were  made  which  provided  a characterization  of  the  correlation 
effects  and,  hence,  could  be  used  to  assess  the  accuracy  of  the  band  models  developed 
under  Task  I,  as  well  as  that  of  the  older  band  model  techniques.  The  emission  of  a 
small  hot  cell  filled  with  mixtures  of  hot  water  and  CO2  was  measured  both  unattenuated 
and  attenuated  by  a long  cell  filled  with  cool  atmospheric  gas.  The  path  conditions  were 
chosen  to  represent  paths  corresponding  to  realistic  plume  observations.  These  measure- 
ments were  made  in  both  the  HgO-COg  band  at  2.7  |xm  and  the  C02  band  at  4.3  (im.  The 
resolution  of  the  spectral  apparatus  was  chosen  to  be  appropriate  for  band  model  methods. 

Task  III.  The  measurements  made  in  Task  II  were  compared  to  the  available  predictions. 
Comparisons  were  made  to  the  band  models  derived  in  Task  I,  to  previously  developed 
band  models,  and  to  line-by-line  calculations  based  on  the  tabulation  of  atmospheric  ab- 
sorption lines  of  Reference  1. 


THE  DEVELOPMENT  OF  BAND  MODELS  TO  TREAT 
CORRELATION  EFFECTS  AND  PATH  INHOMOGENEITIES 

The  band  model  development  described  here  is  patterned  after  the  treatment  of  atmo- 
spheric band  models  by  GoodyJ3^  It  is  the  result  of  a great  deal  of  thought  as  to  how  to  best 
treat  line  correlation  effects.  It  was  originally  hoped  that  an  entirely  different  type  of  statis- 
tical treatment  of  line  correlation  effects  would  be  useful,  but  after  several  unfruitful  efforts, 
the  ideas  described  here  were  pursued.  They  lead  to  a realistic  means  for  the  treatment  of 
correlation  effects  and  are  sufficiently  general  so  that  a large  set  of  physical  conditions  can 
be  treated  provided  suitable  extensions  to  the  model  are  derived.  Only  one  band  model  is 
developed.  It  is  a random  line  model  in  that  the  lines  are  assumed  to  be  randomly  positioned. 

We  originally  planned  to  develop  a model  specifically  for  regularly  spaced  lines.  But 
after  preliminary  efforts,  we  were  convinced  that  no  physical  situation  was  present  where  a 
regular  band  model  would  offer  significant  advantages  over  a random  model.  Our  original 
plan  for  a regular  model  was  based  on  the  fact  that  lines  in  individual  bands  of  C02  are  nearly 
regularly  spaced.  Nevertheless,  the  presence  of  many  individual  C02  bands,  particularly  at 
higher  temperatures,  means  that  the  emission  spectrum  is  a superposition  of  a large  number 
of  bands  with  nearly  regularly  spaced  lines.  Line  spaclngs  vary  among  the  bands.  Each  of 


these  bands  is  nearly  randomly  spaced  with  respect  to  every  other.  Since  the  result  is  a ran- 
dom array  of  lines  in  some  sense,  one  would  expect  that  a random  model  would  accurately 
treat  this  situation  also,  particularly  since  many  lines  are  included  in  an  averaging  interval 
(i.e.,  in  a single  spectral  resolution  element).  Furthermore,  as  the  results  of  Task  n show, 
there  were  only  minor  line  correlation  effects  present  in  COg  under  the  conditions  of  this 
experiment.  This  means  that  it  would  be  impossible  to  distinguish  between  the  performance 
of  the  random  model  and  that  of  a regular  model  in  predicting  line  correlation  effects.  We 
did  not  consider  a regular  model  for  use  with  water  vapor  because  of  the  demonstrated  ran- 
domness of  the  water  vapor  band  (see  Reference  3 and  the  following  sections). 

The  development  given  here  is  restricted  to  collision-broadened  lines.  Collisi'  ki  broad- 
ening is  dominant  in  low  altitude  plumes  viewed  through  paths  which  do  not  extend  above  sev- 
eral tens  of  kilometers  in  altitude.  Furthermore,  the  experimental  apparatus  ueed  in  Task  n 
was  not  sensitive  enough  to  allow  the  hot  cell  pressures  to  be  low  enough  such  that  collision 
broadening  was  not  significant.*  Thus,  the  treatment  here  is  limited  to  collision -broadened 
lines.  The  following  sections  give  the  detailed  development  of  a random  band  model  which  in- 
cludes the  effect  of  line  correlation  and  which  can  treat  path  Inhomogeneities. 

2.1  THE  EQUATION  OF  TRANSFER 

Elementary  considerations  of  the  interactions  of  radiation  with  matter  lead  to  the  so- 
called  equation  of  transfer  (see  Reference  3,  p.  25),  which  treats  all  interaction  processes  as 
either  extinction  or  absorption  and  which  can  be  written  in  the  following  form: 

dL~  = (-L~e~  + J~e~)  da  (1) 

V V V V V 

Here  L~  is  the  spectral  radiance  (called  intensity  by  Goody  and  other  astrophysicists).  Spec- 
tral radiance  is  the  power  emitted  (or  exiting)  from  a surface  per  unit  cross-section  per 
unit  solid  angle,  per  unit  spectral  Interval;  for  wavenumber  units,  it  is  commonly  specified 

2 -i 

in  watts/(cm  -sr-cm  ).  The  quantity  e~  is  the  extinction  coefficient  and  J~  is  called  the 
source  function.  The  Increment  da,  is  the  amount  of  matter  in  an  element  of  path  of  unit 
cross -section  along  which  radiation  is  transferred,  and  v is  the  optical  wavenumber  (or  al- 
ternatively, frequency  or  wavelength). 

The  first  (extinction)  term  assumes  Lambert's  law  of  extinction,  namely  that  the  amount 
of  extinction  is  linear  with  both  radiance  and  wit!,  'he  amount  of  absorbing  matter  in  the  path. 
(Radiation  processes  which  are  nonlinear  in  nature  are  important  in  many  current  problems, 


*Actually,  Doppler  and  collision-broadened  halfwidths  were  of  the  same  order  in  the 
measurements.  Because  of  the  great  difficulty  of  developing  a model  to  treat  mixed  broaden- 
ing, and  because  the  collision  effects  still  dominate  when  the  widths  are  of  the  same  order, 
the  development  of  a model  for  collision -broadened  lines  as  opposed  to  Doppler -broadened 
lines  is  appropriate. 
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but  the  flux  required  is  very  much  larger  than  the  fluxes  to  be  expected  from  thermal  targets 
of  a few  thousand  degrees  Kelvin  such  as  aircraft  and  rocket  plumes.)  The  form  of  the  second 
(emission)  term  also  reflects  such  an  assumption  of  linearity  and  can  be  considered  as  the 
definition  of  the  source  function. 

As  Goody^  states,  this  equation  sets  the  pattern  for  all  formalisms  used  in  the  solution 
of  radiative  transfer  problems,  although  its  physical  content  is  slight.  The  physical  interpre- 
tation of  the  meaning  of  the  equation  is  as  follows:  L~  is  the  spectral  radiance  along  the  line 
of  sight  and  in  Equation  (1)  can  be  considered  to  be  a function  of  the  path  coordinate  a.  The 
quantities  e~  and  J~  can  also  be  considered  functions  of  the  variable*a."  Note  that  the  subscript  j 

implies  a wavenumber  dependence  as  well.  The  variable* a*can  be  taken  as  the  path  coordinate  j 

1 

since  it  is  defined  as  the  amount  of  absorbing  matter  per  unit  path  cross  section  between  some  , 
point  on  the  path  (the  point  of  observation)  and  a general  physical  point  along  the  line  of  sight. 

♦ ~ I 

Thus,  the  variable  a serves  as  a line  of  sight  variable.  Equation  (1)  governs  the  variation  in 
radiance  along  the  line  of  sight  defined  by  values  of* a?  In  order  to  determine  the  radiant  in-  5 

tensity  of  a non -homogeneous  source  such  as  a plume,  at  an  observer  some  distance  from  the  , 

plume,  it  is  necessary  to  solve  Equation  (1)  for  all  lines  of  sight  from  the  observer  passing  ( 

n> 

through  the  plume  and  to  then  integrate  the  results  over  the  apparent  cross  section  of  the  . 

! 

plume.  ; 


In  order  to  facilitate  the  solution  of  Equation  (1),  the  following  algebraic  manipulations 
are  often  made.  We  set  a - 0 at  the  observer  and  express  the  transmittance  from  the  observer 
to  the  point  "a"  as 


•r(£,a)  = exp 


(' 


a 

I 


e~(p,a')da’] 


(2) 


! 

I 


Note  that,  because  of  the  assumed  linear  relationship  between  extinction  and  radiance  and 
amount  of  matter,  this  is  indeed  the  transmittance  of  energy  emitted  at  the  point  "a"  to  the 
observer,  at  wavenumber  v.  After  a derivative  is  taken,  this  definition  can  be  introduced  into 
the  equation  of  transfer.  The  sign  of  the  variable  a in  Equation  (1)  must  be  reversed  to  cor- 
respond  to  the  definition  of  t in  Equation  (2).  In  Equation  (2),  the  variable  a increases  along 
the  line  of  sight  going  away  from  the  observer — i.e.,  in  the  direction  of  increasing  radiance 
when  no  source  is  present  (J~  = 0);  whereas  in  Equation  (1),  a increases  in  the  direction  of 
decreasing  radiance  when  J~  - 0.  The  substitution  of  (2)  into  (1)  yields 

rfMdLyla)  = -L-<a>m^da  * ■>;;<“> 


If  d[7(!5,a)L~(a)]  is  the  differential  of  the  product  of  r and  L~  along  the  line  of  sight,  then 
the  above  equation  can  be  written 
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dtTfoalLjM]  = J-< a) da  (3) 


This  can  now  be  integrated  over  the  line  of  sight  from  the  observer  to  the  far  limit  of  the  line 
of  sight.  This  far  limit  should  extend  to  some  opaque  surface;  however,  in  most  measurement 
systems  radiation  from  portions  of  the  path  behind  the  target  of  interest  is  either  automatically 
subtracted  or  otherwise  compensated  for.  Hence,  the  far  limit  is  usually  taken  at  the  far 
boundary  of  the  target  of  interest.  (For  targets  immersed  in  an  atmosphere  and  viewed  by 
relatively  wide  field  instruments,  radiation  from  the  path  in  front  of  the  target  is  also  often 
automatically  compensated  for.)  Thus  Equation  (3)  can  be  integrated  to  yield  the  quantity 

far  limit 

t(v,  far  limit)L~(far  limit)  - L -.(observer)  = J J~(a)^^- da 

observer 

If  the  far  limit  of  the  path  does  not  yield  significant  radiances,  (i.e.,  L~(far  limit)  = 0),  then 
the  first  term  can  be  ignored,  leaving  the  contribution  from  the  gas  target  alone: 

far  limit 

L~(observer)  = - J J~(a)  da  (4) 

observer 

So  far,  nothing  has  been  said  about  the  source  function.  Under  the  conditions  of  thermo- 
dynamic equilibrium,  which  we  expect  to  hold  for  rocket  and  aircraft  exhaust  plumes  at  rea- 
sonably low  altitudes,  the  source  function  is  the  Planck  blackbody  radiance  function,  a function 
of  temperature  and  frequency  only. 

The  treatment  so  far  has  assumed  that  the  measurement  system  to  be  used  to  measure 
these  radiances  has  perfect  spectral  resolution.  But  every  practical  instrument  measures  an 
average  over  a spectral  region  either  by  means  of  its  spectral  "slit  function"  if  it  is  a dis- 
persing instrument  or  by  other  frequency -selective  mechanisms  in  case  of  radiometric  mea- 
surements. The  case  of  most  interest  to  this  particular  program  is  that  of  a detection  system 
with  moderate  resolution.  Such  an  instrument  takes  an  average  over  a frequency  interval 
large  with  respect  to  individual  molecular  absorption  line  spacings  but  small  with  respect  to 
the  width  of  a molecular  vibration -rotation  band.  Upon  integration  of  Equation  (4)  over  a fre- 
quency band  of  width  Ay,  we  obtain  for  the  average  radiance  seen  by  the  observer: 


far  limit 


L~=  — f f ju,(a)m^idad? 

v Ay  J J v ' da 
Ay  observer 

The  order  of  integration  can  be  reversed  so  that 


far  limit 

§ J jv(a)^d5'da 

observer  Ay 
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If  the  source  function  can  be  assumed  sensibly  constant  over  Ay,  as  it  can  over  the  spectral 
widths  of  25  cm”*  around  4000  cm  * which  are  of  interest  here,  then  J~,(a)  can  be  removed 
from  the  integration  over  At/,  and  we  obtain 

far  limit 


L~  = 
v 


‘ 1 J5c<a) 

“ f~ 

Ay  J 

observer 

far  limit 
r 

L A? 

dr(y  ,a) 

- J h (a) 

c 

da 

observer  c 

■&v' 


da 


da 


(5) 


Equation  (5)  is  similar  to  (4)  except  that  the  transmittance  derivative  involved  is  now  the  spec- 
tral average  of  the  transmittance  derivative  at  perfect  resolution.  ( V is  the  frequency  at  the 
center  of  At/.)  The  integration  of  Equation  (4)  or  (5)  is  often  mapped  into  an  integration  over 
t or  r.  Thus, 

far  limit 

[ J~  (7)  dr  (6) 


L~  = 
v 


observer 


or 


far  limit 

- I ' 

observer 


L~  = - J~(r)dr 

v J v 


(7) 


AH  of  the  above  manipulation  has  been  carried  out  independent  of  the  character  of  the 
emitting  matter.  The  only  restrictions  involved  are 

(1)  Only  linear  processes  are  involved  in  the  extinction  and  emission  processes.  This 
also  implies  that  the  perfectly  resolved  transmittances  will  obey  Beer's  law  (Equa- 
tion 2).  But  the  average  transmittances  are  not  thus  restricted.  An  average  of  a 
number  of  exponentials  like  Equation  (2),  each  of  different  e~,  does  not  necessarily 
lead  to  another  exponential.  In  fact,  most  band  model  expressions  for  spectrally 
averaged  r are  not  exponentials. 

(2)  The  conditions  of  thermodynamic  equilibrium  are  satisfied  everywhere  along  the  path 
of  interest. 

(3)  If  Equation  (5)  or  (6)  Is  to  be  used,  the  Planck  blackbody  function  must  be  essentially 
constant  over  the  spectral  region  Ay. 

The  above  Equations  (4),  (5).  (6),  and  (1)  are  therefore  valid  for  a wide  variety  of  situations 
in  which  the  radiance  seen  at  some  observing  position,  due  to  sources  distributed  along  the 
line  of  sight,  is  to  be  calculated.  The  remainder  of  this  section  applies  these  particular 
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equations  to  calculations  of  gas  radiation  in  the  molecular  v ibratlon  -rotation  bands  of  CO2  and 
HgO  which  are  of  primary  interest  to  this  program. 


2.2  RADIANCE  AT  A DISTANCE  DETERMINED  BY  A MULTIPLICATIVE  TRANSMITTANCE 
In  many  situations  of  this  type,  the  radiance  at  a distance  (through  an  atmospheric  path) 
is  computed  by  first  computing  the  radiance  of  the  gaseous  target  alone  and  then  multiplying 
this  by  the  transmittance  of  the  atmospheric  path.  Such  a process  decouples  the  gas  radiation 
from  the  atmospheric  transmittance.  It  is  useful  at  this  point  to  compare  this  procedure  with 
the  result  of  Equations  (4)  through  (7).  Consider  the  line  of  sight  shown  in  Figure  1(a).  The 
segment  from  0 to  ap  is  the  atmospheric  path  and  that  from  ap  to  ap+t  the  hot  gas  of  interest. 
Assume  that  the  radiance  contribution  from  the  atmospheric  path  (from  0 to  ap)  is  either 
small  or  compensated  for  by  the  detection  system.  The  target  contribution  given  by  Equa- 
tion (5)  is  thus 

Vt 

V ■ | Va')^>da'  (8a) 


The  decoupled  procedure  mentioned  above  corresponds  to  the  evaluation  of  Equation  (5) 
with  the  observer  placed  at  the  boundary  of  target,  as  in  Figure  1(b).  The  result  is  multiplied 
by  the  transmittance  from  0 to  ap.  For  this  case, 

at 

L~  = -30  to  ap)- 1 J~(a’)  ^ da'  (8b) 

0 


Assuming  that  the  target  is  essentially  isothermal,  the  radiances  of  Equations  (8a)  and  (8b) 
are  equal  only  if 

-Ro  to  ap)*[l  - r(0  to  at)]  = t{0  to  ap)  - t^O  to  ap+t) 
or 

t(0  to  at)*7<0  to  ap)  = r(0  to  a fc)  (9) 

The  transmittance  of  the  whole  path  (target  plus  atmosphere)  must  equal  the  product  of 
the  transmittance  of  the  atmosphere  alone  and  the  transmittance  of  the  target  alone. 

There  are  only  a limited  number  of  situations  where  this  can  be  true.  For  instance,  if 
the  transmittances  are  given  by  simple  exponentials  of  an  optical  depth,  relation  (9)  is  easily 
shown  to  hold.  But  the  averaged  transmittances  given  by  band  model  expressions  rarely  are 
represented  by  such  exponentials.  Yet,  there  are  other  cases  where  relation  (9)  holds.  The 
transmittances  are  each  averaged  over  many  lines.  If  the  lines  in  one  portion  of  the  path 
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(a)  Cool  Absorbing  Path  and  Hot  Gas  Plume  Treated  as  a Single  Non-Isothermal  Path 
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(b)  Cool  Absorbing  Path  and  Hot  Gas  Plume  Treated  as  Separate  Decoupled  Paths 


FIGURE  1.  TREATMENT  OF  PLUME -ATMOSPHERE  COMBINATION  BY  TWO  DIFFERENT 

CALCULATIONAL  TECHNIQUES 


(atmosphere)  are  on  the  average  uncorrelated  with  those  in  the  other  portion  (target),  then 
the  average  ot  the  product  of  two  transmittances  will  equal  the  product  of  the  averages.  This 
situation  can  be  expected  to  hold  when  the  atmospheric  species  producing  the  absorption  are 
not  the  same  as  the  emittirg  species  of  the  plume  and  when  the  numbers  of  target  lines  and  of 
atmosphere  lines  in  a resolution  interval  are  large.  The  transmittances  of  different  gases, 
averaged  over  many  lines,  can  be  multiplied  because  the  correlation  between  any  group  of 
lines  of  the  two  molecules  will  be  small. 

Relation  (9)  also  holds  when  either  or  both  of  the  transmittances  is  sensibly  constant  over 
the  integration  interval.  These  transmittances  of  continuum  emitter -absorbers  can  always 
be  safely  multiplied  by  other  transmittances. 

Despite  the  many  cases  where  relation  (9)  holds,  it  cannot  hold  where  the  absorbing  and 
emitting  species  are  identical  and  where  the  band  model  relations  for  the  average  transmit- 
tances deviate  from  simple  exponentials.  Such  is  certainly  the  case  for  water  vapor;  less  so 
for  COg. 

Nor  can  relation  (9)  hold  when  the  resolution  interval  is  so  small  that  it  includes  only  a 
few  lines  from  either  the  emitting  gas  or  the  atmospheric  absorbing  gases.  Thus  the  product 
rule  cannot  be  expected  to  hold  for  a spectral  interval  which  includes  only  a few  lines  from  a 
species  like  HF  (whose  lines  are  widely  separated),  even  when  there  are  many  atmospheric 
lines  in  the  interval. 

2.3  APPLICATION  OF  EQUATIONS  5 AND  6 TO  TARGET  PROBLEMS  OF  INTEREST 

The  procedure  discussed  in  the  previous  section,  in  which  the  target  emission  and  the 

atmospheric  attenuation  are  decoupled,  is  incorrect  in  general  and,  although  it  is  still  com- 
monly used,  it  will  not  be  considered  further  in  this  section.  Comparisons  of  experimental 
measurements  and  calculations  made  using  this  assumption  will  be  discussed  in  a later  sec- 
tion. The  primary  consideration  of  the  present  section  will  be  to  lay  out  the  means  to  be  used 
to  evaluate  Equations  (5)  and  (6)  which  correctly  represent  moderate  resolution  radiances  of 
gaseous  targets  at  a distance.  The  quantity  in  these  equations  which  must  be  determined  is 
the  transmittance  derivative  along  the  path  in  the  case  of  Equation  (5)  and  the  essentially 
equivalent  transmittance  differential  in  Equation  (6). 

2.4  PRIOR  EFFORTS  TO  COMPUTE  TRANSMITTANCE  DERIVATIVES 

Equation  (6)  is  used  in  many  existing  computer  programs,  which  take  first  differences  of 
successive  transmittance  calculations,  each  made  for  a path  from  the  observer  to  a general 
point  in  the  line  sight.  Generally  the  Curtis -Godson  approximation  is  used  to  account  for  path 
non-homogeneities.  Lindquist  and  Simmons^  and  Young^  show  that  the  effective  transmit- 
tance derivative  obtained  by  this  procedure  leads  to  incorrect  and  even  physically  unreasonable 
results  in  some  cases.  Lindquist  and  Simmons  also  determined  a more  accurate  transmittance 
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derivative  in  which  the  Curtis -Godson  approximation  was  applied  to  the  transmittance  deriva- 
tive rather  than  to  the  transmittance  itself.  A heuristic  extension  to  a random  model  was 
also  made.  Young  further  extended  this  formulation  to  several  distributions  of  line  strengths 
but  did  not  account  for  the  fact  that  different  lines  change  strength  with  temperature  (and 
hence  along  the  path)  in  different  ways. 

It  became  apparent  early  in  the  experimental  portion  of  this  contract  that  the  Lindquist  - 
Simmons  approximation  does  not  do  substantially  better  at  predicting  hot  gas  radiance  through 
a cool  path  than  does  the  very  naive  approximation  described  earlier  (hot  gas  radiance  mul- 
tiplied by  cool  gas  transmittance).  Therefore,  development  of  an  expression  for  the  spectrally 
averaged  transmittance  derivative  was  undertaken  using  band  model  procedures.  The  results 
of  this  effort  came  too  late  to  allow  complete  comparisons  to  be  made  of  the  resulting  mathe- 
matical formulations  with  the  experimental  results.  Section  4 describes  those  comparisons 
which  were  made.  The  following  section  describes  the  formulation. 

2.5  TRANSMITTANCE  DERIVATIVE  FOR  A BAND  MODEL 

In  all  of  the  molecular  vibration-rotation  bands  of  interest  here,  and  in  particular  in  the 
2.7  pm  bands  of  HgO  and  C02  and  the  4.3  pm  band  of  COg,  each  band  consists  of  a very  large 
number  of  emission -absorption  lines.  As  the  temperature  increases,  the  population  of  vibra- 
tional states  above  the  ground  state  increases.  This  fact  causes  so  called  "hot  bands"  to  be- 
come important.  These  hot  bands  are  t,ands  which  involve  the  same  change  in  vibrational 
state  as  the  nominal  ground  state  band  but  which  have  an  initial  vibrational  level  above  the 
ground  vibrational  state.  (The  2.7  pm  water  band  results  from  the  and  fundamental 
(Av  = 1)  vibrational  transition,  as  well  as  the  first  overtone  (Av  --  2)  of  i»2  transitions  of  water 
vapor;  tho  4.3  pm  band  of  C02  is  the  fundamental.)  In  addition,  higher  temperatures  re- 
sult in  the  population  of  more  rotational  levels,  which  in  turn  causes  more  lines  of  each 
vibration -rotation  band  to  have  significant  strength. 

Though  "hot  bands"  occur  in  the  same  nominal  spectral  position  as  the  ground  state  bands, 
there  is  no  reason  to  expect  the  lines  of  the  hot  bands  and  ground  state  bands  to  exactly  coin- 
cide in  position.  Thus  the  physical  picture  of  a vibration -rotation  band  presented  by  any 
molecule,  at  temperatures  substantially  above  ambient,  is  that  of  a number  of  arrays  of  lines, 
one  for  each  band  and  all  essentially  uncorrelated  with  each  other.  In  the  case  of  water  vapor, 
the  lines  in  each  band  are  essentially  randomly  distributed.  This  results  from  the  very  com- 
plicated nature  of  the  rotational  structure  of  this  molecule.  Each  band  for  C02  consists  of 
approximately  equally  spaced  lines.  Using  this  picture,  a general  expression  for  the  trans- 
mittance derivative  will  be  derived.  Because  the  different  bands  are  essentially  uncorrelated 
with  each  other,  random  distributions  of  lines  result  for  both  molecules. 

The  following  picture  is  suitable  for  water  vapor.  Consider  a spectral  interval  Lv  having 
within  it  N lines  of  varying  strength  and  randomly  positioned.  Let  us  assume  that  each  of 
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these  lines  is  one  member  of  a random  array  of  equal  strength  lines  of  average  spacing  Ay. 
The  superposition  of  an  appropriate  N of  these  arrays  yields  the  arrangement  of  lines  in  Ay 
described  above.  The  average  spacing  of  these  lines  is 


6 = 


Ay 

N 


Goody  derives  the  transmittance  of  such  an  arrangement  of  lines  to  be 


r = exp 


I 


(10) 


where  W.  is  the  equivalent  width  (in  the  appropriate  wavenumber  or  frequency  unit)  of  the  ith 

* * v> 

line  in  Ap.  The  derivative  of  this  expression  with  respect  to  the  path  variable  a is 


dr  dT 
da  “ da 


r dW 
da 

Note  that  the  first  equality  arises  because  the  derivative  of  an  average  with  respect  to  an  in- 
dependent variable  is  always  the  average  of  the  derivative.  Thus  it  is  now  necessary  to  find 
an  expression  for  the  average  derivative  of  the  equivalent  width  of  the  lines  in  AP. 


At  this  point  it  is  necessary  to  consider  the  type  of  broadening  to  be  expected.  Since  the 
experiments  were  done  in  a region  where  pressure  broadening  is  most  important,  correspond- 
ing to  targets  below  about  20  km  altitude,  Lorentz  line  profiles  will  be  assumed  throughout. 
From  Reference  7, 
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where  and  are  the  absorption  strength  and  the  width  of  the  i line.  The  absorption 

strength  is  given  in  terms  of  the  appropriate  wavenumber  units  over  optical  pathlength,  cra'V 
o 

(molecule/cm  ).  The  halfwidth  is  given  in  wavenumber  units.  From  this  point,  since  we  have 
no  further  need  for  a mathematical  symbol  for  the  resolution  interval,  Ay  will  be  used  as  a 
dummy  variable  to  represent  the  spectral  distance  to  line  center. 
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An  average  for  this  expression  must  now  be  obtained  over  the  N lines  in  the  interval. 
(Note  that  we  have  assumed  that  all  of  the  equivalent  width  is  sensibly  included  in  the  resolu- 
tion interval;  since  line  widths  less  than  0.1  cm"*  are  expected  in  an  interval  of  about  10  to 
30  cm"*,  this  assumption  will  be  satisfied  for  almost  all  lines.)  Young^  has  obtained  the 
required  average  by  applying  the  Curtis -Godson  approximation  to  the  exponential  argument 
(attenuation  term)  so  that  for  each  line  two  strengths  appear,  the  local  strength  S^a)  and  the 
average  strength  Se^(a).  He  has  assumed  all  and  to  vary  with  temperature  in  identical 
fashion  and  has  averaged  over  Equation  (12)  using  different  line  strength  distributions.  In 
this  effort  we  would  like  to  account  for  the  fact  that  the  strengths  of  the  different  lines  in  the 
interval  grow  in  different  ways  with  temperature,  depending  on  their  lower  state  energy. 


The  growth  of  the  strength  of  a line  with  temperature  is  due  solely  to  the  change  in 
equilibrium  populations  of  the  upper  and  lower  states  involved.  This  growth  is  given  by  the 
following  expression: 


(13) 


where  SQ  is  the  strength  at  Tq,  E is  the  lower  state  energy  in  the  appropriate  units,  m is  a 
quantity  depending  on  the  rotational  degeneracy  of  the  molecule  involved,  and  k is  Boltzmann's 
constant.  For  HjO,  m = 3/2;  for  COg,  m = 1.* 

It  is  also  necessary  to  consider  the  Lorentz  halfwidths  of  the  lines  involved,  their  distri- 
bution, and  the  degree  to  which  they  are  dependent  or  independent  of  the  strengths  of  the  lines. 
In  this  report,  the  halfwidths  will  all  be  assumed  identical.  This  is  commonly  done  in  band 
model  treatments,  and  the  introduction  of  varying  line  widths  would  greatly  complicate  the 
formulations. 


The  line  widths  in  a given  band  vary  slowly  across  the  band.  Thus,  for  any  single  band, 
the  assumption  of  constant  line  width  in  any  resolution  interval  is  reasonable.  It  becomes 
more  difficult  to  Justify  the  use  of  constant  line  width  when  many  hot  bands  become  important. 
In  a later  section  of  the  derivation,  the  absorption  in  the  square  root  limit  is  set  equal  to  the 
correct  square  root  limit,  in  which  the  effects  of  variable  line  width  are  accounted  for.  Thus 
the  effects  of  line  width  variation  are  limited  to  the  transition  between  linear  and  square  root 
growth.  This  transition  region  is  also  sensitive  to  a number  of  other  assumptions  of  the  band 
model;  and  since  line  width  variations  (to  the  extent  that  they  are  known)  are  less  than  a fac- 
tor of  3 in  the  wing  of  the  2.7  jxm  band,  it  seems  reasonable  to  neglect  these  variations.  Exist- 
ing band  model  formulations  for  isothermal  transmittance  almost  always  ignore  these  varia- 
tions without  compromising  their  prediction  capability.  However,  to  make  the  same  statement 


*The  vibrational  partition  function  has  been  neglected  in  this  expression.  Its  inclusion 
does  not  change  the  resulting  formulations. 
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for  the  present  formulation  would  require  further  testing.  For  the  present,  we  assume  iden- 
tical halfwidths  for  all  the  lines  in  die  interval  but  allow  them  to  vary  according  to  the  tem- 
perature and  pressure  variation  encountered  along  the  line  of  sight. 

Rather  than  substitute  the  complicated  form  of  Equation  (13)  into  the  expression  for  the 
average  equivalent  width  derivative  at  this  point,  we  will  only  note  the  functional  dependencies 
of  Equation  (13): 
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SQ  and  E are  now  the  randomly  varying  parameters  associated  with  the  N lines  in  a given 
resolution  interval.  The  function  f(SQ,E)  is  a joint  probability  density  function  for  SQ  and  E. 
Note  also  that,  in  this  expression,  the  line  positions  are  assumed  to  be  independent  of  Sq  and 
E,  an  assumption  which  seems  reasonable  considering  the  physical  picture  described  earlier. 
Inserting  Equation  (12)  into  Equation  (14)  yields 

So  E - 00 

f(S0,E)d(Ay)dSQdE  (15) 

In  general  the  integrations  over  Sq,  E,  and  Ay  can  be  interchanged  at  will,  but  the  integration 
over  a'  is  within  the  exponential  and  cannot  be  interchanged.  However,  it  is  this  integral 
which  we  would  like  to  perform  last.  In  other  words,  we  seek  a formulation  in  which  the  in- 
tegrations over  Ay,  SQ,  and  E can  be  performed  once  for  all,  so  that  integrations  over  differ- 
ent physical  paths  can  be  performed  at  will  and  without  requiring  separate  integrations  over 
Ay,  SQ,  and  E.  This  has  been  the  function  of  the  Curtis  -Godson  approximation  in  calculations 
of  transmittance.  Armstrong  has  provided  a general  treatment  of  this  approximation^10^ ; 
and  because  SQ,  E,  and  y all  appear  in  the  integral,  a somewhat  more  complicated  approxima- 
tion to  the  integral  than  discussed  by  Armstrong  must  be  made.  The  Curtis -Godson  approxi- 
mation had  been  intuitively  derived  in  the  past;  Armstrong  showed  that  it  could  be  derived  by 
expanding  the  integrand  within  the  integral  in  a Taylor  series  about  an  appropriate  variable 
and  defining  an  average  for  that  variable  such  that  the  second  term  in  the  expansion  is  zero. 
We  are  going  to  require  averages  of  two  variables  in  order  to  perform  the  desired  integration 
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10.  B.H.  Armstrong,  ''Analysis  of  the  Curtis-Godson  Approximation  and  Radiation  Transmis- 
sion through  Inhomogeneous  Atmospheres,"  Journal  of  the  Atmospheric  Sciences,  Vol.  25, 
p.  312f  (1968). 


30 


and  will  use  an  extension  of  the  treatment  by  Armstrong  to  determine  these  average  values. 
We  will  expand  the  Integral  in  a two-dimensional  Taylor  series  about  the  two  appropriate 
variables  and  define  the  two  averages  so  that  terms  linear  in  the  differences  of  the  two  vari- 
ables are  both  zero.  We  will  also  introduce  the  average  line  strength  over  the  resolution 
interval  S at  any  given  temperature  into  the  formulation.  This  will  make  it  easier  to  relate 
the  results  to  previous  band  model  formulations  and  make  it  possible  to  use  the  results  with 
a minimum  of  changes  to  previous  methods.  The  argument  of  the  exponential  in  Equation  (15) 
becomes: 


r r S(S.E,a’)_ 
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In  order  to  obtain  a form  of  the  Curtis -Godson  approximation  similar  to  that  used  in  the 
past,  all  of  the  integrand  save  the  S(a')  term  will  be  expanded  in  a Taylor  series  as  follows: 
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Note  that  S/S  is  now  treated  as  a parameter  dependent  on  SQ,  E,  and  a'  resulting  from  the 
combined  dependences  of  S and  S.  The  three  terms  shown  are  now  used  as  an  approximation 
to  the  function  h and  substituted  into  Equation  (16)  to  yield 
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We  are  now  free  to  choose  y and  'P/T  to  minimize  the  error  in  the  Integrated  expression. 
In  particular  we  can  choose  y and  Tq/T  so  that  when  the  integration  is  performed  the  terms 
in  the  first  derivatives  vanish.  This  choice  requires  that  y and  T /T  be  given  by 
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With  the  foregoing  expressions  for  y and  Tq/T,  Equation  (18)  becomes 
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The  error  involved  in  this  approximation  can  be  evaluated  for  any  particular  path  of  interest 
by  evaluating  the  higher  order  terms  of  the  Taylor  series.  The  magnitudes  of  the  higher  or- 
der terms  compared  to  Equation  (20)  are  a measure  of  the  accuracy  of  this  approximation. 
Equation  (20)  is  a Curtis -Godson  type  of  approximation  for  the  formulation  used  here.  Note 
that  no  approximation  is  needed  for  the  emission  portion  of  Equation  (15);  the  purpose  of  in- 
troducing the  approximation  is  only  to  allow  the  order  of  integration  over  a to  be  changed,  so 
that  it  can  be  performed  last  rather  than  first  as  required  by  Equation  (15). 

It  is  also  worthwhile  to  note  that  the  variables  about  which  the  integrand  was  expanded 
are  rather  arbitrary;  we  chose  y and  Tq/T  because  the  choice  of  y leads  to  a familiar  form 
of  the  Curtis-Godson  approximation  and  also  because  it  and  Tq/T  are  natural  variables  about 
which  to  perform  this  expansion.  Armstrong  expanded  the  line  shape  function  about  the  broad- 
ening pressure  to  yield  the  Curtis-Godson  approximation  in  terms  of  an  average  pressure.  In 
cases  of  interest  here,  however,  the  large  temperature  variations  will  also  have  some  effect  on 
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the  Lorentz  line  width  so  that  expansion  about  the  width  itself  seems  more  natural.  One  can- 
not determine,  for  general  paths,  that  expansion  variable  (or  variables)  about  which  the  ex- 
pansion would  be  optimum  (that  is,  which  would  lead  to  minimum  error  in  the  path  integral). 
It  can  easily  be  shown  that  expansions  about  variables  which  are  linear  combinations  of  each 
other  are  equivalent;  but,  except  for  this,  different  expansions  are  not  easily  compared  with- 
out introducing  specific  path  examples. 


2.6  DISTRIBUTION  OF  LINE  STRENGTHS 

It  is  now  necessary  to  consider  suitable  forms  for  the  joint  distribution  f(SQ,E)  defined 
in  Equation  (14).  To  begin  with,  it  will  not  do  to  assume  SQ  independent  of  E.  Comparatively 
weak  lines  almost  always  are  weak  because  their  lower  states  are  comparatively  less  popu- 
lated, indicating  a higher  lower  state  energy.  In  simple  approximations,  in  fact,  the  matrix 
element  of  all  transitions  within  a band  are  identical;  the  variations  in  strength  of  the  lines 
in  such  a band  are  due  solely  to  variations  in  lower  state  populations  and  hence  to  changes  in 
E.  In  such  a case  S and  E are  completely  correlated.  The  2.7  urn  Ho0  band  is  made  up  of  a 
number  of  different  vibrational  transitions.  In  addition,  the  hot  bands  for  each  of  these  transi- 
tions become  important  at  the  higher  temperatures.  Since  lines  for  all  vibrational  transitions 
appear  throughout  the  2.7  pm  spectral  region,  every  Av  will  contain  lines  from  all  bands. 

Since  band  strengths  and  lower  states  for  each  of  these  bands  are  different,  one  would  not  ex- 
pect the  individual  line  strengths  and  lower  state  energies  in  a given  AD  to  be  completely  cor- 
related. To  visualize  this  effect,  we  have  in  Figure  2 plotted  S , the  strength  at  296°K,  versus 
the  lower  state  energy  E for  all  the  water  vapor  lines  in  the  40  cm  interval  between  3530 
cm"*  and  3570  cm*  contained  in  the  AFCRL  line  tabulation  of  McClatchey  et  al.  (see  Refer- 
ence 1).  There  are  about  260  lines  in  this  interval;  only  the  160  lines  having  strengths  S 

-24  -1  -2  ° 

greater  than  10  cm  /mol -cm  are  plotted.  The  density  of  points  at  a given  SQ  and  E 

value  is  representative  of  the  value  of  f(SQ,E)  at  that  point.  If  SQ  and  E were  completely  cor- 
related, the  points  would  all  lie  on  a single  curve.  They  do  not,  but  the  points  do  lie  in  a 
crudely  straight  band  running  roughly  from  upper  left  to  lower  right.  The  relation  for  the 
strength  of  an  absorption  line  at  temperature  Tq,  in  terms  of  the  matrix  element  of  the  transi  - 
tion  associated  with  that  line,  is — see  Reference  2 — H 


The  quantities  in  this  expression  that  vary  from  line  to  line  in  any  given  AD  are  the  matrix 
element  R^,  the  lower  state  energy  E,  the  wavenumber  of  the  transition  D,  and  the  rotation- 
quantum -number -dependent  quantities  m.(2j  + l)/g..  AD  is  small  enough  (of  the  order  of  10 

-1  J 1 _i 

cm  ) compared  to  the  wavenumber  of  the  bands  of  interest  here  (~3000  cm  ) so  that  the 
wavenumber  in  the  numerator  can  be  considered  sensibly  constant  for  all  the  lines  in  AD.  If 
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FIGURE  2.  LINE  STRENGTH  VERSUS  LOWER  STATE  ENERGY  FOR  WATER  VAPOR  LINES  BETWEEN  3530 

AND  3571  cm'1 
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the  product  of  the  quantities  m^(2j  + l)/g^  and  the  matrix  elements  were  identical  for  all  lines 
in  Av,  then  SQ  would  be  related  to  E by  a fixed  exponential  relation,  a straight  line  on  the 
semi -log  plot  of  Figure  2.  The  spread  of  the  points  indicates  of  course  that  the  above  men- 
tioned product  is  not  identical  for  all  lines,  but  nevertheless,  the  grouping  of  the  lines  in  a 
region  of  Figure  2 indicates  a limited  variation  in  this  quantity.  We  have  tried  representing 
the  density  of  points  in  this  figure  by  uniform  distributions  over  limited  ranges  of  E and  SQ  to 
try  to  approximate  the  grouping  of  the  points.  The  formulations  became  very  complex,  how- 
ever, and  we  decided  to  approximate  the  distribution  of  points  on  Figure  2 by  a single  straight 
line.  This  leads  to  a formulation  which  is  workable,  which  represents  at  least  the  spirit  of 
the  relationship  of  Equation  (21),  and  which  can  be  chosen  to  represent  the  mean  of  the  points 
of  Figure  2.  It  should  also  provide  a great  improvement  over  previously  made  assumptions 
in  which  all  lines  are  assumed  to  grow  with  temperature  according  to  a common  multiplicative 

function.  The  line  chosen  should  have  a slope  close  to  that  of  Equation  (21),  namely  1/kT  = 

-3  0 

4.86  x 10  cm.  A line  of  this  slope  is  shown  on  Figure  2 also.  It  can  be  seen  the  lines  are 

grouped  roughly  according  to  this  slope.  We  will  use  this  slope  without  further  justification. 

We  have  decided  to  assume  SQ  and  E completely  correlated  as  described  above;  we  must 
then  choose  a distribution  for  SQ.  Using  the  same  set  of  lines  as  used  for  Figure  2,  the  cumu- 
lative distributions  of  SQ  were  computed  and  plotted  on  Figure  3 in  semi -log  coordinates. 

This  presentation  is  similar  to  that  of  Figure  4.6  in  Reference  2.  A straight  line  in  these  co- 
ordinates represents  an  inverse  intensity  distribution  of  line  strengths.  The  points  on  this 
figure  indicate  that  the  choice  of  an  inverse  intensity  distribution  would  represent  the  data 
points  rather  well.  On  the  other  hand,  the  assumption  of  an  exponential  distribution  of  line 
strengths  leads,  in  the  case  of  transmittance  formulations  (as  opposed  to  the  transmittance 
derivative  formulation  sought  here),  to  simpler  formulations.  But  the  points  shown  on  Fig- 
ure 3 do  not  lend  themselves  at  all  to  representation  by  an  exponential  distribution.  This  fig- 
ure shows  many  more  weak  lines  than  an  exponential  distribution  would  predict.  Therefore, 
to  preserve  some  correspondence  to  reality  for  water  vapor,  we  will  use  an  inverse -strength 
distribution  of  the  form 

f(S  ) - £ (22) 

for  some  range  of  strengths  Sm  to  Sg  where  K is  a constant  of  the  distribution,  Sm  is  a mini- 
mum line  strength  and  Sg  is  the  line  strength  when  E is  the  lowest  possible  vibration -rotation 
energy  EQ.  To  summarize,  then,  we  have  assumed  the  lower  state  energy  to  be  completely 
correlated  to  the  line  strength  SQ  (SQ  is  the  line  strength  at  temperature  Tq)  and  the  strengths 
to  be  distributed  according  to  the  inverse  of  the  strength.  The  integrals  over  E in  the  pre- 
ceding expressions  will  thus  disappear. 
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Average  Strength  at  Temperature  T 

Using  the  distribution  described  above,  it  will  be  necessary  to  determine  the  average  line 
strength  as  a function  of  temperature  in  terms  of  the  parameters  of  the  distribution.  This 
will  be  necessary  not  only  because  we  will  have  to  evaluate  these  parameters  from  experi- 
mental absorption  coefficient  data  (which  are  often  presented  in  terms  of  average  line 
strength),  but  also  because  we  have  introduced  this  quantity  into  the  path  integral  in  Equa- 
tion (20). 

The  average  strength  at  any  temperature  is  given  by 

S(T)  = J S(S0,E,T)f(S0)  dSQ  (23) 

where  f(SQ)  is  to  be  taken  from  Equation  (22).  As  discussed  in  the  previous  paragraph,  SQ 
will  be  assumed  to  be  completely  correlated  with  E according  to  the  relation 

S0=Se  exp  [ -(E  - EQ)/kT0]> 


where  S is  the  strength  at  which  the  lower  state  energy  reaches  E . It  is  not  necessary  at 
e ° 

this  point  to  determine  a value  for  Sg.  Inversion  of  this  equation  yields  E in  terms  of  Sq  as 

follows: 


E(So)  = -kTolnsa»Eo 
e 


(24) 


Substitution  of  Equation  (24)  Into  the  expression  for  strength,  Equation  (13),  yields 


sis0,-f)=sj^ 


exp 


k \T 


1 

Ti 


Further  substitution  of  Equations  (22)  and  (25)  into  Equation  (23)  gives 


exp 


j>U  » 
Lk\VT/J 


KdS. 


Since  the  range  of  strengths  to  be  encountered  is  very  large,  Sm/Se  s 0,  so  that 


(25) 
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(t)  seK“>  T^-f) 


S(T)3 


SgK  exp 


In  order  to  evaluate  Equation  (20)  we  need  S/S.  From  Equations  (25)  and  (26)* 

t7t 


s/s  M vi/y 

slv  T / K IsJ 


Substituting  Equation  (27)  into  Equation  (20)  gives 


£)  ^J*,w 


(26) 


(27) 


(28) 


2.7  EXPRESSION  FOR  THE  TRANSMITTANCE  DERIVATIVE 

At  this  point,  we  are  ready  to  construct  the  final  expression  which  represents  the  trans- 
mittance derivative.  It  is  appropriate  first  to  introduce  the  usual  band  model  parameters. 
Except  for  the  variation  of  E with  SQ,  the  band  model  up  to  this  point  is  identical  with  that 
developed  by  Goody  for  an  inverse -strength  distribution  of  lines  (see  Reference  2,  pp.  130- 
136).  To  facilitate  comparisons  with  our  model,  it  is  convenient  to  use  the  same  variables 
and  nomenclature  used  by  Goody.  Thus  we  define 


and 

a 

j*S(a')da’ 

_ o 

U " 27ry6 

(29) 

a 

j*k(a')da' 

- 0 _ 

2ffy 

The  only  difference  between  these  variables  and  those  of  Goody  is  the  inclusion  of  the  factor 
(Tq/T),  the  average  Inverse  temperature  ratio,  in  y.  This  factor,  of  course,  does  not  appear 
on  Goody's  treatment  since  it  arises  from  our  relation  between  E and  SQ.  It  is  included  in  y 
to  provide  exact  correspondence  between  the  transmittance  relation  obtained  here  and 
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Equation  (4.37)  in  Goody's  transmittance.*  Note  that  k(a')  is  the  strength -to -spacing  ratio 
evaluated  at  the  temperature  at  a’,  which  appears  is  all  band  model  treatments.  Equation  (28), 
the  argument  for  the  exponential  for  Equation  (IS),  now  becomes 


M 


f-jf 

o 


7Tlr  3-^-5  st»w  = -2 


'e/  V . t,V 


u /V 
itS\\Se/ 


This  relation,  Equation  (27)  for  S in  terms  of  S and  To/T,  and  Equation  (22)  can  be  substituted 
into  Equation  (15)  for  dW/da  to  give 


e T/T 


W-laJ-0  (A\ 0 

da  ~ it T \SeJ 


2,_v  . ~2 


2u l°o 

TT^ls. 


r (a)  + a?  1 + /Av\  \ e 


=— dAj/dS 
So 


By  making  the  following  changes  of  variables, 


v = — 

y 


the  following  expression  is  obtained  • 


dW  1 _ jyl  1 1 It-1 

27”  v ” rn  - — 1 O I W 


da  S(a)  ” T y n M2  + v2 


„ T7T~ 

exp ^-s-w  0 dwdv 

1 + v 


The  ratio  S /S„  has  been  assumed  to  be  zero.  A further  substitution  of 
m e 


The  variables  u and  y are  very  simply  related  to  the  corresponding  variables  x and  p 
used  in  Reference  7 and  in  earlier  band-model  treatments  by  Simmons.  They  are  u = 4x, 
y = p/8n. 
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Q = 


2uw 


1 tv 


along  with  q(a)  = TQ/T(a)  and  q(a)  = TQ/T(a)  yields 


dW  _1  _ 

da  si[aj  ~ q y * 


f 


1 


exp  (-Q)dQdv 


> 


J 


(33) 


Note  now  that  the  expression  depends  only  on  the  normalized  optical  depth  u,  the  ratio  between 
the  local  and  the  average  inverse  temperature  q/q,  and  the  ratio  of  the  local  and  average  half - 
widths  y/y.  This  corresponding  expression  in  the  Lindquist -Simmons  approximation  depended 
on  only  the  normalized  optical  depth  and  the  halfwidth  ratio.  The  portion  of  Equation  (33)  en- 
closed in  curly  brackets  can  be  replaced  by  a normalized  incomplete  gamma  function,  y*  (see 
Abramovitz  and  StegunJ1^ , p.  260)  as  follows: 


dW  1 

S(a) 


X 

r 


AUrll) 

ffyq  \q) 


J 

-30 


(34) 


A general  solution  for  this  integral  for  a large  range  of  values  of  u,  q/q,  and  y/y  could 
not  be  found.  Extensive  numerical  evaluations  have  not  yet  been  made  but  enough  evaluations 
and  asymptotic  limits  were  found  to  obtain  expressions  which  seem  to  approximate  the  func- 
tion quite  well.  Appendix  A describes  the  evaluations  made  for  this  work. 

The  Interpretation  of  Equation  (34) 

An  expression  which  seems  to  approximate  Equation  (34)  very  well  is  an  expression  of 
the  form 


1 

Vl  +"(au)*> 

where  a and  b are  dependent  on  q/q  and  y/y.  It  is  surprising  that  this  simple  expression 
seems  to  approximate  Equation  (34)  as  well  as  it  does.  Table  1 gives  values  for  a and  b which 
seem  to  work  well  for  q/q  < 1 and  y/y  in  the  vicinity  of  0.1  to  10.  Figure  4 shows  plots  of  the 


11.  M.  Abramovitz  and  I.A.  Stegun,  Handbook  of  Mathematic  Functions,  AMS-55,  National 
Bureau  of  Standards,  Washington,  District  of  Columbia,  1964. 
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TABLE  1.  COEFFICIENTS  USED  IN  APPROXIMATING 
THE  RELATION  OF  EQUATION  (34)* 


1 

Vi  + (au)b 


where 


a = s-  for  q/q  < 1/2 

(y/yj 

a = ?isZaj.  {or  1 > q/q  > 1/2 

(y/W 


and 


b - 2(q/q)  for  q/q  < 1/2 
b = 1 for  q/q  > 1/2 


*See  Figure  4. 
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FIGURE  4,  DIMENSIONLESS  AVERAGE  EQUIVALENT  WIDTH  DERIVATIVE  VERSUS  DIMEN- 
SIONLESS OPTICAL  DEPTH  FOR  TWO  HALFWIDTH  RATIOS,  FOR  VARIOUS  TEMPERATURE 

RATIOS 
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approximation  for  y/y  = 1.0  and  2.0,  for  various  values  of  q/q.  Figure  4(b)  also  shows 

several  computed  values  of  the  complete  function  of  Equation  (34).  Since  these  values  were 
calculated  numerically  by  hand,  they  are  of  unknown  accuracy.  Agreement  with  the  approxi- 
mation seems  to  be  of  the  order  of  5%.  Appendix  A discusses  the  means  used  to  obtain  the 
approximations.  It  would  have  been  desirable  to  perform  extensive  tests  of  this  approxima- 
tion compared  to  Equation  (34),  but  time  and  funds  were  not  available  to  do  so.  The  result  is 
important  enough  to  be  presented  as  we  have,  with  minimal  confirmation,  since  (as  discussed 
in  Section  4)  it  forms  the  basis  of  a much  improved  ability  to  predict  plume  radiation  through 
an  atmospheric  path. 


The  range  of  parameters  which  are  important  for  plume  calculation  is  as  follows.  The 
parameter  q/q  is  the  ratio  between  the  inverse  of  the  ocal  temperature  to  the  average  in- 
verse temperature.  For  the  experiments  presented  in  this  work 


1/1200 
q”  17300 


0.25 


This  is  the  reason  for  the  use  of  q/q  = 0.25  in  Figure  4.  For  a plume  submerged  in  an  atmo- 
sphere at  250°K,  q/q  could  be  expected  to  vary  between,  say, 


1/600 

17250 


- 0.42  and 


1/1600 

1/250 


0.16 


Values  of  q/q  > 1 correspond  to  targets  cooler  than  ambient,  not  an  interesting  case  for  plume 
studies.  As  the  path  progresses  from  the  near  boundary  of  the  target  to  the  far  boundary,  the 
value  of  q/q  approaches  unity.  In  the  experiments  described  in  Section  3,  the  value  of  q/q 
went  from  0.25  at  the  near  boundary  at  all  frequencies  and  increased  to  about  0.3  at  the  far 
boundary  in  the  band  center.  At  the  edges  of  the  band,  where  there  is  a great  difference  be- 
tween the  values  for  S/6  at  ambient  temperature  and  at  1200°K,  the  value  of  q/q  at  the  far 
boundary  became  somewhat  greater  than  0.4. 


y/y  has  a rather  special  interpretation  in  this  model.  If  this  ratio  is  interpreted  strictly 
as  a ratio  of  halfwidths,  it  should  be  a function  only  of  the  local  pressure  and  temperature 
and  the  average  pressure  and  temperature.  But  earlier  in  the  development  the  variables  K 
and  y were  lumped  together  into  the  variable  y.  This  variable  controls  the  approach  of  the 
transmittance  curve-of -growth  to  the  square -root  region;  to  ensure  correspondence  with 
reality,  its  value  is  either  taken  from  experimental  data  or  computed  using  the  special 
strength  weighting  given  by  Goody  (Reference  3,  Equation  (4.41)).  In  the  same  fashion,  the 
parameter  y/y  controls  the  approach  of  the  transmittance  derivative  into  the  region  that  cor- 
responds to  the  square -root -region  for  transmittance.  Hence,  the  parameter  y/y  must  also 
be  taken  !-;ra  the  variable  y.  Using  the  definition  given  in  Equation  (29), 


y = 


-Bl 

5q 


y = % 
y 6q 
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Note  that  K is  a fixed  quantity  for  the  whole  path  since  it  is  a parameter  of  the  distribution 
for  SQ.  Also  6 is  invariant  over  the  path;  since  the  model  assumes  that  all  lines  are  present 
at  all  path  conditions,  only  their  strengths  vary.  Thus  y/y  should  be  taken  as 

t_x.a 

Y fq 

This  will  ensure  proper  correspondence  with  the  square -root  region.  This  ratio  corresponds 
directly  to  the  ratio  p/8  in  the  Lindquist -Simmons  approximation  except  for  the  temperature 
factor. 

Youngt^]  indicates  that  values  of  y/y  ranging  from  10*7  to  200  can  be  expected  in  plume 
calculations.  Thus  the  investigation  of  Equation  (34)  given  here  is  nowhere  near  complete. 

We  have  investigated  the  behavior  of  Equation  (34)  only  for  values  of  y/y  of  1 and  2.  A great 
deal  of  effort  remains  to  be  accomplished  before  the  behavior  of  Equation  (34)  can  be  said  to 
be  well  known. 

_ _ -1/2 

A strange  feature  of  Equation  (34)  is  that,  for  q/q  < 1/2,  d,W/da  1/S  does  not  go  as  u 
for  larger  u,  but  rather  goes  as  u"*^.  This  is  apparent  in  the  approximations  given  in 
Table  1 and  in  the  values  calculated.  This  behavior  is  puzzling  and  contrary  to  what  one 
would  expect.  It  means  that  so  many  weak  lines  become  strong  that  square -root  growth  can 
not  be  achieved.  This  is  another  facet  of  Equation  (34)  that  requires  further  investigation. 

From  the  evaluations  and  analysis  that  were  done  of  Equation  (34),  it  is  apparent  that  the 
character  of  the  function  changes  greatly  at  q/q  = 1/2.  Unfortunately,  we  have  not  yet  done 
any  careful  evaluations  at  q/q  = 1/2  in  order  to  uncover  exactly  why  the  function  changes  be- 
havior at  this  value.  Further  investigation  of  this  problem  should  be  undertaken  before  ex- 
tensive use  is  made  of  Equation  (34). 

2.8  USE  Or  THE  BAND  MODEL 

We  have  imcorporated  the  bandmodel  represented  of  the  preceding  sections  into  our  cal- 
culation procedures.  Much  of  the  procedure  is  similar  to  that  used  previously;  only  those 
portions  dealing  with  the  evaluation  of  the  transmittance  derivative  have  been  modified.  This 
routine  performs  a numerical  evaluation  of  Equation  (5),  using  Equation  (11)  to  obtain  values 
for  dr/da.  Since  Equation  (11)  and  all  subsequent  equations  through  Equation  (34)  assume  a 
fixed  frequency  Interval,  the  derivative  in  Equation  (11)  corresponds  exactly  to  the  partial 
derivative  given  in  Equation  (5).  By  means  of  Equation  (11)  it  is  easy  to  evaluate  dr/ da  once 
Equation  (34)  has  been  evaluated: 


i2.  S.J.  Young,  Band  Model  Calculations  of  Atmospheric  Transmittance  for  Hot  Gas  Line 
Emission  Sources,  SAMSO-TR-74 -248,  The  Aerospace  Corporation,  El  Segundo, 
California,  December  1974, 
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d?_  r dW  -SfldW' 
da  “ 5 da  “ T‘  6[_f  da  _ 


(35) 


where  the  quantity  in  the  brackets  is  the  quantity  given  by  Equation  (34)  and  S/6  is  the  strength  - 
to -spacing  ratio.  Thus  at  any  position  a in  the  path,  an  expression  for  dr/da  can  be  obtained 
if  band  model  parameters  along  the  whole  path  are  available.  The  band  model  parameters 
required  are  k = S/6  and  y as  a function  of  temperature  (and  frequency). 


The  computational  procedure  used  here  operates  as  follows.  Given  a non-isothermal, 
non -iso baric  path  of  known  temperature  and  composition,  it  is  first  divided  into  a number  of 
segments,  each  of  which  can  be  considered  sensibly  isothermal.  The  number  of  integrals  into 
which  a path  must  be  broken  for  accurate  representation,  which  has  been  investigated  by 
Young^1^  , depends  e*t  the  path  composition.  In  our  routine,  this  break-up  can  be  performed 
either  automatically  or  manually. 

Once  the  path  is  broken  up  into  n segments,  the  following  variables  are  determined  for 
each: 


*1  - 1'o(Tl)-pbi 
A(ak)t  = klTp-pej-l. 


(36) 


where  k is  the  average  strength -to -spacing  ratio  for  temperature  T^  pe_  is  the  effective  ab- 
sorbing pressure  of  the  active  gas  (either  COg  or  HgO)  in  atmospheres  at  STP,  and  f . is  the 
length  of  the  segment.  The  variable  y is  dimensionless;  furthermore,  the  reference  pressure 
for  the  value  of  y in  y is  one  atmosphere  — so  that  the  pressure  in  atmospheres  can  be  used 
to  preserve  the  dimensionless  character  of  y.  To  retain  identity  with  previous  work,  we  have 
used  optical  path  length  in  atmosphere -cm.  Thus, 


"ei  " Hia 


27305 
ri 


where  pJa  is  the  actual  partial  pressure  of  the  active  gas  in  segment  i.  This  temperature  ad- 
justment on  the  pressure  is  required  to  insure  correspondence  between  the  optical  path  length 
in  atm-cm  and  the  number  of  absorbing  molecules  in  a unit  cross-section  of  path.  The  above 
definition  for  pgl  allows  easy  comparison  with  the  AFCRL  line  tabulation^ . Since  the  vari- 
able k is  the  strength -to -spacing  ratio,  it  will  have  the  units  of  the  reciprocal  of  the  optical 
path.  (Strength  has  the  units  of  wavenumber  per  optical  path  unit;  spacing  is  in  units  of  wave- 
number.)  Thus  the  products  A(ak)^  are  ’-nitless. 

P^l>  which  is  the  broadening  pressure,  is  given  by 

"bl  ■ "at 


where  the  sum  over  j is  a sum  over  all  the  other  species  present  and  is  the  broadening 
coefficient  for  species  J with  the  active  species  a.  p^  is  the  pressure  of  species  j in  path  seg- 
ment i.  Since  few  accurate  values  for  such  broadening  coefficients  are  actually  available,  we 
assume  that  all  other  species  broaden  like  air  and  use  the  values  for  HgO  and  C02  valid  for 
air -broadening 

o(H20,  air)  = 0.2 
o(C02,  air)  = 0.77 

yo(Ti>  is  the  2nd  band-model  parameter  at  one  atmosphere  pressure  and  temperature  ob- 
tained from  band-model  tabulations. 


From  the  optical  depth,  we  take  an  average  value  y and  define  a dimensionless  quantity  u: 


The  average  inverse  temperature  is  also  calculated*. 

% 

q,  • ^ (38) 

where  q^  = -r. 


The  transmittance  derivative  is  now  computed  at  both  euds  of  every  path  segment.  This 
is  accomplished  from  Equation  (35): 


dr  _ - S /l  dw\ 
da  " T ' da  j 


= r(a)k(a)f 


At  the  near  edge  of  segment  l, 


46 


unear,i  “ ui-l 


a\  I =jLjL 

tt/neair,i  ?!-!  Ql-l 


near,!  i-1 


and 


si  = VA'k-r  =)  1 

'near,i  lLl  a/near,i  Q/near,iJ 


For  the  far  boundary  of  each  segment,  the  corresponding  values  are 
ufar,i  = ui 


__h% 

far,i  yi  Qi 


and 


dr 

da 


) -T.k.f  L?)  ;&)  1 

'far,!  1 L 1 «/far,i  vfar.ij 


The  transmittances  r used  in  these  expressions  were  obtained  from  the  curve -of -growth  for 
this  band  model.  This  expression  has  been  shown  to  be  identical  to  that  given  by  Goody  for 
an  inverse  strength  distribution  of  lines  (see  Reference  3,  Equation  4.37). 


?t  * fc!v[e  ‘V“i 


-u 
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) + 2uj  e (yuj)  + IjfUj))  “ 1 


■] 


(39) 


IQ  and  Ij  are  Bessel  functions  of  the  first  kind  with  imaginary  arguments.  Values  for  IQ  and 
Ij,  were  taken  from  [11] . 

The  effective  transmittance  derivative  for  a given  segment  is  taken  to  be  the  average  of 
the  near  and  far  values  for  that  segment: 


‘vsK’lfs),  *s)  ,k 

'i  L far,l  'nearjiJ 


In  order  to  ensure  consistency  between  these  derivatives  and  the  r values  used  for  farther 
segments,  the  transmittance  at  the  far  boundary  of  each  segment  is  modified,  after  the  trans- 
mittance derivative  calculations,  to  be  equal  to  the  transmittance  at  the  near  boundary  of  that 
segment  minus  the  transmittance  change  determined  for  that  segment: 
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h = Vi  ‘ A?i 

Thus,  only  Tj  remains  unchanged  from  its  value  obtained  directly  by  Equation  (39). 

The  values  of  Ar^  are  then  used  directly  in  a form  analogous  to  Equation  (6)  to  obtain  the 
radiance 

h s^L$Ti),An 

These  procedures  provide  means  by  which  the  band  model  developed  here  is  used  to  compute 
radiances  of  non -isothermal  paths. 

It  is  not  clear  how  the  number  of  segments  into  which  the  path  is  divided  affects  the  re- 
sults of  the  calculation.  This  is  the  most  Important  aspect  yet  to  be  analyzed  in  the  calcula- 
tions themselves.  It  is  not  sufficient  to  divide  the  path  into  segments  of  the  order  of  the 
temperatu  re -scale -length  because  more  variables  than  temperature  are  involved.  Since 
dr/da  is  not  constant  within  a segment  even  though  the  segment  is  isothermal,  it  appears  that 
a careful  treatment  of  the  number  of  segments  is  required.  Young  states  that,  for  near  homo- 
geneous paths,  the  Lindquist -Simmons  approximation  requires  a larger  number  of  segments. 
The  process  described  here  provides  a treatment  for  the  first  path  element  which  is  identical 
to  that  provided  by  an  isothermal  treatment  and,  hence,  should  do  a good  job  for  the  first  seg- 
ment. On  the  other  hand,  when  the  second  and  third  and  subsequent  segments  are  of  similar 
composition  to  the  first,  the  treatment  of  their  elements  by  the  present  procedure  is  not  so 
easily  examined.  These  unknown  properties  of  the  current  calculation  procedures  should  be 
carefully  examined. 


3 

MEASUREMENTS  OF  THE  HOT  WATER  VAPOR  AND  C02  RADIANCE 
AS  VIEWED  THROUGH  COOL  WATER  VAPOR  AND  C02 

tto  definitive  set  of  measurements  existed  which  accurately  characterizes  the  amount  of 

correlation  actually  present  between  hot  and  cold  water  vapor  and  CO«  lines  in  the  bands  of 

fll 3 * *  & 

interest  for  plume  measurements.  Burch1  J made  several  measurements  which  served  only 
to  demonstrate  the  effect.  We  are  aware  of  no  other  controlled  measurements  of  this  corre- 
lation. Therefore  a portion  of  the  present  program  was  devoted  to  the  measurement  of  the 
radiance  of  hot  water  and  CO2  as  viewed  through  cool  mixtures  of  the  same  gases.  These 
measurements  can  be  looked  upon  as  measurements  of  the  degree  to  which  the  product  rule 
discussed  in  Section  2 holds,  or  as  measurements  of  the  transmittance  derivative,  or  they  can 
be  compared  directly  to  band  model  techniques  of  various  types. 


48 


We  originally  planned  to  measure  the  transmittance  derivative  directly.  The  measure- 
ments were  to  involve  the  radiance  of  a very  short  isothermal  hot  cell  viewed  through  a long 
cold  cell.  This  situation  would  closely  approximate  a measurement  of  the  integrand  of  Equa- 
tion (5).  Actually,  however,  a hot  cell  60  cm  long  was  required  to  obtain  a large  enough  radi- 
ance signal  to  be  measured  accurately.  In  the  frequency  bands  of  interest,  at  conditions  appli- 
cable to  plume  studies,  this  cell  length  produces  optical  paths  extending  into  the  square  roof 
region  at  many  wavelengths.  The  measurements  are  valid  measurements  of  the  effect  of  line 
correlation,  but  band  model  calculations  are  required  to  extract  transmittance  derivatives 
from  them. 

The  measurements  are  of  emission -absorption  spectra  of  a two -element  optical  path. 

The  two  elements  consist  of  (a)  a heated  cell  filled  with  a mixture  of  COg,  HgO,  and  nitrogen  j 

(and,  in  a few  cases,  other  gases  to  simulate  the  combustion  products  present  in  a rocket  ex- 
haust) and  (b)  a long  cell  at  ambient  temperature  filled  with  mixtures  of  HgO,  COg,  and  nitro- 
gen to  simulate  atmospheric  paths.  The  cool  cell  is  a multiple  pass  evacuable  cell  capable  of 
operation  at  a variable  number  of  path  lengths.  The  two  paths  were  arranged  in  series  so  that 

< 

the  instruments  viewed  the  hot  cell  through  the  cool  cell. 

i 

The  detection  apparatus  consisted  of  a commercial  grating  monochromator  and  detector  ! 

set  up  to  give  spectral  resolution  appropriate  to  band  model  methods.  A heated  blackbody 
source  was  used  for  the  absorption  measurements.  The  instrument  was  always  operated  as  a 
dual -beam  instrument  in  that  the  measurements  were  always  referred  to  identical  measure-  j 

ments  made  through  a nitrogen-flushed  reference  beam.  The  latter  was  identical  in  configura- 
tion to  the  sample  beam  except  that  a long  cool  cell  was  present  only  in  the  sample  beam. 

Every  spectral  scan  was  taken  as  the  ratio  of  the  energy  in  the  sample  beam  to  the  energy 
received  in  the  reference  beam.  This  served  to  compensate  for  source  variations  and  for 
residual  ambient  atmospheric  absorption  not  eliminated  in  the  system  by  other  means. 

3.1  DETAILED  EXPERIMENT  DESCRIPTION 

Hot -through -cold  radiative  transfer  measurements  were  performed  with  the  experimental 
system  shown  schematically  in  Figure  5.  The  hot  gas  sample  was  placed  in  a 60  cm  alumina 
cell  with  sapphire  windows  (H/C)  and  maintained  at  the  desired  temperature  by  an  electrically 
heated  oven  (0).  The  long  ambient -temperature  absorbing  path  was  obtained  using  a multiple 
pass  White  cell  (W/C)  capable  of  path  lengths  up  to  200  meter The  White  cell  could  be 
put  in  or  out  of  the  optical  path  by  means  of  a sliding  carriage  on  which  mirrors  Fj,  Fg,  and 
F3  were  mounted.  Spherical  mirrors  Sj  and  Sg  were  used  in  conjunction  with  flats  Fj,  F^ 
and  Fg,  Fg  respectively  to  match  the  F -number  of  the  White  cell  with  that  of  the  optical 
system. 

13.  G.H.  Lindquist,  et  al..  Investigations  of  Chemical  Laser  Processes,  Report  No.  191300-1-P, 
Environmental  Research  Institute  of  Michigan,  Ann  Arbor,  Michigan,  February  1973. 
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Reference 


USED  FOR  RADIATIVE  TRANSFER  MEASUREMENTS 


i 

Radiation  from  the  blackbody  source  (BB)  was  split  into  two  beams:  a reference  beam  | 

chopped  at  420  Hz  and  a sample  beam  chopped  at  780  Hz.  After  passing  through  separate  legs  f 

of  the  system,  these  beams  were  recombined  at  a Perkin-Elmer  Model  89G  spectrometer  | 

(Sp),  detected  with  a single  nitrogen-cooled  detector  (D),  and  synchronously  amplified  using  j 

two  lock-in  amplifiers  operating  at  420  and  780  Hz  respectively.  Signals  from  these  amplifiers  j 
were  digitized  and  ratioed  using  a Hewlett-Packard  digital  voltmeter.  They  were  then  punched  \ 

I 

on  computer  cards  for  data  reduction.  Background  atmospheric  absorption  in  the  paths  tra-  \ 

versed  by  both  beams  was  minimized  by  flushing  the  entire  system  with  dry  nitrogen  obtained 
from  the  boil -off  of  a liquid  nitrogen  tank.  Figure  6(a)  shows  an  overall  view  of  the  measure-  ! 

ment  system. 

Both  emission  and  absorption  scans  were  performed  with  the  same  system.  For  an  ab- 
sorption measurement  all  shutters  (Bj,  B2,  Bj,  B^)  were  left  open  and  chopper  Cj  was  oper- 
ated with  chopper  C2  stationary.  In  this  configuration  any  radiation  emitted  by  the  hot  gas 
was  unchopped  and  would  not  be  detected.  For  an  emission  scan,  shutter  Bj  was  closed  and 
chopper  C2  operated  with  chopper  stationary.  This  eliminated  blackbody  radiation  from 
the  sample  side  of  the  system  and  chopped  all  emitted  radiation.  As  will  be  seen  below,  both 
absorption  and  emission  scans  are  required  to  characterize  the  emission  and  absorption  of 
the  path  in  the  system. 

Measurements  were  made  both  on  the  2.7  pm  water  vapor-C02  band  and  the  4.3  pm  C02 
band.  The  system  configuration  was  somewhat  different  for  each  of  the  two  spectral  bands. 

Table  2 summarizes  these  differences.  \ 

The  line  correlation  effect  to  be  observed  is  temperature  dependent;  hence,  measure- 
ments at  a variety  of  hot  cell  temperatures  would  be  desirable.  However,  a temperature  near 
1200°K  was  required  to  obtain  the  necessary  signal  levels.  This  temperature  was  very  close 
to  the  highest  which  the  cell  could  sustain  without  softening  of  the  window  seals.  Because  of 
these  two  considerations  the  measurements  were  limited  to  1200°K. 

To  ensure  even  heat  distribution,  a 60  cm  alumina  hot  cell  was  enclosed  in  a Monel  tube 
which,  in  turn,  was  surrounded  by  seven  individually  controllable  furnace  elements.  The 
j alumina  hot  cell  is  shown  in  Figure  6(b).  The  thermocouples  placed  along  the  Monel  tube 

I monitored  the  cell  temperature.  The  furnace  elements  were  adjusted  manually  to  achieve 

isothermal  conditions  along  the  hot  cell. 

3.2  HOT  CELL  CONSTRUCTION 

The  construction  of  the  hot  cell  posed  substantial  problems.  A fused  quartz  cell  was 
available  for  this  program  which  would  sustain  the  high  temperatures  required.  It  was  used 
for  preliminary  measurements  at  2.7  pm,  but  did  not  transmit  beyond  about  3 pm.  Sapphire 
appeared  to  be  the  only  window  material  capable  of  both  withstanding  the  high  temperatures 
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TABLE  2.  SPECTROMETER  PARAMETERS 

2.7  am  Measurements 

4.3  am  Measurements 

Detector 

Cooled  PbS 

Cooled  InSb 

Grating 

64  mm  x 64  mm  Bausch  & Lomb  replica  100  grooves/mm,  5.18  am  blaze 

Slit  Width 

500  am 

150  am 

Spectral 

Resolution 

14  cm”* 

3.2  cm"* 
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(a)  Overall 

FIGURE  6.  PHOTOGRAPHS  OF  EXPERIMENTAL  APPARATUS 


(b)  Alumina  Hot  Cell 

StAPHS  OF  EXPERIMENTAL  APPARATUS 


and  of  transmitting  beyond  3 jim.  However,  in  previous  programs  we  had  never  succeeded  in 
maintaining  a seal  between  a sapphire  window  and  a cell  material  at  temperatures  near  1200°K. 

At  the  beginning  of  the  program  we  obtained  information  about  a sealing  process,  from 
National  Beryllia  Corporation,  which  was  demonstrated  to  maintain  a vacuum-tight  seal  be- 
tween sapphire  and  alumina.  It  was  our  plan  to  purchase  a custom-made  cell  from  National 
Beryllia  consisting  of  a 60 -cm  long,  1-inch  diameter  alumina  tube  with  sapphire  windows 
sealed  to  it  with  the  National  Beryllia  process.  However,  a furnace  capable  of  firing  such  a 
long  cell  was  not  available'. 

A number  of  unsuccessful  attempts  were  made  to  seal  the  windows  to  the  cell  by  obtaining 
optical  contact  between  sapphire  windows  and  the  ends  of  an  alumina  tube.  We  did  succeed  in 
making  seals  at  room  temperature  with  such  a technique;  but  when  the  cell  was  heated,  differen- 
tial tharr'~\  expansion  of  the  alumina  tube  always  caused  the  seal  to  fail  over  at  least  a por- 
tion of  the  seal. 

A sample  of  the  frit  used  in  the  National  Beryllia  sealing  process  had  been  obtained 
earlier  along  with  instructions  for  its  use.  After  several  unsuccessful  attempts  to  find  fur- 
naces capable  of  performing  the  required  firings,  we  finally  built  a simple  propane -gas -fired 
brick  kiln  of  the  appropriate  dimensions  to  enclose  the  cell.  The  sapphire  windows  were  suc- 
cessfully fired  onto  the  cell  with  this  kiln  using  the  National  Beryllia  frit.  An  NBS  -calibrated 
optical  pyrometer  was  used  to  measure  the  critical  firing  temperatures.  The  windows  were 
held  to  the  alumina  tube  during  firing  by  applying  a very  slight  vacuum  to  the  fill  tube  which 
had  been  previously  fired  to  the  alumina  tube. 

This  process  was  ultimately  successful  in  obtaining  a rugged  gas-tight  cell  with  sapphire 
windows  capable  of  withstanding  temperatures  of  1200°K  - see  Figure  6(b).  We  consider  the 
successful  construction  of  this  cell  to  be  a significant  technological  accomplishment,  particu- 
larly in  view  of  the  many  previously  unsuccessful  attempts  to  construct  such  a cell,  both  on 
this  and  several  previous  contracts. 

3.3  COLD  CELL  CONSTRUCTION 

The  cold  cell  wa6  an  evacuable  5 meter -long  multiple -pass  White  cell  constructed  for 
this  effort.  In  a cell  of  this  type  the  energy  entering  the  cell  at  the  entrance  slit,  where  the 
incoming  energy  is  focussed,  is  reimaged  by  a spherical  mirror  at  the  far  end  of  the  cell  onto 
a second  spherical  mirror  at  the  near  end  of  the  cell.  This  mirror  directs  the  energy  onto  a 
second  spherical  mirror  at  the  far  end  of  the  cell  which,  in  turn,  redirects  the  energy  back  to 
another  position  on  the  front  mirror.  This  process  continues,  alternating  between  the  two 
spherical  mirrors  at  the  far  end  of  the  cell  until  the  energy  focussed  back  on  the  front  mirror 
coincides  with  the  position  of  the  exit  slit.  Figure  7 is  a schematic  of  the  cell  showing  the 
light  path. 
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FIGURE  7.  5 METER  "WHITE"  CE 


The  cell  was  made  of  thick  wall  hard  aluminum  tubing.  The  mirror  mounts  were  designed 
to  be  as  immune  as  possible  to  vibration.  Vibration  effects  become  significant  at  total  path 
lengths  near  200  m (40  passes).  Most  measurements  were  made  at  a 100  m path  length  (20 
passes). 

3.4  GAS  HANDLING  SYSTEM 

Since  measurements  were  to  be  made  on  mixtures  of  a number  of  gases,  a manifold  sys- 
tem was  constructed  whereby  any  desired  mixtures  of  gases  could  be  applied  to  the  cold  cell 
and  the  hot  cell  individually.  The  basic  pressure  measurements  for  the  filling  process  were 
made  uft*ng  a capacitance  manometer  capable  of  reading  from  0.01  to  100  mm  of  Hg.  Higher 
pressures  were  measured  with  a Wallace  and  Tieraan  high  accuracy  pressure  gauge  in  the 
case  of  the  hot  cell,  and  a mercury  manometer  in  the  case  of  the  cold  cell.  The  water  vapor 
was  obtained  from  a constant  temperature  bath  of  distilled  water,  which  also  served  as  a 
check  on  the  accuracy  of  the  capacitance  manometer. 

The  filling  process  for  either  cell  was  carried  out  by  first  filling  the  cell  with  the  gas  to 
be  present  at  the  lowest  partial  pressure.  Subsequent  gases  were  added  in  the  oi  ier  of  in- 
creasing desired  partial  pressures.  This  procedure,  in  addition  to  the  excellent  reproducibil- 
ity and  accuracy  of  the  capacitance  manometer,  made  accurate  fills  possible. 

Adsorption  effects  were  noted  for  water  vapor,  particularly  in  the  hot  cell.  The  account 
for  these  effects,  the  water  vapor  pressure  was  allowed  to  equilibrate  until  the  pressure  read 
on  the  capacitance  manometer  agreed  with  the  equilibrium  vapor  pressure  of  the  water  bath. 
Additional  gases  were  added  to  the  hot  cell  only  after  this  equlibration process  was  complete. 

No  adsorption  effects  were  observed  for  any  of  the  other  gases  used.  As  a result  of  this  ob- 
servation, the  assumption  was  made  that  the  presence  of  the  other  gases  did  not  alter  the 
equilibrium  between  the  water  vapor  in  the  cell  and  the  water  adsorbed  on  the  cell  walls. 

Based  on  the  published  accuracy  of  the  capacitance  manometer  system  and  on  calibrations 
made  with  constant  temperature  baths,  and  considering  our  attempts  to  assure  cell  equilibrium, 
it  is  our  opinion  that  the  cell  fills  are  accurate  to  within  several  percent. 

3.5  THE  EFFECT  OF  THE  CELL  WINDOWS  ON  THE  MEASUREMENTS 

Since  the  hot  gas  cell  is  maintained  at  an  elevated  temperature  the  emissions  of  the  cell 
windows,  as  well  as  their  absorption  and  reflection,  must  be  accounted  for  if  a true  gas  emis- 
sion is  to  be  determined.  For  an  optically  thick  material  it  may  be  shown  that  the  emission 
at  a given  wavenumber  v is  given  by 

L~  = l!o~  (40) 

v v v 

where  L~  is  the  spectral  radiance  of  the  material,  a its  absorptance  and  L^  the  spectral  ra- 
diance of  a blackbody  at  the  same  temperature.  For  a thin  optically  transparent  material 
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like  the  windows  of  the  cell  used  here,  multiple  reflections  must  be  taken  Into  account.  In 
this  case  the  spectral  radiance  is  given  by 


I~  = l![p<o  •<**)  +«]/[!  -P2(l  -a)2) 


L~  = -p(l  -o)J  (41) 

Here  a and  p are  the  hot  cell  window  abeorptance  and  reflectance.  For  simplicity,  the  spec- 
tral dependences  of  a and  p are  ignored. 

At  1200°K,  for  the  sapphire  windows  used,  a ~ 0.18  and  p ~ 0.07  in  the  4.3  pm  band,  and 
a - 0.02  and  p ~ 0.07  at  2.7  pm.  Although  this  implies  an  error  of 

t f 1.061  at  4.3  pm 


1 • [1.074  at  2.7  pm 

if  multiple  reflections  are  ignored,  it  may  be  shown  that  the  correction  to  the  hot  gas  radiance 
or  hot-through-cold  radiance  is  given  by 


‘"VI1  'll 


Therefore,  neglecting  multiple  reflections  introduces  an  error  of  ~1%  at  4.3  pm  and  ~0.1% 
at  2.7  pm.  Because  of  the  added  complexity  in  accounting  for  the  multiple  reflections,  and 
because  we  do  not  expect  accuracies  better  than  about  5%  from  the  system  as  a whole,  mul- 
tiple reflections  have  been  ignored. 

As  a result  of  the  complications  introduced  by  the  alumina  cell,  measured  radiances  do 
not  give  the  derived  hot  through  cold  radiance  directly.  This  may  be  appreciated  by  writing 
an  explicit  expression  for  the  observed  radiance  in  terms  of  the  optical  parameters  of  the  ex- 
perimental system.  Using  the  notation  of  Figure  8,  and  assuming  that  the  hot  cell  is  isother- 
mal throughout,  we  can  set 

L*  > L*  L* 
w g 

Again,  the  spectral  dependence  is  suppressed  for  simplicity.  The  measured  radiance  Y of  the 
hot  cell  and  its  gas  sample  as  seen  through  the  long  cold  path  is,  therefore, 

Y - rr-Vw/cTg(1  •p)2<1  "“w*1  • (1 

LBB' 


"»  (1  ‘ p)2(1  -°w>] 


FIGURE  8.  NOTATION  USED  IN  DERIVING  THE  HOT  THROUGH  COLD  RADIANCE  MEASUREMENT  EXPRESSION 


again  ignoring  multiple  reflections  within  the  cell  windows.  This  equation  was  obtained  by 
summing  the  contributions  from  the  rear  window,  the  gas  in  the  cell,  and  the  front  window, 
and  then  separating  out  those  terms  containing  a . 

D 

The  desired  output  is  the  radiance  of  the  hot  gas  as  viewed  through  the  White  cell.  In 
terms  of  the  variables  of  Figure  8,  this  is  given  by  inversion  of  Equation  (42)  as 


L*ot  r = 
g g 


y-^u-pV,Vc[‘*(‘-*>2(‘-vI 

bB* j 


_LBB' 


Since  the  bracketed  term  in  the  denominator  is  on  the  order  of  1.01  with  sapphire  windows  at 
the  experimental  conditions  used,  it  may  be  taken  equal  to  unity  with  negligible  error.  Addi- 
tional measurements  are  necessary  in  order  to  account  for  the  terms  (1  - pVjL-T  / 1 + 

0 1 22  ww/CL 

(1  - p)  (1  - a ) and  (1  - p)  (1  - a ) . The  following  measurements  are  convenient  for  this: 

an  emission  measurement  with  the  White  cell  and  hot  cell  empty  (a  = 0,  r = 1),  which  we 

D O 

call  a "cell"  scan,  given  by 

CeU  = -i-(l  - phv[l  * (1  - P)2d  - %>]  Tw/c 
lbb' 

and  another  emission  measurement  similar  to  the  "cell"  scan  except  with  the  blackbody  shut- 
ter B1  open  to  the  sample  side  of  the  system  ("Beam"  scan) 


Beam 


^•(1  -p)4(l  -«w)2 


+ cell 


In  terms  of  these  measurements,  the  hot  through  cold  radiance  is 


I-  Y - ceU  t 

LgTg  = L “gTB  = Beam  - cell 


The  blackbody  radiance  L*  B may  be  calculated  using  the  Planck  function  once  its  temper- 
ature is  determined  with  an  optical  pyrometer  and  the  hot  cell  window  reflectance  p is  calcu- 
lated from  refractive  index  data.  Consequently,  the  only  unknown  in  this  equation  is  the  trans- 
mission of  the  cold  gas  ~'"nple  in  the  White  cell  (r^).  This  parameter  may  be  evaluated  using 
the  "Beam"  scan  above  and  an  identical  scan  made  with  the  White  cell  filled.  In  this  case, 

- Beam  (White  cell  filled)  /45\ 

Tg  “ Beam  (White  ceil  empty) 

where  it  is  assumed  that,  because  there  is  no  correlation  between  the  transmission  of  the 
White  cell  itself  and  the  cold  gas  sample,  the  total  transmission  of  this  system  is 


i . .vcf .. ..  .*.i 
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Determination  of  hot  through  cold  radiance  therefore  requires  four  individual  scans: 

Y = emission  measurement  of  hot  cell  and  hot  gas  viewed  through  a long  cold  path 
Cell  = emission  measurement  of  empty  hot  cell 

Beam  (White  cell  empty)  = emission  measurement  of  empty  hot  cell  with  blackbody  radia- 
tion background,  White  cell  empty 

Beam  (White  cell  filled)  = emission  measurement  of  empty  hot  cell  with  blackbody  radia- 
tion background,  White  cell  filled 


3.6  SYSTEM  CHECKS 

During  data  taking  the  oven  temperature  was  continuously  monitored  with  ten  thermo- 
couples placed  along  the  length  of  the  alumina  cell.  The  internal  consistency  of  the  data  can 
be  tested  by  calculating  an  apparent  temperature  for  the  emitting  gas  and  comparing  this  value 
with  the  oven  thermocouple  temperatures.  The  radiance  and  the  absorption  of  the  hot  cell 
alone  (cold  cell  evacuated)  as  measured  by  the  system  is  used  in  the  inverse  Planck  function 
to  determine  an  apparent  temperature  T*. 


2~3 

2hc  v a 


*+l 


(46) 


where  is  the  measured  gas  spectral  radiance  obtained  from  the  scans  of  Equation  (42)  with 
the  cool  cell  empty,  h is  Planck’s  constant,  c the  speed  of  light,  v the  spectral  position  in 
cm"1,  and  a the  absorption  of  the  gas  at  wavenumber  v. 

Two  additional  scans  are  required  to  evaluate  a independent  of  L : (a)  an  absorption 

O w 

scan  (ABSO)  with  the  hot  cell  filled  and  the  cold  cell  empty  and  (b)  an  absorption  background 

scan  (ABSB)  identical  to  the  absorption  scan  but  with  both  cells  empty.  From  these  two  scans  the 

hot  gas  absorption, a can  be  obtained: 

© 

„ _ ABSB-ABSO  /47\ 

"g ~ ABSB  (47) 


and  consequently  T*  evaluated.  The  discrepancy  between  T*  and  the  oven  temperature  deter 
mined  from  the  thermocouple  readings  is  indicative  of  errors  in  the  data  taking  or  reduction 
procedures.  These  checks  were  performed  on  all  measurements. 


A second  check  was  performed  several  times  during  the  program  to  verify  the  pyrometric 
temperature  reading  of  the  blackbody.  In  this  case  the  hot  cell  is  filled  with  a gas  sample 
which  is  opaque  in  the  spectral  region  of  interest  (at  = 1).  An  emission  scan  made  with  this 
gas  fill  may  then  be  used  to  determine  the  blackbody  radiance  as 


i 


i 

i 

i 


j 


i 

i 

i 

I 

1 

! 


I 


61 


Lm  a 

T0g 


Y -cell 
Beam -ceil 


:](1  -P)5 


(48) 


where  is  the  Planck  function  evaluated  at  the  temperature  of  the  oven;  Y,  ''beam,''  and 
o 

"cell"  are  the  data  scans  mentioned  earlier;  and  p is  again  the  reflection  loss  at  a single 
sapphire -air  boundary.  Once  this  spectral  radiance  is  known,  the  blackbody  temperature  may 
be  calculated  from  the  inverse  Planck  function  as  in  Equation  (5),  and  compared  with  the  py- 
rometric  temperature  measurement. 


3.7  WAVELENGTH  CALIBRATION 

Wavelength  calibration  was  performed  in  the  2.7  pm  region  by  use  of  a small  burner 
placed  in  the  reference  beam,  which  burned  SFg  and  Hg*  HF,  a primary  combustion  product 
of  this  reaction,  has  an  absorption  band  in  this  region.  The  HF  line  positions  are  well  known 
and  the  lines  well  isolated  so  that  an  accurate  wavelength  calibration  was  easily  obtained. 

In  the  4.3  pm  region,  the  calibration  was  performed  by  observing  the  CO  fundamental  at 
2147  cm”*  (4.66  pm)  using  the  cold  cell  and  a path  length  of  100  meters.  The  CO  lines  are 
well  separated  and  their  positions  well  known  so  a calibration  could  be  easily  established. 

Constants  relating  the  optical  frequency  to  the  spectrometer  drum -number  were  used  in 
the  data  reduction  program  to  relate  the  digitized  experimental  data  to  the  optical  frequency. 
Data  were  collected  at  equal  drum -number  increments  so  that  knowledge  of  the  beginning  drum- 
number  of  each  scan  completely  defined  the  optical  frequency. 

3.8  DATA  REDUCTION 

In  addition  to  the  digitally  collected  data,  analog  chart  recordings  were  made  as  weU. 

These  were  used  as  preliminary  data  checks  to  see  that  the  system  was  operating  as  expected. 
They  were  particularly  useful  in  noting  drift  in  the  system.  A sample  analog  chart  is  shown 
in  Figure  9. 

The  data  scans  were  first  examined  for  drift.  Once  noted,  the  data  reduction  program 
automatically  applied  drift  corrections  based  on  the  noted  drift  punched  on  a header  card. 

The  initial  drum -number  of  the  scan  was  also  punched  on  a header  card;  subsequent  drum- 
number  changes  were  checked  by  coding  one  data  point  every  drum  revolution. 

The  appropriate  data  scans  — CELL,  BEAM,  ABSB,  ABSO,  Y (hot  through  cold),Y  (hot 
alone),  and  BEAM  (cold  cell  filled)  — were  read  into  the  University  of  Michigan's  360  system 
and  edited  where  necessary. 

The  data  reduction  program  then  matched  all  scans  for  drum -number  and  computed  the 
required  results  using  Equations  (44)  through  (47),  Results  were  obtained  for:  (a)  hot  cell 
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radiance,  (b)  cold  cell  transmittance,  (c)  hot  cell  radiance  as  viewed  through  the  cold  cell, 

(d)  hot  cell  absorptance,  and  (e)  hot  cell  apparent  temperature.  Both  data  input  and  data  out- 
put have  been  saved  on  tape.  The  outputs  were  automatically  plotted  on  the  plotting  system  at 
the  University  of  Michigan  Computing  Center. 

The  apparent  temperature  plot  was  checked  for  consistency  with  the  thermocouple  read- 
ings. When  apparent  temperatures  fell  within  about  25°C  of  the  thermocouple  temperature  in 
the  middle  of  the  band  (indicating  radiance  errors  of  less  then  7%)  and  remained  fairly  uniform 
across  the  band  (although  showing  deviations  at  both  band  edges),  the  data  were  assumed  to  be 
internally  consistent.  When  these  conditions  did  not  exist,  the  raw  data  were  re-examined  for 
unnoticed  drifts,  and  appropriate  changes  were  made.  If  the  data  could  not  be  made  to  satisfy 
these  consistency  criteria,  the  run  was  rejected  and  the  tests  repeated. 

Note  that  while  these  consistency  checks  verify  that  measured  emission  values,  absorption 
values,  and  temperature  values  all  make  sense,  they  have  no  bearing  at  all  on  the  accuracy  of 
the  cell  fill.  The  consistency  checks  monitor  the  radiometric  accuracy  of  the  measurements, 
but  the  manometer  system  is  relied  upon  for  the  accuracy  of  the  celi  filling  procedures. 

3.9  RESULTS 

During  the  program,  tests  were  made  with  many  different  cell  fills  and  the  results  of  15  tests 
are  included  in  this  report.  Table  3 describes  these  tests.  After  successful  operation  has  been 
confirmed  using  typical  water  vapor-nitrogen  fills  (Tests  8,  9,  and  11),  priorities  for  other  mix- 
tures were  assigned  by  mutual  agreement  between  ERIM  and  the  Air  Force,  and  subsequent  runs 
were  made  in  the  order  of  these  priorities.  There  were  insufficient  funds  to  obtain  good  data  on 
all  of  the  tests  planned.  Those  placed  at  the  top  of  the  priority  lists  simulated  real  amine-fuel 
rocket  exhausts  (Tests  15  through  20).  They  were  mixtures  of  water  vapor,  COg,  CO,  Hg,  and 
Ng  at  a total  pressure  of  0.1  atmosphere.  Measurements  were  made  at  both  2.7  /xm  and  4.3  /xm 
for  3 different  atmospheric  paths.  Subsequent  measurements  were  made  of  pure  water  vapor  at 
2.7  jam,  and  pure  CO^  at  2.7  /xm  and  4.3  /xm.  Detailed  test  conditions  are  shown  with  the  data  in 
Appendix  B.  The  lack  of  funds  eliminated  some  water  vapor  data  with  longer  hot  cell  optical 
paths  (Tests  12  through  14).  In  addition,  tests  to  check  the  assumption  that  the  correlation  be- 
tween lines  of  different  absorbing  gases  is  zero  (Tests  4 and  10)  were  planned  but  not  completed. 

The  results  obtained  from  a typical  water  vapor  run  at  2.7  /xm  are  shown  in  Figure  10(a-e). 
Five  plots  are  shown.  The  first  (a)  is  a plot  of  apparent  temperature  (T*)  versus  frequency.  It 
can  be  seen  that  the  apparent  temperature  was  within  about  ±15°K  of  the  hot  cell  temperature 
(1202°K)  over  all  of  the  band  where  significant  absorption  occurred;  and  at  the  ends  of  the 
band,  where  the  cell  began  to  become  transparent,  the  apparent  temperature  makes  large  ex- 
cursions to  either  side  of  the  cell  temperature.  This  is  the  proper  appearance  for  these  T* 
determinations.  The  temperature  deviation  of  ±15°K  represents  a radiance  uncertainty  of 
±5.6%,  which  is  a reasonable  accuracy  for  measurements  of  this  type. 
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TABLE  3.  SUMMARY  OF  MEASUREMENT  CONDITIONS 
(See  Appendix  B for  detailed  physical  parameters) 
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*Not  completed  because  of  time  and  funding  limitations. 
^Completed  before  priority  assignment. 

*See  tables  in  Appendix  B. 
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FIGURE  10.  TYPICAL  WATER  VAPOR  DATA  (Continued) 
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FIGURE  10.  TYPICAL  WATER  VAPOR  DATA  (Continued) 


The  second  (b)  and  third  (c)  curves  show  the  hot  cell  absorptance  and  the  cold  cell  trans- 
mittance. The  fourth  curve  (d)  shows  the  measured  hot  cell  radiance  (the  cold  cell  evacuated) 
on  the  upper  curve,  as  well  as  the  hot  cell  radiance  as  viewed  through  the  cold  cell  (lower 
curve).  This  plot  demonstrates  the  attenuation  produced  by  the  cold  cell.  The  fifth  plot  (e) 
demonstrates  the  degree  to  which  the  water  vapor  lines  in  the  hot  and  cold  cells  are  indeed 
correlated.  The  upper  curve  is  the  product,  wavenumber  by  wavenumber,  of  the  measured 
hot -cell  radiance — upper  curve  on  (d) — and  the  cold  cell  transmittance  (c).  The  lower  curve, 
which  is  the  measured  hot  cell  radiance  as  viewed  through  the  cold  cell,  is  the  same  as  the 
lower  curve  in  (d).  The  difference  between  these  two  curves  shows  the  degree  to  which  the 
multiplicative  transmittance  approach  fails.  In  the  center  of  the  band  the  correlation  effect 
produces  a radiance  which  is  lower  by  a factor  of  as  much  as  two  compared  to  what  one  would 
expect  using  the  multiplicative  property.  Detailed  discussion  of  these  effects  is  delayed  until 
Section  4 where  comparison  of  various  computation  techniques  is  made. 

Figure  ll(a-e)  show  typical  C02  data  at  4.3  4m  in  the  same  five  types  of  curves.  The 
apparent  temperature  curve  (a)  has  a behavior  similar  to  that  for  water  vapor  except  that  it 
does  not  show  the  fluctuations  at  the  edge  of  the  band.  This  is  expected  because  of  the  gener- 
ally smoother  character  of  the  spectra  in  this  region.  Figure  11(b)  through  (d)  shows  the  hot 
cell  absorptance  (b),  the  cold  cell  transmittance  (c),  as  well  as  the  hot  cell  radiance  (d,  upper) 
and  the  observed  hot  cell  through  cold  cell  radiance  (d,  lower).  Figure  11(e)  again  shows  the 
degree  to  which  line  correlation  effects  occur.  One  of  the  primary  features  of  this  and  all  the 
4.3  |im  C02  spectra  is  that  the  multiplicative  transmittance  approach  seems  to  work  well  for 
most  of  the  band  wings.  The  large  amount  of  absorption  in  the  simulated  paths  precludes 
judgments  concerning  band  center.  The  product  of  the  hot -cell  radiance  and  the  cold  cell  trans- 
mittance agrees  very  well  with  the  hot  cell  through  cold  cell  radiance  except  for  a very  small 
region  on  the  long  wavelength  side  of  the  "blue  spike"  at  2390  cm'1.  Thus,  it  appears  that 
correlation  effects  can  indeed  be  ignored  for  the  C02  4.3  am  band.  This  does  not  hold  for  the 
2.7  pm  COg  band,  however.  The  correlation  effects  observed  in  this  band  are  midway  between 
those  observed  for  water  vapor  and  those  observed  for  C02  at  4.3  pm.  Again,  detailed  discus- 
sion of  the  data  is  deferred  to  Section  4. 

Appendix  B contains  plots  of  all  of  the  data  taken  on  this  program  for  use  in  future  com- 
parisons with  calculation  techniques.  In  addition,  the  results  have  been  saved  on  computer  tape 
for  future  use. 


4 

COMPARISON  OF  MEASUREMENTS  WITH  PREDICTION  METHODS 

In  order  to  judge  the  ability  of  plume  calculation  methods  to  predict  the  line  correlation 
effects  that  are  present  (or  not  present)  in  the  measured  hot -through -cold  radiances,  it  is 
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WAVENUMBER  CCH-il 

(a)  Measured  Apparent  Temperature 
FIGURE  11.  TYPICAL  COg  DATA  (Continued) 


WAVENUMBER  CCM-1) 

(b)  Measured  Hot  Cell  Absorptance 
FIGURE  11.  TYPICAL  C02  DATA  (Continued) 


WAVENUMBER  CCH-ll 

(e)  Measured  Hot-Through-Cold  Radiance  (Lower)  and  Product  of  Hot  Cell  Radiance  and  Cold  Cell  Transmittance  (Upper) 


worthwhile  to  compare  the  measurements  made  to  the  band  model  development  of  Section  2 
and  to  other  Information  available  for  computing  these  effects. 

It  would,  of  course,  be  most  desirable  If  there  were  a single  computational  scheme  which, 
though  difficult,  could  be  relied  upon  to  give  absolutely  correct  results.  Although  the  line 
tabulation  of  McClatchey  might  appear  to  provide  the  basis  for  such  a technique,  It  has  been 
shown  to  lack  many  lines  whose  states  are  Important  at  high  temperatures^1^ . Therefore, 
for  the  time  being  the  tabulation  Is  considered  only  as  an  additional  basis  of  calculations. 

We  had  also  planned  to  compare  each  of  the  measurement  data  sets  to  Individual  line-by- 
llne  calculations  using  the  McClatchey  parameters.  These  calculations  represent  a substantial 
expenditure  of  funds,  and  when  It  appeared  that  they  could  not  be  considered  to  provide  defini- 
tive results,  we  sought  another  less  expensive  means  of  comparing  the  measurements  to  the 
line  tabulation. 

The  method  chosen  was  to  construct  a band  model  from  the  line  tabulation.  Goldman  and 
Kyle^1^  state  that  such  a model  fits  the  llne-by-line  calculation  to  within  5%.  This  is  within 
the  accuracy  to  be  expected  from  the  measurements  themselves.  The  band  model  parameters 
were  computed  for  the  random  model  with  an  inverse -strength  distribution  of  line  strengths, 
as  is  used  in  Section  2.  These  band  model  parameters  were  then  checked  with  line-by-line 
calculations  over  a single  small  spectral  interval  for  one  of  the  water  vapor  cases. 

The  other  set  of  emission -absorption  information  with  which  the  data  are  compared  con- 
sists of  model  calculations  based  on  an  extensive  tabulation  of  band  model  parameters  devel- 
oped by  Ludwig,  et  al.  and  most  recently  tabulated  in  Reference  1. 

In  addition  to  the  emission-absorption  parameters,  there  is  an  additional  variable  which 
affected  the  agreement  of  calculation  and  experiment.  This  is  the  curve -of -growth  used  with 
the  various  calculational  techniques.  Since  the  curve -of -growth  to  be  used  is  a direct  outcome 
of  the  band  model  to  be  used  (e.g.,  Equation  (39)  of  Section  2),  one  would  not  expect  it  to  be  an 
independent  variable.  We  consider  it  as  such  since  it  affects  the  fit  between  the  calculated 
and  measured  hot  cell  radiance  substantially. 

There  are  three  ways  in  which  the  two  sets  oi  parameters  may  be  used  in  computational 
schemes  for  non-isothermal  paths.  The  first  means  is,  of  course,  the  band  model  scheme 
developed  in  Section  2.  The  second  is  the  Lindquist -Simmons  approximation  developed  in  the 
pastel . The  third  is  the  multiplicative -transmittance  approach  described  in  the  early  portion 
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14.  G.H.  Lindquist,  High  Resolution  Measurements  of  Water  Vapor  Spectrum,  ongoing  con- 
tract on  ARPA  Order  2656,  Environmental  Research  Institute  of  Michigan,  Ann  Arbor, 
Michigan,  April  1974 -August  1975. 

15.  A.  Goldman  and  T.G.  Kyle,  "A  Comparison  Between  a Statistical  Model  and  Line-by-Line 
Calculations  with  Application  to  the  9.6  pm  Ozone  and  the  2.7  pm  Water  Vapor  Bands,'1' 
Applied  Optics,  Vol.  7,  p.  1167f  (1968). 
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of  Section  2 in  which  the  radiance  of  the  hot  gas  source  is  simply  multiplied  by  the  cool  path 
transmittance  over  small  wavenumber  intervals.  The  measurements  show  that  this  latter 
technique  is  not  adequate  for  water  vapor  if  a factor -of -2  error  is  unacceptable.  On  the  other 
hand,  the  measurements  show  the  technique  to  be  satisfactory  for  the  CO2  wings  at  4.3  pm. 

An  additional  comparison  which  should  be  made  is  one  in  which  the  statistical  treatment 
of  Section  2 is  compared  with  a realistic  ensemble  of  lines.  In  this  case,  of  course,  line-by- 
line calculations  are  necessary.  Here  the  line  tabulation  serves  as  a realistic  array  of  vary- 
ing strength,  width,  position,  and  temperature  dependence,  all  falling  in  a narrow  frequency 
interval.  The  statistical  treatment  of  Section  2,  used  in  obtaining  a transmittance  derivative 
for  this  array  of  lines,  can  be  compared  to  the  true  transmittance  derivative  calculated  using 
a line-by-line  procedure. 

The  following  sections  describe  the  comparisons  which  we  performed  using  these  varia- 
bles. Because  of  the  large  number  of  variables,  we  did  not  obtain  all  possible  comparisons 
but  tried  to  obtain  enough  of  them  to  ascertain  the  degree  to  which  band  model  calculations 
could  be  used  to  reproduce  the  measurements. 

4.1  THE  MEASUREMENTS  COMPARED  TO  THE  GENERAL  DYNAMICS  PARAMETERS 

The  parameters  are  most  easily  checked  with  respect  to  emission  and  absorption  data  by 
comparing  measured  and  computed  hot  cell  radiance  values  and  cold  cell  transmittance 
values.  Each  path  (hot  and  cold)  alone  represents  an  isothermal  path  whose  emission  or  ab- 
sorption depends  only  on  the  parameters  used. 

The  band  model  tabulations  of  General  Dynamics^  are  designed  for  use  in  band  model 
calculations  of  hot  gases.  For  COg  these  tabulations  are  the  result  of  theoretical  calculations 
which  used  the  known  spectroscopic  constants  for  C02  (to  give  accurate  band  positions)  but 
used  harmonic  oscillator  approximations  to  obtain  hot  band  strengths  from  the  experimentally 
determined  strengths  of  the  ground  state  bands.  These  tabulations  are  generally  assumed  to 
be  accurate  for  COg. 

For  water  vapor,  parameter  tabulations  were  determined  experimentally.  A substantial 
number  of  measurements  of  the  isothermal  absorption  of  hot  water  vapor  were  made  using  a 
variable  length  strip-burner  as  the  source  of  water  vapor.  These  provided  absorption  data 
above  1200°K.  Data  below  this  temperature  were  based  both  on  extrapolation  and  comparison 
with  other  low  temperature  measurements.  Values  as  published  did  not  fit  th  .xtension  am- 
bient temperature  band  model  tabulation  of  Anding,  but  under  a later  effort  the  300°K  values 
were  modified  to  agree  with  the  model  of  Anding.  The  modification  was  not  entirely  success- 
ful, as  will  be  seen  by  comparison  with  the  cold  cell  measurements. 

The  General  Dynamics  tabulation  contains  values  of  the  average  strength-to-spacing 
ratio  and,  as  the  second  band  model  parameter,  the  average  line  spacing  which  could  be 
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combined  with  the  average  line  width  to  give  a second  parameter.  The  line  strength -to- 
spacing  ratio  is  the  first  band  model  parameter  in  most  band  model  formulations,  including 
the  one  developed  in  Section  2 and  the  one  which  we  have  commonly  used  prior  to  this  work. 
Under  previous  efforts,  the  average  line  spacing  information  had  been  combined  with  average 
half  widths  to  obtain  the  dimensionless  parameter  0: 


where  the  parameter  1/5  was  obtained  from  the  GD  tabulations.  The  parameter  y used  in 
Section  2 is  related  by  0 by 


General  Dynamics  extracted  a value  oi  1/ 5 from  calculations  and  absorption  measurements 
using  a curve  of  growth  of  a form  equivalent  to 


W/  6 = (3 


V1 +lx 


where  x is  defined  as  Jkda//3.  Values  of  k and  f?  were  used  directly  in  the  comparison  of  the 
data  with  the  General  Dynamics  tabulations.  Some  dependence  on  the  curve -of -growth  was 
apparent,  as  will  be  discussed  in  Section  4.3. 


During  the  measurement  program,  only  three  different  hot  cell  fills  were  used:  (1)  a 0.1 
atmosphere,  50-50  mixture  of  water  vapor  and  nitrogen,  (2)  a 0.1  atmosphere,  50-50  mixture 
of  COg  and  nitrogen,  and  (3)  a 0.1  atmosphere  mixture  simulating  a rocket  exhaust.  Compar- 
ison of  these  hot  cell  radiances  with  computed  values  using  the  General  Dynamics  parameters 
are  shown  for  typical  tests  in  the  upper  curves  of  Figures  12  through  16.  Comparisons  for 
additional  tests  would  only  yield  information  on  measurement  reproducibility,  which  was 
found  to  be  quite  good. 


Figure  12  compares  radiances  resulting  from  the  General  Dynamics  parameters  to  the 
measurements  made  on  pure  water  vapor  at  2.7  /i m.  In  general,  prediction  by  these  param- 
eters is  very  good  at  1200°K.  The  predicted  radiance  is  slightly  low,  but  the  choice  of  a 
curve  of  growth  affects  these  values  greatly  because  the  optical  depths  near  band  center  fall 
intermediate  between  the  linear  and  square  root  region,  where  the  different  curves  of  growth 
have  their  greatest  variability.  (See  Section  4.3  for  an  explanation  of  the  curves  of  growth 
used.) 


Figure  13  shows  the  hot  cell  radiance  measured  for  a COj-Ng  mixture  at  4.3  (im  com- 
pared to  the  predicted  results  using  the  General  Dynamics  parameters  for  Test  5.  Although 
the  shape  of  the  observed  spectrum  agrees  with  the  predicted  shape,  the  observed  radiance  at 
the  peak  of  the  band  is  substantially  less  than  the  predicted  radiance.  The  observed  and 
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FIGURE  12.  COMPARISON  OF  MEASURED  HOT  CELL  RADIANCE  WITH  RADIANCE  CALCULATED  USING  GENERAL 

DYNAMICS  PARAMETERS,  TEST  8R,  2.7  pim  WATER  VAPOR 
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FIGURE  13.  COMPARISON  OF  MEASURED  HOT  CELL  RADIANCE  WITH  RADIANCE  CALCULATED  USING  GENERAL 

DYNAMIC 3 PARAMETERS,  TEST  5,  4.3  nm,  C02 


predicted  curves  cross  in  the  long  wavelength  vnrg  so  that  the  far  wing  is  underpredicted 
rather  than  overpredicted.  This  latter  feature  produces  severe  disagreement  between  the 
observed  and  predicted  hot -through -cold  radiance  using  these  parameters.  It  is  apparent 
that  the  General  Dynamics  parameters  for  the  COg  region,  although  producing  gross  agree- 
ment with  experiment,  give  substantial  variances  in  the  detailed  shape  of  the  band,  so  much 
that  the  red  spike  region  is  distorted.  This  is  surely  the  result  of  the  very  heavy  dependence 
on  theory  in  obtaining  these  parameters. 

Figure  14  shows  the  hot  cell  radiance  measured  for  a COg-Ng  mixture  at  2.7  pm  com- 
pared to  predicted  results  using  the  GD  tabulation  for  Test  3.  It  can  be  seen  that  the  General 
Dynamics  tabulation  predicts  the  observed  hot  cell  radiances  very  poorly.  Reasons  for  this 
are  not  known  except  that  perhaps  they  are  again  the  result  of  the  heavy  dependence  upon 
theory  in  obtaining  the  parameters.  The  shape  of  the  spectra  are  well  predicted,  but  the  ab- 
solute values  are  not.  We  believe  that  the  observed  spectra  accurately  characterize  the  true 
gas  radiance. 

Figures  15  and  18  show  comparisons  for  the  hot  cell  fill  simulating  a rocket  plume.  At 
4.3  |im  (see  Figure  15)  the  prediction  is  within  25%  on  the  long  wavelength  side  of  the  band, 
but  the  agreement  at  band  center  is  much  better  than  was  observed  for  the  COg-Ng  runs. 

This  would  appear  to  indicate  that  the  treatment  of  broadening  is  not  correct  in  our  calcula- 
tion procedure.  The  broadening  coefficients  are  poorly  known  for  and  essentially  unknown 
for  water  vapor  and  hydrogen,  and  it  is  difficult  to  identify  the  reason  for  the  band  center  dis- 
agreement between  Figures  13  and  15  without  making  measurements  on  pure  COg.  Such  mea- 
surements are  necessary  to  separate  errors  in  the  treatment  of  foreign  gas  broadening  from 

■*r 

errors  in  the  band  model  parameters  themselves.  It  appears  that  a study  of  foreign  gas 
broadening  of  COg  is  needed,  as  well  as  a better  determination  set  of  band  model  parameters, 
if  the  agreement  in  the  4.3  g m band  is  to  be  improved. 

The  comparison  for  the  simulated  rocket  plume  mixture  at  2.7  pm  is  shown  in  Figure  16. 
The  water  vapor  features  of  the  band  are  reasonably  well  treated,  commensurate  with  the  re- 
sults for  water  vapor-Ng.  The  COg  features  are  also  reasonably  well  predicted.  That  we  can 
make  this  statement  is  very  surprising  in  view  of  the  very  poor  agreement  observed  in  the 
COg-Ng  runs.  Again  this  indicates  that  the  broadening  of  COg  is  poorly  treated  using  the 
available  broadening  coefficients. 

Three  basic  atmospheric  paths  were  used  in  the  program  simulating  10,  50,  and  200  km 
paths  at  15  km  altitude.  Figures  17  through  22  show  the  measured  cold  cell  transmittances 
compared  to  GD  predictions.  From  Figures  17  and  22,  it  can  be  seen  that  predictions  for 
water  vapor  are  reasonably  good,  perhaps  underestimating  the  absorption  slightly.  On  the 
•)th*>r  hand,  the  results  for  COg  are  not  very  good  at  all.  The  band  appears  slightly  shifted  in 
some  cases,  and  at  intennediate  optical  depths  substantial  disagreement  is  apparent.  Because 


FIGURE  17.  COMPARISON  OF  MEASURED  COLD  CELL  TRANSMITTANCE  WITH  TRANSMITTANCE  CALCULATED 
USING  GENERAL  DYNAMICS  PARAMETERS,  TEST  20,  2.7  nm,  10  km  PATH  AT  15  km  ALTITUDE 
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FIGURE  21.  COMPARISON  OF  MEASURED  COLD  CELL  TRANSMITTANCE  WITH  TRANSMITTANCE  CALCULATED 
USING  GENERAL  DYNAMICS  PARAMETERS,  TEST  16,  4.3  lira,  50  km  PATH  AT  15  km  ALTITUDE 
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the  GD  parameters  for  water  vapor  at  these  low  temperatures  are  to  some  extent  extrapola- 
tions while  for  COg  they  are  based  primarily  on  theoretical  calculations,  poor  agreement  at  low 
temperatures  may  not  be  surprising. 


4.2  THE  MEASUREMENTS  COMPARED  TO  THE  AFCRL  LINE  PARAMETER  TABULATION 
IN  A BAND  MODEL 

In  order  to  facilitate  the  comparison  between  the  line  tabulation  of  AFCRL  and  the  mea- 
surements of  Section  3,  a band  model  was  developed  from  the  line  tabulation  using  the  tech- 
niques described  in  [4] . The  band  model  was  developed  for  use  with  a curve  of  growth  like 


W/d  = fl.  (49) 

where  again  x = j kda/0.  (See  footnote  p.  39.)  The  parameters  k and  fi  were  determined  as  a 
function  of  temperature  and  for  one  atmosphere  pressure  of  pure  absorber  using  the  following 
relations: 

N 

k(T>  * V 


fy i(T> 

■ ^V-  150 

where  is  the  spectral  resolution  used  for  the  band  model.  A value  of  25  cm  1 was  used 
fjr  A?  in  the  2.7  4m  region  and  10  cm  * was  used  in  the  4.3  4m  region.  In  order  to  deter- 
mine the  strengths  as  a function  of  temperature,  a harmonic  oscillator  vibrational  partition 
function  was  used  in  addition  to  the  rotational  partition  function  given  in  Reference  2.  This 
proved  to  be  quite  important  in  obtaining  reasonable  results,  especially  for  CC>2  which  has  a 
rather  low-frequency  bending  vibration. 

The  second  parameter  (i  was  determined  using  the  relation 


where  and  signify  values  at  one  atmosphere  pressure.  yQ.(T),  which  is  the  halfwidth 
for  pure-gas  broadening  has  an  inverse  square  root  temperature  dependence.  Values  were 
obtained  from  the  air -broadened  halfwidths  contained  in  the  AFCRL  tabulation  by  division  by 
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the  appropriate  broadening  coefficient — 0.2  in  the  case  of  water  vapor  and  0.77  for  CC^.  The 
coefficient  2n  in  Equation  (51)  differs  from  the  coefficient  of  1/4  given  in  Reference  3,  Equa- 
tion (4.41),  in  order  that  the  proper  square  root  limit  is  approached  using  the  curve  of  growth 
in  Equation  (49). 

Comparisons  of  the  hot  cell  radiance  using  these  hand  model  parameters  are  shown  in 
Figures  23  through  25.  In  the  2.7  pm  water  vapor  region,  (see  F'gure  23)  the  agreement  with 
hot  cell  radiance  measurements  is  reasonably  good  near  band  center,  except  near  3700  cm’1 
where  the  line  parameters  predict  substantially  more  radiance  than  observed.  In  the  wings 
of  the  band,  however,  the  agreement  is  not  so  good.  In  the  long  wavelength  wing  between  3000 
and  3300  cm  1 agreement  is  poor-  Presumably  this  is  because  the  line  tabu.ation  has  not 
been  carried  out  to  high  enough  rotational  quantum  numbers  for  the  1200°K  cell  temperature. 
The  agreement  at  band  center  is  particularly  satisfying,  however,  and  is  better  at  some  fre- 
quencies than  that  obtained  using  the  GD  parameters. 

The  band  model  parameters  obtained  from  the  line  tabulation  yield  good  results  at  the 
short  wavelength  side  of  the  4.3  pm  CC>2  band  (see  Figure  24)  but  do  very  poorly  at  predicting 
the  middle  and  long-wavelength  portions  of  the  bard.  This  again  is  probably  due  to  the  omis- 
sion from  the  tabulation  of  many  lines  which  are  important  at  U mperatures  near  1200°K. 

Throughout  the  2.7  pm  band  of  C02  (see  Figure  25)  agreement  is  poor.  Radiance  is  un- 
derpredicted on  the  basis  of  the  line  tabulation  to  about  the  sanu  degret  as  :t  is  overpredicted 
based  on  the  GD  parameters.  The  large  amount  of  disagreement  in  this  band  is  surprising 
and  unexplainable  at  present. 

The  agreement  of  the  band  model  parameters  from  the  CR1.  t.itmlation  wit1,  the  cool  cell 
transmittance  values  is  basically  good,  as  one  would  expect.  Fi  -um  .-  2'  through  31  demon- 
strate this.  Several  features  of  the  comparison  are  puzzling  however.  First,  m the  2.7  pm 
COg  band  the  cold  ceil  transmittance,  like  the  hot  c<  il  radi  u..  e.  is  not  acci  afh  predicted. 
Again,  measurements  on  pure  C02  paths  are  needed  to  sort  1 d • s*'  ut.i  < it.u  ties.  Second, 
as  seen  in  Figure  27,  a portion  of  the  water  vajxjr  band  is  not  u>  !!  pr<  dieted,  also  for  unknown 
reasons.  Third,  the  measured  short  wavelength  shoulder  o!  ti  ■ i 3 .in  TO.,  band  shows  suh- 

_t 

Stantial  absorption  below  the  shortest  wavelength  band  head  at  23.'.  m ‘ Flu  calculations 
show  an  abrupt  enci  to  the  absorption  at  this  wavelength.  Vhe  observe.!  .ih*.  opt  mi.  could  in- 
produced  by  the  integrated  effect  of  the  line  wings  of  tin  stiong  lines  in  t!.t  nano.  Band  model 
methods  do  not  easily  treat  these  effects,  except  in  as:  empirical  fashion:  neitiu  r will  line -by- 
line calculations  treat  them  unless  the  far  line  shapes  are  .tccui.itt  1\  kin-wn.  Thus  m accept 
them  as  real,  unaccounted-for  effects  for  which  an  • mpiruai  r »rrectioi.  m st  1 <■  auded  over 
long  atmospheric  paths. 
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FIGURE  23.  COMPARISON  OF  MEASURED  HOT  CELL  RADIANCE  WITH  RADIANCE  CALCULATED  USING  PARAMETERS 
OBTAINED  FROM  AFCRL  TABULATION,  TEST  8,  2.7  Mm,  WATER  VAPOR 
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FIGURE  94.  COMPARISON  OF  MEASURED  HOT  CELL  RADIANCE  WITH  RADIANCE  CALCULATED  USING 

OBTAINED  FROM  AFCRL  TABULATION  TEST  7,  4.S  am,  C02 
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4.3  EFFECT  OF  THE  CURVE  OF  GROWTH  ON  THE  AGREEMENT  BETWEEN  THE 

MEASUREMENTS  AND  THE  CALCULATIONS 

We  did  not  expect  the  calculations  to  be  greatly  dependent  on  the  curve  of  growth.  How- 
ever, the  hot  cell  radiances  fell  in  regions  of  the  curve  of  growth  intermediate  between  the 
linear  and  square  root  regions  at  many  of  the  frequencies.  Where  this  was  the  case  the  com- 
puted results  were  affected  by  the  curve  of  growth  used.  For  the  General  Dynamics  compar- 
ison shown  in  Figures  12  through  22  the  curve  of  growth  of  Equation  (49)  was  used.  This 
curve  of  growth,  which  is  valid  for  a random  model  with  an  exponential  distribution  of  lines, 
is  the  form  used  to  obtain  the  General  Dynamics  parameters  from  experimental  data  in  the 
first  place.  Other  curves  of  growth  that  could  have  been  used  are  (a)  the  function  valid  for  an 
inverse -strength  distribution,  Equation  (39),  and  (b)  the  Ladenburg-Reiche  function,  valid  for 
band  of  lines  of  equal  strengths.  The  curve  of  growth  of  Equation  (49)  produces  the  best  fit 
between  the  General  Dynamics  parameters  and  the  measured  data.  On  the  one  hand,  this  is 
not  surprising  because  of  its  use  by  GD  in  fitting  the  data.  On  the  other  hand,  if  General 
Dynamics  were  indeed  able  to  accurately  determine  the  square  root  limit  for  each  frequency, 
then  the  GD  parameter  determination  should  have  been  independent  of  the  curve  of  growth,  in 
which  case  the  best -fit  curve  of  growth  for  these  measurements  should  be  that  which  describes 
the  line  distribution  best.  Figure  3 showed  that,  for  water  at  3550  cm'1  at  least,  the  line 
strength  distribution  is  well  represented  by  an  inverse -intensity  function.  As  a result,  one 
would  expect  that  Equation  (39)  would  provide  the  best  curve  of  growth.  Since  it  does  not,  it 
can  only  be  surmised  that  GD  was  not  able  to  determine  the  actual  square  root  limit  for  water 
vapor  but  only  approximated  it  by  extrapolation  using  the  curve  of  growth  of  Equation  (49). 

The  parameters  obtained  from  the  AFCRL  tabulation  do  include  the  proper  square  root 
limit  by  definition,  provided  that  the  line  strengths  and  widths  are  complete  and  accurate. 

(They  appear  to  be  accurate  near  band  center.)  Hence  we  would  expect  the  best  curve -of  - 
growth  to  be  that  which  corresponds  most  closely  to  the  true  statistical  character  of  the  band, 
again  that  of  Equation  (39).  This  curve -of -growth  is  indeed  best  for  some  portions  of  the 
band,  but  not  for  others.  Figure  32  presents  predictions,  using  parameters  obtained  from  the 
AFCRL  tabulations  and  based  on  the  curves  of  growth  of  both  Equations  (39)  and  (49),  com- 
pared to  measured  data. 

The  discussion  here  only  serves  to  point  out  the  problems  involved  in  the  choice  of  a 
curve-of-growth.  It  is  not  possible  to  separate  the  effects  of  possible  failure  of  the  assump- 
tions of  the  band  model  (e.g.,  non-random  spacing,  coincidental  line  occurrences,  etc.)  with 
the  possible  effects  of  missing  lines  in  the  case  of  the  AFCRL  tabulation  or  with  effects  of 
errors  in  parameter  determination  in  the  case  of  the  General  Dynamics  parameters. 

Thus,  the  conclusion  to  be  made  from  these  sections,  at  least  for  isothermal  paths,  is 
that  set  of  parameters  and  curve-of-growth  should  be  used  which  give  the  best  agreement 
with  measurement. 
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4.4  HOT -THROUGH -COLD  RADIANCE  AND  THE  EFFECT  OF  THE  CALCULATIONAL 

METHOD 

In  Section  2 a band  model  was  developed  specifically  to  handle  the  very  non -isothermal 
conditions  expected  in  computing  the  hot  gas  radiance  as  viewed  through  cold  gases.  A pre- 
vious model^  was  also  developed  to  treat  such  cases,  but  it  was  aver  tested  on  real  situa- 
tions. In  this  section  the  two  calculation  techniques  will  be  compared  to  experimental  data 
* for  hot -through-cold  radiances.  We  had  hoped  to  make  this  comparison  independent  of  the 

choice  of  spectroscopic  parameters  (GD  parameters  or  AFCRL  line  tabulation),  but  there  was 
enough  variability  in  the  degree  to  which  these  parameters  fit  the  hot  cell  radiance  measure- 
ments to  produce  variations  in  the  way  the  two  calculational  methods  predicted  the  radiances. 

Figure  33  shows  comparisons  of  the  Lindquist -Simmons  approximations^  of  bottom 
data  points,  the  band  model  of  Section  2 (middle  data  points),  and  the  multiplicative  - 
transmittance  approach  (upper  data  points),  along  with  measured  hot -through -cold  radiances 
(curve)  for  a typical  water  vapor  run.  The  General  Dynamics  parameters  and  the  curve -of - 
growth  of  Equation  (49)  were  used  for  all  of  these  comparisons  (even  though  the  curve -of - 
growth  of  Equation  (39)  applies  strictly  to  the  band  model  of  Section  2).  The  case  chosen  was 
that  representing  the  presence  of  the  largest  line  correlation  effects.  It  can  be  seen  that  the 
band  model  of  Section  2 agrees  best  with  measurements,  although  the  values  it  produces  are 

!■  systematically  low. 

Figure  34  presents  the  same  data,  but  based  on  the  band  model  parameters  obtained  from 
* the  AFCRL  line  tabulations.  Here,  near  band  center,  the  band  model  of  Section  2 represents 

the  measured  data  very  closely.  Near  the  band  edge  none  of  the  computed  curves  fits  the 
data  well,  but  this  is  again  representative  of  the  line  tabulations  failure  to  predict  band  edges 
because  of  the  lack  of  upper  state  lines. 

Figure  35  shows  a curve  similar  to  Figure  33  for  the  4.3  pm  COg  band.  Here,  the  cor- 
relation effect  is  known  to  be  negligible  as  discussed  in  Section  3.  The  three  predictions  are 
similar  except  for  minor  variations  as  might  be  expected  based  on  the  known  lack  of  correla- 
tion effects.  The  results  do  not  agree  well  with  the  data  however  because  of  the  lack  of  agree- 
ment with  the  GD  band  model  parameters  as  outlined  in  Section  4.1. 

Figure  36  shows  curves  similar  to  Figure  34  for  the  4.3  pm  COg  band.  Here  the  AFCRL 
parameters  are  used  for  the  predictions.  This  curve  demonstrates  how  poorly  the  line  tabu- 
l lation  predicts  the  red  spike.  The  three  calculations  show  some  effect  of  line  correlation  (the 

> upper  set  of  points  is  well  separated  from  the  lower  set  of  points).  This  indicates  that  the 

lines  in  the  tabulation  are  not  so  densely  spaced  as  to  eliminate  correlation  effects  at  1200°K, 
whereas  the  measurements  indicate  no  correlation  effect  in  these  regions  of  the  band.  This 
is  further  evidence  for  missing  lines.  The  disagreement  between  measurement  and  each  of 
the  predictions  is  large  enough  that  one  cannot  judge  the  best  calculational  technique. 
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Multiplicative  Transmittance 
(upper)  Band  Model  of  Section  2 
(lower)  Lindquist  - Simmons  Approx. 


C AJTTrr 


i f— ) 


01  X (t  W3-HS- W3/M) 

1“  c 


3QNUIQUU  TUUi03dS 


FIGURE  35.  COMPARISON  OF  THE  HOT -THROUGH- COLD  RADIANCE  AS  MEASURED  AND  PREDICTED  BY 
3 METHODS  FOR  THE  4.3  fin  C02  BAND,  USING  THE  GENERAL  DYNAMICS  PARAMETERS  (TEST  7) 


Figures  37  and  38  show  curves  for  the  2.7  pm  COg  band.  Since  correlation  effects  are 
important  here,  the  three  predictions  again  give  different  results.  Here  again,  the  disagree- 
ment between  the  measurements  and  the  available  parameters  is  so  great  that  a judgment 
concerning  the  calculational  techniques  cannot  be  made. 

From  the  foregoing  figures  it  can  be  seen  that  the  water  vapor  region  at  2.7  /im  is  the 
only  region  where  both  (1)  the  band  model  parameters  are  consistent  with  the  data,  and  (2) 
where  correlation  effects  are  present.  It  is  only  under  these  conditions  that  the  three  calcu- 
lational techniques  (the  band  model  of  Section  2,  the  Lindquist -Simmons  approximation,  and 
the  multiplicative  transmittance  approach)  can  be  fairly  compared.  In  this  region  the  band 
model  of  Section  2 does  indeed  do  a better  job  of  prediction  than  the  other  two  methods.  How- 
ever, the  agreement  between  this  new  method  and  the  experimental  data  is  not  as  good  as  it 
might  be.  We  currently  attribute  this  uncertainty  to  the  uncertainties  in  the  emission - 
absorption  parameters  compared  to  the  measured  hot  gas  radiance  values. 


4.5  COMPARISON  OF  BAND  MODEL  TO  AN  ARRAY  OF  WATER  VAPOR  LINES 

A much  more  satisfying  way  of  evaluating  the  band  model  of  Section  2 would  be  to  com- 
pare it  to  a specific  array  of  absorption  lines.  Such  an  array  need  not  correspond  exactly  to 
a real  band  in  order  to  provide  a test  of  the  model.  We  have  performed  such  an  evaluation 
for  a spectral  resolution  element  25  cm”1  wide  centered  at  3550  cm”1. 

It  was  done  by  computing  the  transmittance  derivative  with  infinitely  fine  resolution  using 
the  AFCRL  line  tabulation,  smoothing  the  result  to  25  cm”1,  and  comparing  the  result  to  the 
transmittance  derivative  of  Equation  (34),  which  was  evaluated  using  the  band  model  param- 
eters determined  from  the  line  tabulation. 

The  line-by-line  calculations  of  the  transmittance  derivative  were  made  using  the  same 
homogeneous  path  approximations  as  were  made  in  the  development  of  Section  L This  pro- 
cedure made  the  comparisons  free  of  any  effects  introduced  by  non -homogeneous  paths  and 
limited  the  comparison  to  the  statistical  properties  of  the  arrangement  of  the  array  of  lines. 
The  following  equation  was  used  in  the  line-by-line  calculation  of  the  transmittance  derivative. 


N S,(T)y,(T,p) 


’ ' fcfr? <t,p>  * [*  - *0 J2  [rfr.P>  * (I  - \) 


S-ODy^T.^a 


(52) 


where  T,p  are  the  temperature  and  broadening  pressure  of  the  isothermal  cool  path,  and  a is 
its  absorbing  length.  T,p  are  the  local  temperature  and  broadening  pressure  of  the  hot  emit- 
ting path. 

Equation  (52)  was  used  to  compute  the  transmittance  derivative  for  several  cases  repre- 
senting different  optical  path  lengths  corresponding  to  several  values  of  the  quantity  u,  as 
well  as  for  a single  combination  of  T = 1200,  T = 300,  P = 0.0379  atmosphere  of  water  vapor, 
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and  P = 0.003  atmosphere  of  water  vapor.  No  foreign  gas  broadening  was  used  in  the  calcula- 
tion. This  gave  the  conditions  q/q  = 1/4,  y/y  = 8 so  that  y/y  = 2 in  the  band  model  of  Sec- 
tion 2.  The  transmittance  derivative  computed  by  Equation  (52)  was  smoothed  with  a square 
slit  function  25  cm”*  wide  to  correspond  to  the  averaging  width  used  to  compute  band  model 
parameters  from  the  AFCRL  tabulation. 

Figure  39  shows  the  average  equivalent  width  derivative  as  predicted  by  the  Lindquist - 
Simmons  approximation,  the  band  model  of  Section  2,  and  these  line-by-line  calculations. 

The  results  from  the  line-by-line  calculations  fall  between  the  two  band-model  calculation 
methods,  and  they  are  a result  of  the  fact  that  neither  set  of  assumptions  is  completely  satis- 
fied by  the  group  of  lines  contributing  to  the  line-by-line  calculations.  Calculations  at  addi- 
tional frequencies  are  sure  to  give  different  results,  and  more  extensive  comparisons  are 
needed  before  judgments  are  made.  (Because  the  development  of  Section  2 was  not  completed 
until  the  very  end  of  the  effort,  comparisons  could  not  be  exhaustive.)  Any  substantial  dis- 
agreement between  the  band  model  of  Section  2 and  the  line-by-line  calculations  which  persist 
is  most  probably  the  result  of  the  assumed  relationship  between  E and  SQ.  This  relationship 
could  be  adjusted  to  provide  optimum  agreement  with  line-by-line  calculations  once  a defini- 
tive comparison  is  made.  Such  a task  should  receive  high  priority  in  future  work. 


5 

CONCLUSIONS  AND  RECOMMENDATIONS 

Because  there  was  a mixture  of  experimental  and  theoretical  efforts  in  this  work,  there 
are  a great  many  conclusions. 

The  evidence  for  these  conclusions  is  not  complete  in  some  cases,  but  several  significant 
statements  can  be  made  concerning  plume  calculation  techniques.  We  have  organized  these 
conclusions  into  four  categories. 

5.1  THE  PRESENCE  OF  CORRELATION  EFFECTS 

It  is  concluded  from  the  experimental  results  that  in  the  2.7  pm  water  bands,  the  accu- 
rate calculation  of  plume  radiation  at  a distance  requires  that  line  correlation  effects  be 
accounted  for  if  errors  oi  the  order  of  a factor  cf  2 are  to  be  avoided.  Based  on  these  mea- 
surements, we  expect  that  line  correlation  effects  will  be  important  for  this  band  independent 
of  the  optical  path  involved.  Therefore,  the  calculation  techniques  used  should  be  able  to 
treat  these  effects.  Although  its  accuracy  has  not  been  completely  verified  during  this  pro- 
gram, the  band  model  of  Section  2 does  appear  to  do  a creditable  job  of  handling  correlation 
effects. 
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The  C09  band  at  4.3  pm  is  not  so  easily  dismissed.  At  the  moderate  altitude  conditions 
used  in  this  program  (about  0.1  atmosphere)  it  should  be  valid  to  ignore  line  correlation 
effects  in  the  long  wavelength  wing.  The  measurements  confirm  this  and  the  band  model 
parameters  predict  it.  Thus  at  0.1  atmospheres  or  more,  even  the  multiplicative -transmittance 
approach  provides  realistic  results  for  plume  radiation  through  an  atmospheric  path.  The 
reason  for  this  cannot  be  that  the  emission  and  absorption  lines  are  uncorrelated.  Rather,  it 
must  be  that  the  lines  are  wide  enough  and/or  so  densely  spaced  that  either  the  emission  or 
absorption  spectrum  (or  both)  takes  on  a continuous  appearance.  That  this  is  almost  certainly 
the  case  can  be  assessed  from  the  band  model  parameter  ft  (or  equivalently  y),  which  is  a 
measure  of  approach  to  square -root  growth  where  line  overlap  is  important.  The  General 
Dynamics  tabulations  for  the  4.3  pm  C02  band  indicate  values  for  ft  of  the  order  of  30  to  50  at 
1200°K  and  100  or  larger  at  higher  temperatures.  The  fact  that  these  values  are  so  large  im- 
plies line  overlapping.  At  temperatures  below  1200°K  the  values  are  not  large  and  are  gener- 
ally less  than  10  at  sea  level  pressure.  Thus,  as  temperature  and  pressure  are  decreased, 
the  GD  parameters  predict  that  correlation  effects  will  become  more  important.  The  current 
measurements,  limited  in  temperature  and  pressure  by  the  sensitivity  of  the  system  as  it  was 
operated  during  this  program,  were  not  able  to  include  hot -through -cold  radiance  measure- 
ments at  lower  temperatures  than  1200°K,  nor  at  optical  depths  smaller  than  that  produced 
by  60  cm  of  C02  at  a 0.05  atmosphere  pressure.  Thus,  the  possible  correlation  effects  in  the 
4.3  pm  band  occur  at  lower  temperatures  (and  pressures)  than  could  be  produced  in  these  ex- 
periments and  remain  essentially  unexplored.  Measurements  specifically  aimed  at  the  con- 
ditions where  correlation  effects  might  be  important  could  be  made  witli  the  instrumentation 
used  in  this  program,  with  addition  of  a grating  which  gives  optimum  dispersion  around  4.3 
pm  and  by  sacrificing  some  of  the  spectral  resolution.  (The  4 cm  spectral  resolution  used 
in  this  region  is  finer  than  needed.) 

Since  the  band  model  parameters  extracted  from  tns.  line  tabulation  do  not  predict  the  ob- 
served hot -through -cold  radiances  well  in  the  4.3  pm  region,  they  are  probably  not  useful  for 
predicting  when  high  temperature  correlation  effects  will  be  important.  This  is  most  likely 
the  result  of  omission  of  most  of  the  ’'hot”  lines  from  the  AFCRL  tabulation.  At  lower  tem- 
peratures, where  the  tabulation  should  be  valid,  the  AFCRL  tabulation  predicts  ft  values  at 
only  2 to  3,  even  at  mid-band,  so  that  one  would  indeed  expect  correlation  effects  under  these 
conditions. 

The  short  wavelength  wing  of  the  4.3  pm  C02  band  consists  primarily  of  the  higher  rota- 
tional lines  of  the  v 3 fundamental  band,  the  so-called  "blue  spike.”  One  would  expect  correla- 
tion effects  to  be  important  here  except  right  at  the  band  head,  and  the  measurements  do  in- 
deed show  this. 

Measurements  in  the  2.7  pm  C02  band  show  moderate  line  correlation  effects,  but  the 
poor  agreement  of  both  sets  of  band  model  parameters  with  the  measurements  makes  it 
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difficult  to  assess  the  degree  to  which  these  effects  can  be  handled  using  present  calculation 
techniques.  One  would  expect,  however,  that  the  development  of  % set  of  parameters  which 
adequately  treats  this  band  is  the  primary  missing  link;  the  calculation  techniques  which  work 
for  the  2.7  pm  water  band  and  the  4.3  pm  COs,  band  should  also  be  suitable  for  treating  cor- 
relation effects  in  the  2.7  pm  COg  band. 

To  summarize,  correlation  effects  are  important  in  the  2.7  pm  water  band.  They  are  ob- 
served in. all  the  measurements  throughout  the  band.  Therefore,  special  attention  should  bs 
given  to  plume  calculation  techniques  to  be  sure  that  correlation  effects  are  included,  as  hi 
the  model  of  Section  2.  Correlation  effects  were  observed  only  in  the  "blue  spike"  region  of 
the  4.3  pm  C02  band,  but  parameter  values  indicate  that  correlation  techniques  could  be  im- 
portant throughout  the  band  at  lower  temperatures  and  pressures.  Correlation  effects  were 
observed  throughout  the  2.7  pm  CC>2  band,  but  they  were  less  severe  than  those  observed  in 
water  vapor. 


5.2  CONCLUSIONS  CONCERNING  THE  CONSISTENCY  OP  THE  MEASUREMENTS  WITH 

BAND  MODEL  PARAMETERS 

The  band  model  parameters  of  General  Dynamics  do  reasonably  well  in  calculating  the 
radiance  expected  from  the  isothermal  hot  cell  alone.  The  agreement  for  water  vapor  is  rea- 
sonably good  throughout  the  band. 

For  the  4.3  pm  C02  band,  the  long  wavelength  wing  is  not  as  well  predicted  as  the  re- 
mainder of  the  band;  this  effect  produces  substantial  disagreement  in  the  prediction  of  hot- 
through-cold  radiances  in  this  portion  of  the  band. 

The  band  model  tabulations  for  the  2.7  pm  COg  band  produce  radiance  predictions  sub- 
stantially different  from  those  observed.  These  disagreements  appear  to  be  real  and  substan- 
tial, particularly  when  the  tabulations  are  used  to  compute  radiances  of  hot  targets  viewed 
through  cool  atmospheric  paths.  Therefore,  it  appears  that  a series  of  isothermal  measure- 
ments are  needed  to  resolve  the  dlfficiencies  both  in  the  long  wavelength  wing  of  the  4.3  pm 
band  and  throughout  the  2.7  pm  band.  It  appears  that  the  calculations  made  by  Malkmus,  from 
which  the  existing  General  Dynamics  C02  parameters  were  obtained,  are  quite  consistent  with 
existing  measurements,  but  differ  enough  in  the  details  of  the  long  wavelength  wing  to  produce 
non -negligible  inaccuracies  in  hot -through -cold  radiances. 

The  General  Dynamics  parameters  do  not  predict  the  cold  cell  transmittances  well.  This 
is  probably  the  result  of  the  fact  that  they  were  obtained  primarily  by  extrapolation  from  high 
temperature  measurements. 

The  band  model  parameters  derived  from  the  AFCRL  line  tabulation^  do  an  adequate 
job  of  predicting  the  measurements  only  when  all  lines  which  are  significant  at  the  tempera- 
tures involved  are  included  in  the  tabulation.  This  is  the  case  for  the  cold  cell  transmittances 
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and,  it  appears,  for  regions  in  the  center  of  the  water  vapor  band  at  2.7  pm.  At  the  edges  of 
the  2.7  pm  band,  however,  the  parameters  from  the  AFCRL  tabulation  do  not  predict  the 
observed  absorption  accurately.  That  the  band  center  is  well  treated  while  the  band  edges 
are  poorly  treated  indicates  that  lines  are  missing  mostly  from  the  band  edges.  This,  in 
turn,  probably  indicates  that  important  higher  level  rotational  lines  are  missing.  Fcr  C02, 
the  CRL  line  parameters  agree  well  on  the  short  wavelength  side  but  poorly  on  the  long  wave- 
length side.  The  short  wavelength  wings  contain  the  band  heads  of  all  the  R-branches  of  the 
bands;  if  the  tabulation  extends  to  the  band  heads  or  nearly  so,  this  could  be  the  reason  for 
the  good  agreement  in  the  R-branch.  Missing  rotational  lines  might  then  be  the  reason  for  the 
poor  agreement  in  the  long  wavelength  wing  (P -branches). 

We  have  not  examined  the  line  tabulation  in  order  to  sort  out  the  reasons  for  this  lack  of 
agreement,  except  to  note  that  the  tabulation  of  the  fundamental  extends  out  nearly  to  the 
Vg  fundamental  band  head  which  produces  the  "blue  spike"  at  4.17  pm.  The  tabulation  ends 
only  several  lines  short  of  this  band  head. 

At  its  present  state,  the  AFCRL  line  tabulation  is  surely  not  adequate  to  accurately  pre- 
dict plume  radiation  above  temperatures  of  1200°K.  The  reasons  for  the  disagreement  be- 
tween the  line  tabulation  and  the  measurements  made  in  the  program  need  to  be  determined, 
and  the  important  missing  lines  added  to  the  tabulation  where  possible.  This  task  is  now 
underway,  but  the  substantial  disagreement  between  the  line  tabulation  and  the  measurements 
described  herein  indicates  that  there  is  much  to  be  done. 

Previous  discussion  has  dealt  mostly  with  comparisons  of  the  hot  cell  radiance  measure- 
ments with  band -model  calculation  results  using  different  sets  of  parameters.  Since  measure- 
ments were  made  with  rather  short  high  temperature  optical  paths,  the  comparisons  apply 
more  to  the  first  band  model  parameter  (the  line -strength -to -spacing  ratio)  rather  than  the 
second  (which  governs  the  approach  to  the  square -root  region).  But  the  hot -through -cold 
radiance  measurements,  having  been  made  through  long  cold  paths,  depend  heavily  on  the 
ratio  of  a local  value  of  the  2nd  band  model  parameter  in  the  hot  path  to  its  average  value  up 
to  the  hot  path.  This  dependence  is  demonstrated  in  both  the  Lindquist -Simmons  approxima- 
tion and  the  band  model  developed  in  Section  2 by  the  importance  of  the  variable,  y/y.  Agree- 
ment of  either  of  these  methods  with  the  observed  data  depends  heavily  upon  the  accuracy  of 
the  second  band  model  parameter,  thus  creating  a problem  which  is  illustrated  by  the  wings 
of  the  2.7  pm  water  vapor  band. 

The  extreme  wings  of  this  band  in  the  hot  cell  radiance  are  not  significantly  attenuated 
by  the  cold  cell,  as  shown  by  the  measured  hot -through -cold  radiance  data  and  by  cold  cell 
transmittance  measurements.  Yet  both  the  Lindquist -Simmons  and  the  band-model  of  Sec- 
tion 2 predict  substantial  attenuation  in  these  wings  and  produce  much  poorer  percentage 
agreement  with  measurements  in  this  region  than  in  any  other.  Since  the  hot  cell  radiance  is 
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correctly  predicted,  this  poor  agreement  is  either  due  to  a gross  departure  o!  the  wing  of  the 
2.7  pm  band  from  the  band  model  assumptions  or  to  errors  in  the  second  -band  -model  param- 
eter in  the  hot  cell.  It  appears  that  values  for  the  second  parameter  in  these  wings  are  the 
result  of  an  extrapolation  in  the  General  Dynamics  worlJ^ , perhaps  because  in  their  mea- 
surements the  square  root  region  could  not  be  approached.  This  could  easily  account  for  the 
observed  disagreement.  The  accurate  calculation  of  plume  radiances  through  cool  atmospheric 
paths  thus  requires  accurate  values  for  the  second  band  model  parameter  in  the  hot  path.  It 
is  possible  that  a formulation  could  be  found  which  would  allow  the  determination  of  the  second 
band  model  parameter  in  a hot  path  from  a hot -through -cold  radiance  measurement.  Such  a 
technique  would  be  very  useful  in  eliminating  uncertainties  in  these  parameter  values.  It  is 
much  easier  to  produce  experimentally  a long,  cool  path  having  a short  hot  segment  at  one  end 
than  a long,  hot  isothermal  path. 

5.3  CONCLUSIONS  CONCERNING  CALCULATIONAL  TECHNIQUES 

Several  different  techniques  for  calculating  the  radiance  of  a plume  through  a cool  atmo- 
sphere have  been  discussed.  First,  the  band  model  of  Section  2 has  been  developed  and  a 
method  described  for  its  use.  Second,  the  naive  assumption  that  the  plume  radiance,  calculated 
without  accounting  for  the  atmosphere,  can  be  corrected  for  atmospheric  attenuation  by  mul- 
tiplication by  the  atmospheric  transmittance  has  been  discussed.  Third,  comparisons  have 
been  made  using  the  Lindquist -Simmons  approximation  and  two  other  methods.  The  band 
model  developed  in  Section  2 represents  the  line  arrangement  within  the  plume  and  atmosphere 
more  accurately  than  the  other  two;  and,  because  a simple  approximation  to  the  average  equiv- 
alent width  derivative  was  found,  it  is  simple  to  use.  However,  the  degree  to  which  this  band 
model  agrees  with  the  measurements  is  obscured  by  the  disagreement  between  calculations 
based  on  the  available  band  model  parameters  and  the  measured  radiance  values.  Further- 
more, the  difficulty  of  evaluating  the  expression  for  the  average  equivalent  width  derivative  — 
Equation  (34),  made  it  necessary  to  use  an  approximation  without  adequate  evaluation  of  the 
regions  of  validity  of  the  approximation  or  of  the  error  involved  in  it. 

5.4  RECOMMENDATIONS 

Based  on  the  foregoing  conclusions,  a number  of  recommendations  for  necessary  addi- 
tional work  can  be  made. 

(1)  The  function  described  in  Equation  (34)  must  be  carefully  evaluated  and  appropriate 
approximations  determined,  including  regions  of  validity  and  errors  involved.  The 
departure  from  inverse  square-root  limit  for  high  temperature  targets  should  also 
be  examined. 

16.  C.B.  Ludwig,  '‘Measurements  of  the  Curves -of -Growth  of  Hot  Water  Vapor,"  Applied 

Optics,  Vol.  10,  p.  1057 f (1971). 
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(2)  Similar  band  models  for  other  distributions  of  line  strengths  (exponential,  exponen- 
tially tailed  inverse,  and  uniform)  should  be  evaluated  in  order  to  examine  the  effect 
of  curve  of  growth  on  the  band  model. 

(3)  The  effect  of  path  segment  size  on  the  error  involved  in  making  plume  radiance  cal- 
culations should  be  carefully  investigated  for  this  new  model. 


(4)  From  the  GD  and  the  AFCRL  tabulations,  a new  set  of  band-model  parameters  should 
be  developed  which  will  provide  the  low  temperature  accuracy  of  the  AFCRL  tabula- 
tion and  the  high  temperature  accuracy  of  the  GD  tabulation. 

(5)  Work  should  begin  on  a program  to  determine  a new  set  of  band  model  parameters 
for  the  2.7  pm  C02  band  and  the  long-wavelength  wing  (the  ’’red  spike")  of  the  4.3 
/zm  C02  band.  This  program  should  receive  high  priority  since  none  of  the  available 
parameters  predict  the  2.7  /zm  C02  radiances  accurately  enough  for  useful  model 
comparisons.  The  same  statement  is  true  to  a lesser  extent  in  the  red  spike  region. 
The  apparatus  used  for  the  current  measurements  would  be  ideal  for  such  a program. 

(6)  The  calculation  method  developed  in  Section  2 should  be  extended  to  include  Doppler 
broadening,  and,  following  that,  mixed  Lorentz  and  Doppler  broadening.  This  will 
extend  the  usefulness  of  these  calculation  techniques  to  higher  altitude  cases.  Cor- 
relation effects  are  likely  to  become  more  important  as  the  pressure  is  reduced  so 
that  the  availability  of  calculational  techniques  which  adequately  treat  line  correla- 
tion effects  becomes  particularly  important  as  altitude  is  increased. 


(7)  Only  minimal  line  correlation  effects  were  observed  in  the  4.3  pm  C02  "red  spike" 
region  at  the  conditions  under  which  measurements  were  made  (1200°K,  15  km  simu- 
lated altitude).  Correlation  effects  are  expected  to  become  more  important  as  tem- 
perature is  reduced.  (Many  plumes  of  interest  are  at  lower  temperature.)  Thus, 
hot -through -cold  measurements  should  be  made  at  lower  temperatures  than  1200°K 
to  investigate  the  correlation  effects  expected.  These  measurements  can  be  per- 
formed with  existing  facilities  by  using  a different  grating  with  dispersion  optimized 
for  the  4.3  pm  region  and  by  degrading  the  spectral  resolution  to  an  appropriate 
value  of  10  cm"*.  (The  present  program  utilized  a 4 cm"1  resolution.) 


(8)  The  degree  to  which  the  lines  of  different  gases  are  truly  uncorrelated  remains  to 
be  explored.  Some  measurements  of  these  effects  between  HgO  and  C02  should  be 
made  in  the  2.7  pm  region,  if  for  no  other  reason  than  to  check  for  non -correlation. 
Further,  for  gases  having  widely  spaced  lines,  such  as  HF,  the  assumption  that  lines 
of  different  gases  are  uncorrelated  could  prove  to  be  invalid  for  certain  narrow  spec- 
tral intervals.  If  such  gases  are  important  plume  species,  measurements  or  calcula- 
tions to  assess  these  efft  cts  should  be  undertaken. 
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Appendix  A 

EVALUATION  OF  EQUATION  (34) 

In  order  to  be  able  to  use  the  results  of  the  band  model  development  of  Section  2,  it  is 
necessary  to  evaluate  the  expression  of  Equation  (34),  namely 
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Values  of  g are  needed  for  the  following  ranges  of  variables: 

0.001  < 1000. 

0.2  x 10 "7  < y/y  < 200. 

0.16  < q/q  < 1.0. 

Several  approaches  were  used  in  obtaining  the  approximations  for  g given  in  Table  1. 
Evaluation  by  Series 

Initial  efforts  to  evaluate  this  function  centered  around  a series  approximation'"  to  the 
normalized  gamma  function,  y*.  From  Reference  11, 

y‘(a,x)  = fjsjg(a'^)n.  (53) 
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The  integral  can  now  be  performed  using  the  formula  3.259.3  of  Gradshteyn  and  RyzhuJ17^  to 
yield  the  following  two  equivalent  expressions. 

*We  wish  to  acknowledge  the  contributions  of  B.D.  Spieth  of  AFRPL  to  this  portion  of  the 
evaluation. 


17.  I.S.  Gradshteyn  and  I.W.  Ryzhik,  Table  of  Integrals,  Series  and  Products,  Academic  Press, 
New  York,  New  York,  1965. 
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The  two  expressions  (55)  and  (56)  are  advantageous  for  different  situations.  Since  the 
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hypergeometric  series  in  Equation  (55)  is  a power  series  in  [1  - l/(y/y)  ] it  will  converge 
more  rapidly  than  (56)  when  y/y  > 1.  Conversely,  Equation  (56)  will  converge  more  rapidly 
than  (55)  when  y/y  < 1.  Using  this  rationale,  after  writing  out  and  simplifying  the  appropriate 
series  for  each  y/y  regime,  we  get  the  following.*  For  y/y  < 1: 
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For  y/y  > 1: 
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For  an  isothermal  path  (q/q  = y/y  = 1),  the  integral  reduces  to  the  following: 


*The  notation  of  Gradshteyn  and  RyzhllJ1^  has  been  used,  i.e., 
(2n  - 1).’.’  = 1-3-5-7  2n  - 1 


(2n).'.’  = 2*4‘6‘8  2n  = 2n-nl 
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(59) 


GO 


rfu  1 1^  V ("D02”'1”  (2°  " !)•  •' 

(Hyir 


n=0 


ThL«  series  should  be  identical  to  the  derivative  of  the  curve -of -growth  of  Equation  (39)  since 
the  band  model  is  identical  to  that  leading  to  Equation  (39)  for  an  Isothermal  path.  This  is 
shown  to  be  true  in  the  following  way. 


«(u,U)  = s 


1 dW 
S(a)  ^ 


_ 6 d(W/6)du 
' §{a)  du  da 


From  Equation  (39), 

j = 2jry[e‘uIo(u)  + 2u  e^fyu)  + Ij(u)]  - lj 
and,  from  Equation  (29),  du/da  = S/2iry6  so  that 

g(u,l,l)  = ^|e_UI0(u)  + 2u  e^yu)  + Ij(u)]  - lj 


= e'u[Io(u)  + I1(u)] 


(60) 


Although  we  did  not  show  the  equivalence  of  Equations  (59)  and  (60)  analytically,  we  did 
determine  that  they  were  numerically  identical  for  u < 1. 

It  was  at  this  point  that  difficulties  were  encountered.  As  u - 0,  g - 1 regardless  of  the 
value  of  y/y  or  q/q.  The  value  of  g begins  to  depart  substantially  from  unity  for  u greater 
than  1,  so  that  large  values  of  u are  of  primary  interest.  The  alternating  series  representa- 

O 

tions  developed  above  satisfy  the  convergence  tests  for  all  u.  The  inner  series  in  (y/y)  con- 

verges  easily  since  the  term  in  (y/y)  never  exceeds  unify.  But  the  outer  series  involves 

powers  of  u,  and  although  the  terms  of  this  series  eventually  approach  zero,  they  increase  in 

magnitude  (although  they  alternate  in  sign)  for  the  first  few  terms.  For  u = 10,  these  terms 
20 

become  so  large  (~10  ) before  they  decrease  that  there  is  no  hope  of  evaluating  the  series  on 
the  available  computer  system  without  severe  round-off  error.  Furthermore,  we  are  most 
interested  in  the  asymptotic  behavior  of  g for  large  u.  Thus,  a series  representation  of  the 
type  developed  above  is  not  very  useful  for  present  purposes. 

We  then  sought  a series  representation  for  g valid  for  large  u.  However,  since  the  argu- 
ment of  the  normalized  gamma  function  in  Equation  (34)  takes  on  all  values  between  0 and  u 
during  the  integration,  a series  representation  for  y*  valid  only  for  large  arguments  will  not 
suffice.  More  analytical  effort  is  required  here. 

Numerical  Evaluation 

In  order  to  obtain  a limited  set  of  values  for  the  function  of  Equation  (34),  the  function  was 
evaluated  numerically  after  first  making  the  following  transformation: 
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;.'; ^ /. At  - ,A‘.i  -*  J„  uL!}*  , * , , 1*  ^ *>. 


■/«**■'  •-  *^s’  r*  »t*m^ 


r**  w *,"  n '■r  • 


v = tan  0/2 


(61) 


Then 


«(U,  rff,  q/5 ■ fr/7  ^r(V®  f 

5 COS  0[(y/y)Z  - 1]  + (y/y)2  + 1 


(62) 


This  function  was  evaluated  graphically  after  plotting  the  integrand.  Values  of  y*  were  com- 
puted from  expressions  6.5.29  and  6.5.32  in  Reference  11.  Values  were  computed  for  q/q  = 
1/4  (corresponding  to  the  value  encountered  in  the  measurements)  and  for  y/y  = 1.8  and  2.0. 
These  corresponded  to  the  values  needed  at  3550  cm"*  in  order  to  compare  with  the  measure- 
ments. Values  were  computed  for  u = 1, 10,  100,  and  1000.  The  product  of  q/qr(q/q)  and  the 
integrand  of  Equation  (62)  is  plotted  in  Figure  A-l  for  the  four  u values,  from  which  results 
the  graphical  integrations  were  performed.  Figure  A -2  shows  the  computed  values.  These 
evaluations  provided  some  idea  of  the  behavior  of  the  function. 

Asymptotic  Values 

In  order  to  investigate  the  behavior  of  the  function  it  is  necessary  to  know  its  behavior 
at  various  limits.  At  u = 0,  the  function  is  easily  shown  to  be  unity.  For  the  isothermal 
case,  (y/y  = q/q  = 1),  the  asymptotic  values  can  be  found  from  Equation  (60). 

lim  [g(u,  1,  1)]  = l/~  (63) 

u-oo  ’ 


For  the  Lindquist -Simmons  approximation,  the  limit  is  shown  to  be  (see  Reference  9) 

lim  g (a,  y/y)  - (64) 

u—oo  y T 


Thus  inverse  square  root  behavior  for  large  u might  be  expected  as  all  of  the  lines  approach 
the  square  root  region. 

A similar  asymptotic  limit  for  Equation  (34)  was  obtained  by  operating  on  Equation  (15) 

in  Section  2 (i.e.,  before  the  reduction  in  terms  of  the  incomplete  gamma  function).  A term 

o 

of  unity  was  neglected  in  comparison  to  v in  both  the  denominator  of  the  integrand  and  the 
denominator  of  the  exponential.  Then  the  expression  (y/y2)  - l/v2  was  assumed  to  be  small 
for  all  v so  that 

(y/y)2  - 1 

2 

v J 

We  obtained  the  result: 


1 . Wrt-  i 5 exp 
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* V 
*r 


V u'^du 

i _ i[u' + (y/y)2  - 1] 

llm  g(u,y/y,q/q)  =— y/y  q/q T7i= — 

u-oo  V^r  (2u)q/q 


1/2 


This  could  be  evaluated  only  for  integer  q/q  to  yield 


(65) 


lim  g(u,  y/y,  q/q)  = 
u-oo 


JLt| 1 

if  y q V2u(q/q  - 1/2) 


(66) 


for  large  u.  This  expression  has  a singularity  at  q/q  = 1/2.  Although  the  expression  was 
strictly  valid  for  integer  q/q,  it  appears  to  be  useful  for  non-integer  values  larger  than  1 as 
well  (although  these  values  are  not  of  significant  interest).  The  integral  of  Equation  (65)  was 
evaluated  for  q/q  = 1/2  to  yield 


lim  g(u,  y/y,  1/2)  = ^ y/y  (1/2) 
u--oo 


VT 


(67) 


In  order  to  obtain  corresponding  relations  for  q/q  < 1/2,  it  became  apparent  during  the 
numerical  integration  that  for  almost  all  of  the  range  of  important  values,  as  well  as  for 
small  q/q 


We  introduce  this  to  obtain 


oo 


2\q/q 


g<», y/f.  <i/i)  - U dv 

y q \q/  ^ y/y  + v 


(68) 


The  product  of  the  integrand  of  this  function  and  q/qr(q/q)  is  plotted  in  Figure  A-S  after 
transformation  into  the  variable  <t>,  along  with  the  exact  function  for  u = 1000,  q/q  1/4,  and 
y/y  = 2.  It  can  be  seen  that  the  approximation  agrees  well  with  the  exact  representation  of 
Equation  (34)  except  near  <f>  = it  and  should  not  contribute  substantial  error  to  the  integral. 

This  agreement  is  expected  to  remain  true  for  small  q/q  and  large  u.  It  can  be  seen  from 
Equation  (68)  that,  where  the  approximation  holds,  the  transmittance  derivative  will  follow  a 
u"Vq  dependence.  This  somewhat  surprising  result  indicates  that  square  root  growth  is  never 
fully  achieved  for  some  lines  and  may  be  the  result  of  allowing  the  integration  over  strength 
SQ  (i.e.,  w)  to  extend  completely  to  zero,  whereas  the  inverse  strength  distribution  cannot  be 
normalized  when  this  occurs.  The  integral  in  Equation  (68)  cannot  be  performed  in  closed 
form  except  for  some  special  cases.  For  y/y  = 1,  q/q  < 1/2, 


i 


! 


i 
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to  g<«,wa « |IW5XW V?-  m 

For  (y/y)2  large  enough,  unity  in  the  numerator  of  Equation  (68)  can  be  ignored  compared 
2 

to  v to  yield 

lim  g(u,  y/y,  q/q)  = |r  2u) cosec  rr(q/q  + 1/2)  (70) 

The  integrand  for  this  additional  approximation  is  also  shown  in  Figure  A -3.  The  approxima- 
tion becomes  more  valid  as  u becomes  larger.  The  limit  given  by  Equation  (70)  is  plotted  on 
Figure  A-l  along  with  the  computed  values. 

The  asymptotic  limits  were  compared  in  order  to  obtain  a simple  generalized  form  which 
can  be  useful  for  representation  of  a large  range  of  values  of  q/q,  which  does  not  have  a sin- 
gularity at  q/q  = 1/2.  Figure  A -4  shows  the  asymptotic  limits  of  those  terms  of  Equations  (66), 
(69),  and  (70)  governing  their  behavior  with  q/q.  Thus,  for  Equation  (66)  we  have  plotted 

\^!(q/q -1/2) 

For  Equation  (69)  we  have  plotted 


For  Equation  (70)  we  have  plotted 
cosec  ir(q/q  + 1/2) 

q \q/ 

The  approach  to  the  singularity  at  q/q  = 1/2  is  apparent.  By  examination  of  this  figure,  we 
have  chosen  unity  as  being  representative  of  these  functions  for  q/q  < 1/2  without  the  singu- 
larities, and  V2/ir(q/q)  as  being  representative  for  1/2  < q/q  < 1.  (Values  for  q/q  > 1 are 
not  required.)  These  limits  are  shown  on  Figure  A -4. 

As  a result  of  the  preceding  paragraphs,  a set  of  asymptotic  limits  which  appear  to  be 
suitable  for  most  situations  that  we  expect  to  encounter  has  been  synthesized.  To  summarize: 

lim  g(u,  y/y,  q/q)  = (y/y)2^^”^  for  | < 1/2 

u— co  q 

Approximate  Functional  Form  for  Equation  (34) 

Curves  of  growth  of  the  type  given  in  Equations  (39)  and  (49)  have  been  represented  in 
the  form  (see  Reference  12) 


Since  Equation  (34)  for  an  isothermal  case  can  easily  be  shown  to  follow  such  a law,  we 
attempted  to  fit  such  a function  to  the  points  computed  for  q/q  « 1/4,  y/y  = 2 and  the  asymp- 
totes given  in  Equation  (71).  The  result,  also  shown  in  Figure  A-l,  agreed  quite  well  with  the 
computed  data  points.  On  this  basis,  and  in  view  of  the  necessity  of  approximating  Equa- 
tion (34)  without  the  time  or  funds  to  perform  a more  detailed  examination  of  its  properties, 
this  form  of  expression  was  used  to  represent  the  function  for  all  q/q  < 1 and  all  y/y,  using 
the  asymptotic  limits  given  in  Equation  (71).  This  procedure  leads  directly  to  the  specifica- 
tions of  Equation  (72)  given  in  Table  1 of  Section  2. 

Discussion 

We  have  evaluated  Equation  (34)  exactly  for  only  a few  values  of  u,  q/q,  and  y/y.  Asymp- 
totic limits  for  large  u were  examined  for  most  combinations  of  q/q  and  y/y.  Because  of 
time  and  budget  limitations,  the  bounds  within  which  these  asymptotic  forms  were  valid  were  not 
determined  carefully.  Therefore,  the  rather  simple  approximation  used  may  not  represent 
the  function  well  in  certain  domains.  The  crossovers  of  these  functional  forms  (see  Fig- 
ure 4,  Section  2)  are  indications  of  this.  More  detailed  investigations  of  the  behavior  of  this 
function,  including  bounds  for  the  validity  of  the  approximations  as  well  as  approximations  for 
values  of  u near  zero,  still  need  to  be  investigated.  Since  the  band  model  of  Section  2 is  a 
substantial  improvement  over  previous  band  models,  such  additional  effort  would  be 
worthwhile. 


Appendix  B 

COLLECTION  OF  MEASURED  DATA  ' 

In  order  that  the  data  acquired  in  this  program  be  available  to  others  wishing  to  test  cal- 
culation techniques,  this  appendix  includes  all  of  the  hot  cell  radiances,  the  cold  cell  trans- 
mittances,  and  the  corresponding  cell  radiances  as  measured  through  the  cold  cell.  Curves 
are  also  included  showing  the  product  of  the  hot  cell  radiance  and  the  cold  cell  transmittance. 
The  difference  between  this  latter  function  and  the  measured  hot -through -cold  radiance  shows 
the  degree  to  which  line  correlation  effects  are  important. 


The  organization  of  the  data  set  for  each  test  is  as  follows.  First  a table  is  shown  which 
gives  the  detailed  physical  parameters  corresponding  to  that  test.  Second,  a figure  is  pre- 
sented showing  two  curves,  the  upper  being  the  measured  unattenuated  hot  cell  radiance  and 
the  lower  the  measured  hot  cell  radiance  as  attenuated  by  the  cold  cell.  Third,  a figure  show- 
ing the  corresponding  measured  cold  cell  transmittarce  is  given.  Fourth,  we  present  a figure 
showing  the  hot  cell  radiance  as  measured  through  the  cold  cell  again  (lower  curve  in  every 
case)  and  the  product  of  the  unattenuated  hot  cell  radiance  and  the  cold  cell  transmittance 
(upper  curve).  The  difference  between  the  two  curves  on  this  figure  demonstrates  the  degree 
of  line  correlation  effects  present  in  each  case.  Last,  a plot  is  presented  of  the  measured 
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apparent  temperature  as  a function  of  wavenumber.  The  degree  to  which  these  values  deviate 
from  the  hot  cell  temperature  given  in  the  table  are  a measure  of  the  radiometric  consistency 
of  the  data. 
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TABLE  B-l.  TEST -l  PARAMETERS 


Spectral  Region 

2.7  4 m 

Hot  Cell 

Simulated  Altitude:  16  km 
3.0  atm -cm  C02 

Total  Absolute 
Pressure 

76.0  mm 

p 

H2° 

0.0  mm 

p 

pco2 

38.0  mm 

?N2 

38.0  mm 

L » 0.6  m 

T = 1202°K 

Cold  Cell 

Simulated  200 -km  Path 
(o'  15  -tan  Altitude 

Total  Absolute 
Pressure 

98.0  mm 

V 

0.0  mm 

pco2 

88.63  mm 

PN2 

9.37  mm 

L = 100  m 

T = 295°K 

WAVENUMBER  (CM-D 

FIGURE  B-I(c).  MEASURED  HOT  THROUGH  COLD  RADIANCE  (LOWER)  AND  PRODUCT  OF  HOT  CELL  RADIANCE 

AND  COLD  CELL  TRANSMITTANCE  (UPPER)  - Test  1 


WAVENUMBER  CCM-1) 

FIGURE  B-l(d).  MEASURED  APPARENT  TEMPERATURE  - Test  1 


TABLE  B-2.  TEST -2  PARAMETERS 


Spectral  Region 

2.7  nm 

Hot  Cell 

Simulated  Altitude:  16  km 
3.0  atm -cm  C02 

Total  Absolute 
Pressure 

76.0  mm 

PH2° 

0.0  mm 

pco2 

38.0  mm 

\ 

38.0  mm 

L = 0.6  m 

T = 1202°K 

Cold  Cell 

Simulated  50 -km  Path 
& 15 -km  Altitude 

Total  Absolute 
Pressure 

98.0  mm 

ph2o 

0.0  mm 

pco2 

22.15  mm 

?N2 

75.86  mm 

L = 100  m 

T = 2970K 
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WAVENUMBER  CCM-U 

FIGURE  B-2(c).  MEASURED  HOT  THROUGH  COLD  RADIANCE  (LOWER)  AND  PRODUCT  OF  HOT  CELL 

TRANSMITTANCE  (UPPER)  - Test  2 


FIGURE  B-2(d).  MEASURED  APPARENT  TEMPERATURE  - Test  2 


TABLE  B-3.  TEST -3  PARAMETERS 


Spectral  Region 

2.7  Mm 

Hot  Cell 

Simulated  Altitude:  16  km 
3.0  atm -cm  CO2 

Total  Absolute 
Pressure 

76.0  mm 

V 

0.0  mm 

Pc02 

38.0  mm 

P“j 

38.0  mm 

L = 0.6  m 

T = 1200°K 

Cold  CeU 

Simulated  10 -km  Path 
0 15 -km  Altitude 

Total  Absolute 
Pressure 

98.0  mm 

p 

h2o 

0.0  mm 

PC02 

4.43  mm 

?N2 

93.57  mm 

L ---  100  m 

T = 297°K 

143 


144 


146 


146 


FIGURE  B-3(c).  MEASURED  HOT  THROUGH  COLD  RADIANCE  (LOWER)  AND  PRODUCT  OF  HOT  CELL  RADIANCE 

AND  COLD  CELL  TRANSMITTANCE  (UPPER)  - Test  3 


FIGURE  B-3(d).  MEASURED  APPARENT  TEMPERATURE  - Test  3 


TABLE  B -4.  TEST-5  PARAMETERS 


Spectral  Region 

.4.3  jm 

Hot  CeU 

Simulated  Altitude:  16  km. 

0.598  atm -cm  COg 

Total  Absolute 
Pressure 

76.0  mm 

P"2° 

0.0  mm 

pco2 

7.58  mm 

Pf,2 

68.43  mm 

L = 0.6  m 

T = 1200°K 

Cold  Cell 

Simulated  200 -km  Path 

@ 15 -km  Altitude 

Total  Absolute 
Pressure 

98.0  mm 

p 

h2° 

0.0  mm 

*0 

o 

o 

to 

88.63  mm 

?N2 

9.37  mm 

L = 100  m 

T = 297 °K 
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WAVENUMBER  CCM-1) 

FIGURE  B -4(a).  MEASURED  HOT  CELL  RADIANCE  (UPPER)  AND  HOT -THROUGH -COLD  RADIANCE  (LOWER) 

Test  5 


•I 4— 

or  so* 
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WRVENUMBER  (CM-li 

FIGURE  B-4(c).  MEASURED  HOT  THROUGH  COLD  RADLANCE  (LOWER)  AND  PRODUCT  OF  HOT  CELL  RADIANCE 

AND  COLD  CELL  TRANSMITTANCE  (UPPER)  - Test  5 


TABLE  B-5.  TEST -6  PARAMETERS 


Spectral  Region 
Hot  Cell 


Cold  Cell 


4.3 

Simulated  Altitude:  16  km 
0.598  atm -cm  C02 

Total  Absolute 


Pressure 

76.0  mm 

p . 

0.0  mm 

-h2o 

7.58  mm 

-co2 

?N2 

68.43  mm 

L = 0.6  m 

T = 1200°K 

Simulated  50 -km  Path 
@15 -km  Altitude 


Total  Absolute 


Pressure 

98.0  mm 

p 

0.0  mm 

H2° 

p.  _ 

22.15  mm 

-co2 

?N2 

75.85  mm 

L = 100  m 

T = 297°K 
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WAVENUMBER  CCH-ll 

FIGURE  B-5(a).  MEASURED  HOT  CELL  RADIANCE  (UPPER)  AND  HOT -THROUGH -COLD  RADIANCE  (LOWER) 
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FIGURE  B-5(d).  MEASURED  APPARENT  TEMPERATURE  - Test  6 


TABLE  B-6.  TEST -7  PARAMETERS 


spectral  Region 
Hot  CeU 


Cold  Cell 


4.3  n m 

Simulated  Altitude:  16  km 
0.598  atm -cm  C02 

Total  Absolute 


Pressure 

76.0  mm 

ph2° 

0.0  mm 

pco2 

7.58  mm 

68.43  mm 

L = 0.6  m 

T = 1200°K 

Simulated  10 -km  Path 
@ 15 -km  Altitude 

Total  Absolute 

Pressure 

98.0  mm 

ph2o 

0.0  mm 

pco2 

4.43  mm 

Pn2 

93.57  mm 

L = 100  m 

T = 297°K 
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WAVENUMBER  CCM-U 

FIGURE  B-6(a).  MEASURED  HOT  CELL  RADIANCE  (UPPER)  AND  HOT -THROUGH -COLD  RADIANCE  (LOWER) 
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TABLE  B-7.  TEST-8R  PARAMETERS 


Spectral  Region 
Hot  Cell 


Cold  Cell 


2.7  ftm 

Simulated  Altitude:  16  km 

3.0  atm -cm  HgO 


Total  Absolute 

Pressure  76.7  mm 


L = 0.6m 


38.0  mm 
38.7  mm 
0.0  mm 
T - 1215°K 


Simulated  86.0 -4cm  Path 
& 18 -km  Altitude 


Total  Absolute 
Pressure 


•o 


L = 100  m 


53.3  mm 
2.28  mm 
0.0  mm 

51.0  mm 
T = 296°K 


163 


WAVENUMBER  CCM-i) 

FIGURE  B-7(a).  MEASURED  HOT  CELL  RADIANCE  (UPPER)  AND  HOT -THROUGH -COLD  RADIANCE  (LOWER) 

Test  8R 
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WAVENUMBER  CCM-1) 

FIGURE  B-7(b).  MEASURED  COLD  CELL  TRANSMITTANCE  - Test  8R 


WAVENUMBER  CCM-l) 

FIGURE  B-7(d).  MEASURED  APPARENT  TEMPERATURE  - Test  8R 


TABLE  B-8.  TEST -9R  PARAMETERS 


Spectral  Region 

2.7  n m 

Hot  CeU 

Simulated  Altitude:  16  km 
3.0  atm -cm  HgO 

Total  Absolute 
Pressure 

76.0  mm 

p 

H2° 

38.0  mm 

pco2 

0.0  mm 

38.0  mm 

L = 0.6  m 

T = 1202°K 

Cold  CeU 

Simulated  28.8 -km  Path 
(S'  18 -km  Altitude 

Total  Absolute 
Pressure 

53.2  mm 

ph20 

0.76  mm 

pco2 

0.0  mm 

PN, 

52.46  mm 

L = 100  m 


T = 293°K 
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WAVENUMBER  (CM-1) 

FIGURE  B-8(c).  MEASURED  HOT  THROUGH  COLD  RADIANCE  (LOWER)  AND  PRODUCT  OF  HOT  CELL  RADIANCE 

AND  COLD  CELL  TRANSMITTANCE  (UPPER)  - Test  9R 
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FIGURE  B-8(d).  MEASURED  APPARENT  TEMPERATURE  - Test  9R 


TABLE  B-9.  TEST-1 1MR  PARAMETERS 


Spectral  Region 
Rot  Cell 


Cold  Cell 


2.7  |im 

Simulated  Altitude:  16  km 
3.0  atm -cm  HgO 

Total  Absolute 

Pressure  76,0  mm 


ph20 

38.0  mm 

pco2 

0.0  mm 

Pn2 

38.0  mm 

L = 0.6  m 

T = 1202°K 

Simulated  288 -km  Path 
(o'  18 -km  Altitude 

Total  Absolute 
Pressure 

53.2  mm 

ph20 

7.6  mm 

PC02 

0.0  mm 

Pt,2 

45.6  mm 

L = 100  m 

T = 296°K 
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WAVENUMBER  CCM-1) 

FIGURE  B-9(a).  MEASURED  HOT  CELL  RADIANCE  (UPPER)  AND  HOT -THROUGH -COLL  RADIANCE  (LOWER) 

Test  11MR 


FIGURE  B-9(b).  MEASURED  COLD  CELL  TRANSMITTANCE  - Test  11MR 
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WAVENUMBER  (CM-1) 

FIGURE  B-9(c).  MEASURED  HOT  THROUGH  COLD  RADIANCE  (LOWER)  AND  PRODUCT  OF  HOT  CELL  RADIANCE 

AND  COLD  CELL  TRANSMITTANCE  (UPPER  - Test  11MR 


FIGURE  B-9(d).  MEASURED  APPARENT  TEMPERATURE  - Test  11MR 


TABLE  B-10.  TEST -15  PARAMETERS 


Spectral  Region 
Hot  Cell 


Cold  Cell 


4.3  nm 

Simulated  Amine  Exhaust  Mixture 
© 16 -km  Altitude 

Total  Absolute 

Pressure  76.0  mm 


PH2° 

35.26  mm 

pco2 

7.58  mm 

pco 

1.44  mm 

P«2 

4.64  mm 

?N2 

27.09  mm 

L - 0.6  m 

T = 1201°K 

Simulated  200 -km  Path 

(S'  15 -km  Altitude 

Total  Absolute 
Pressure 

98.0  mm 

\° 

8.64  mm 

pco2 

88.63  mm 

PN, 

0.73  mm 

L = 100  m 


T = 297°K 
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WflVtNUMBER  (CM— 1 ) 

FIGURE  B -10(a).  MEASURED  HOT  CELL  RADIANCE  (UPPER)  AND  HOT -THROUGH -COLD  RADIANCE  (LOWER) 

Test  15 


WflVtNUMBER  (CM-1) 

FIGURE  B -10(c).  MEASURED  HOT  THROUGH  COLD  RADIANCE  (LOWER)  AND  PRODUCT  OF  HOT  CELL  RADIANCE 

AND  COLD  CELL  TRANSMITTANCE  (UPPER)  - Test  15 


FIGURE  B -10(d).  MEASURED  APPARENT  TEMPERATURE  - Test  15 


TABLE  B-ll.  TEST -16  PARAMETERS 


Spectral  Region 
Hot  Cell 


Cold  Cell 


4.3  11m 

Simulated  Amine  Exhaust  Mixture 
(S'  16 -km  Altitude 

Total  Absolute 


Pressure 

76.0  mm 

ph20 

35.26  mm 

pco2 

7.58  mm 

p 

*co 

1.47  mm 

PH 

H2 

4.64  mm 

z 

27.06  mm 

L = 0.6  m 

T = 1200°K 

Simulated  50 -km  Path 

(S'  15 -km  Altitude 

Total  Absolute 
Pressure 

98.0  mm 

ph20 

2.17  mm 

pco2 

22.15  mm 

Pn2 

73.68  mm 

L = 100  m 

T = 298°K 
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WAVENUMBER  (CM-t) 

FIGURE  B-ll(c).  MEASURED  HOT  THROUGH  COLD  RADIANCE  (LOWER)  AND  PRODUCT  OF  HOT  CELL  RADIANCE 

AND  COLD  CELL  TRANSMITTANCE  (UPPER)  - Test  16 
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WAVENUMBER  (CM~l) 

FIGURE  B -11(d).  MEASURED  APPARENT  TEMPERATURE  - Test  16 


TABLE  B-12.  TEST -17  PARAMETERS 


Spectral  Region 
Hot  Cell 


Cold  Cell 


4.3  4m 

Simulated  Amine  Exhaust  Mixture 
@ 16 -km  Altitude 

Total  Absolute 


Pressure 

76.0  mm 

ph20 

35.26  mm 

PC°2 

7.59  mm 

p 

*co 

1.47  mm 

PH 

"2 

4.67  ram 

PN 

W2 

27.21  mm 

L - 0.6  mm 

T = 1202°K 

Simulated  10 -km  Path 

15 -km  Altitude 

Total  Absolute 

Pressure 

98.0  mm 

ph2o 

0.43  mm 

pco2 

4.43  mm 

?N2 

93.14  mm 

L = 100  m 

T = 296°K 
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WAVENUMBER  tCM-1) 

FIGURE  B -12(a).  MEASURED  HOT  CELL  RADIANCE  (UPPER)  AND  HOT -THROUGH -COLD  RADIANCE  (LOWER) 

Teat  17 


190 


o 


33NUIQUU  THU103dS 


191 


102 


FIGURE  B -1 2(d).  MEASURED  APPARENT  TEMPERATURE  - Test  17 
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TABLE  B-13.  TEST-18  PARAMETERS 
Spectral  Region  2.7  (im 


Hot  Cell 


Cold  Cell 


Simulated  Amine  Exhaust  Mixture 
& 16-km  Altitude 


Total  Absolute 
Pressure 


76.0  mm 


p 

h2o 

35.26  mm 

p 

pco2 

7.58  mm 

pco 

1.47  mm 

?H2 

4.62  mm 

p 

M 

W2 

27.06  mm 

L = 0.6  mm 

T = 1202°K 

Simulated  200  -km  Path 
(?'  15 -km  Altitude 

Total  Absolute 
Pressure 

98.0  mm 

ph20 

8.64  mm 

pco2 

88.63  mm 

\ 

0.73  mm 

L = 100  m 

T = 301°K 
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FIGURE  B-13(b).  MEASURED  COLD  CELL  TRANSMITTANCE  - Test  18 


TABLE  B-14.  TEST-19  PARAMETERS 


Spectral  Region 
Hot  Cell 


Cold  Cell 


2.7  pm 

Simulated  Amine  Exhaust  Mixture 
0 16 -km  Altitude 

Total  Absolute 


Pressure 

76.0  mm 

ph20 

35.26  mm 

pco2 

7.58  mm 

p 

CO 

1.47  mm 

P»2 

4.64  mm 

27.06  mm 

L = 0.6  m 

T = 1202°K 

Simulated  50 -km  Path 

@ 15 -km  Altitude 

Total  Absolute 


Pressure 

98.0  mm 

ph2o 

2.17  mm 

pco2 

22.17  mm 

\ 

73.65  mm 

L = 100  m 

T = 301°K 
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WAVENUMBER  (CM-1) 

FIGURE  B-14(a).  MEASURED  HOT  CELL  RADIANCE  (UPPER)  AND  HOT -THROUGH -COLD  RADIANCE  (LOWER) 

Test  19 


WAVENUMBER  CCM-1) 

FIGURE  B -14(b).  MEASURED  COLD  CELL  TRANSMITTANCE  - Test  19 


TABLE  B-15.  TEST-20  PARAMETERS 


Spectral  Region 
Hot  Cell 


Cold  Cell 


2.7  \i  m 

Simulated  Amine  Exhaust  Mixture 
& 16 -km  Altitude 


Total  Absolute 


Pressure 

76.0  mm 

PH2° 

35.26  mm 

pco2 

7.58  mm 

p 

CO 

1.47  mm 

Pn2 

27.06  mm 

P«2 

4.64  mm 

L = 0.6  m 

T = 1202°K 
Simulated  10 -km  Path 

& 15 -km  Altitude 


Total  Absolute 


Pressure 

98.0  mm 

p 

H2° 

0.43  mm 

pco2 

4.43  mm 

PN2 

93.14  mm 

L = 100  m 

T = 297°K 
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WAVENUMBER  CCM-1) 

FIGURE  B-15(a).  MEASURED  HOT  CELL  RADIANCE  (UPPER)  AND  HOT -THROUGH -COLD  RADIANCE  (LOWER) 

Test  20 
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